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INTRODUCTION 


The problem of determining differential equations by information on solutions 
of these equations is called the inverse problem for differential equations. 


Examples of these kind of problems are the well-known Sturm-Liouville 
problems. Many important application problems connected to elastic 
displacement, electromagnetic oscillation, diffuseand other processes in nature 
and society lead to inverse problems. The extent ofthese problems is expanding 
constantly. At present a substantial amount of literature is devoted to these 
problems. One of the first investigations into these kind of questions was the 
inversion of kinematic problems of seismics, the essence of which consists in 
the determination ofthe velocity of propagation of elastic waves by time of their 
movement. A one-dimensional case of one of such problems was considered by 
Herglots (1905). He obtained the formula of inversion, which later became the 
basis of the solution of many important problems in geophysics, in particular the 
basis of structure determination of the Earth’s crust and the Earth’s mantle. 


As a second classic direction in the theory of inverse problems one can mention 
the inverse problem of the theory of potential, which consists in a form of 
description of body shape and density of this body on a known potential. The 
uniqueness theorem of the solution of one of these problems was proven by 
Novikov (1938) for the first time. 


Problems linked to the Sturm-Liouville equation and its generalization are a 
third direction in the theory of inverse problems. The sense of these problems 
is the following: we know the spectral function or scattering data ofa differential 
operator; it is required to define this operator. The first uniqueness theorems of 
the solution were obtained in papers by Ambarzumjan (1929) and Borg (1945). 


Inverse problems are usually nonlinear and are separated into one-dimensional 
and multidimensional problems, depending on whether the sought function (or 
functions) is a function of one variable or of many. These problems, especially 
multidimensional problems, are often ill-posed in the classic sense. In this 
sense, questions of the uniqueness of the solution and a search of minimal 
information, which makes the inverse problem determined, have particular 
actuality. The theorem of uniqueness of the solution of the complex 
multidimensional inverse problem for the Schródinger equation in the class of 
piecewise-analytical functions was established forthe first time by Berezanskij 
(1958). 


Multidimensionality of inverse problems has particluar value at present, 
because practice shows that many investigating processes are described by an 
equation, of which the coefficient essentially depends on many variables. 
Inverse problems often bring integrals to first kind operator equations. Some 
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inverse problems for hyperbolic equations are, for example, reduced to the 
investigation of integral equations of Volterra type ofthe first kind. In turn, this 
sometimes allows one (chiefly in one-dimensional inverse problems) to find 
equations of the second kind with operators possessing sufficiently good 
properties, for example compact operators, which gives a method of investigation. 
The basis of such a convergence is often a formula for the solution of direct or 
inverse problems. In many cases, particularly when information on the solution 
ofan equation is given by only a part ofthe boundary ofthe considered domain 
(practice is needed in those problems), such a convergence of the inverse 
problem to an integral equation of the second kind often appears to be 
impossible. One ofthe reasons for this is the ill-posedness of these problems. 
These questions require new approaches. The general theory of operator 
equations of the first kind and their applications was developed in papers by 
Tikhonov (1943, 1963), Tikhonov and Arsenin (1977), Lavrent'ev (1955, 
1959) and Ivanov (1963). The new methods developed there, have found wide 
implementations, in particular in the theory and practice of ill-posed problems. 
Many inverse problems are closely connected to problems of integral geometry. 
Thereby, it appears necessary to investigate new problems ofintegral geometry, 
when manifolds in which the sought function (or functions) is integrated, are 
complicated by their structure. Important results and applications in the case of 
linear and other manifolds are obtained from papers by Radon (1917), Courant 
and Hilbert (1962), John (1955), Khachaturov (1954), Kostelyanec and 
Reshetnyak (1954), Gel'fand.(1960), Gel'fand et al. (1966), Helgason (1959), 
Semyanistyi (1966) and others. At present these results on integral geometry 
and other applications have found applications in inverse problems, in particular 
in tomography. Multidimensional problems of integral geometry in the case of 
complex manifolds andtheir connectionto inverse problems are formulated and 
investigated in articles by Lavrent'ev et al. (1970, 1986), Romanov (1987), 
Anikonov (1987b), Anikonov and Pestov (19902), and Bukhgeim (1983, 1986, 
1988). A significant contribution to the theory of inverse problems has also been 
made by Prilepko and Kostin (1993), Isakov (1990), Blagoveshchenskii (1986), 
Anikonov (1975), Anikonov et al. (1993), Lorenzi (1992), Yamamoto (1990, 
1992), Anger (1990), Sabatier (1990a,b) and others. 


This monograph is devoted to statements of multidimensional inverse problems, 
in particular to methods of their investigation. Questions of the uniqueness of 
solution, solvability and stability are studied. Methods to construct a solution 
are given, and in certain cases inversion formulas are given as well. Concrete 
applications of the theory developed here are also given. Where possible, we 
have stopped to consider the method of investigation of the problems, thereby 
sometimes losing generality and quantity of the problems, which can be 
examined by such a method. 


This work was supported by the Russian Foundation for Fundamental Research, 
grant No. 93-011-181. 


CHAPTER 1 


Operator Equations and Inverse Problems 


1.1 DEFINITION OF QUASIMONOTONICITY, THE UNIQUENESS 
THEOREM 


The investigation of inverse problems for differential equations frequently reduces to the 
investigation of the first type operator equations. Such operators may come up as func- 
tionals of the solutions of some differential equation correspond to this equation. The 
other words, the operator is defined as the set of equations and its values are functionals 
of the solutions of these equations. It is often found that the equation is defined simply 
by one real function belonging to some class. So the operator of the inverse problem is 
defined at the set of the functions. The operators of inverse problems may be of complex 
nature; they are nonlinear as a rule. Some of them have the property that if one of the 
functions is larger than another in a subset of the domain of definition of these functions, 
then images of these functions are different. 

It was this property which was used for the proof of the solution uniqueness of one- 
dimensional inverse problems of electroprospect in the paper of (Tikhonov, 1949). In 
the multidimensional case (Berezanskii, 1958) applied this property for the proof of the 
solution uniqueness of the inverse problem of the Schrodinger equation in the class of the 
piecewise analytical functions. We pick out the operators by such a definition. 

Let E be a set of elements z, and let (A)g be some set of real functions A(z) in 
E. Suppose that an element 7 belonging to some set (7) corresponds to every function 
A(z) € (A]g by r = MA, where M is some operator. 

As a covering of set E we understand totality (w) of sets w € E, union of which is 
equal to E. 


Definition. Operator M is called quasimonotonic with respect to the covering {w} of set 
E, if it follows from inequality A,(z) > A2(z) holding for all x belonging to, at least, one 
non-empty set w € {w} that MA, Z M), Ai(z) € (4), i=1,2. 


Lemma 1.1. If for every To € (7) the equation MA = To has a unique solution, then the 
operator M is quasimonotonic with respect to any covering {w} of set E. 


This statement is obvious. 
Our main purpose is to show that in special cases it follows from being quasimonotonic 
with respect to a fixed covering that the solution of the equation MA = 7 is unique. 
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Let R”+! be an (n + 1)— dimensional real Euclidean space for (z, y), zx € R^, —oo < 
y < oo, n 2 0, let E and E be semispaces with y > 0 and y > 0. Denote by (A) the 
set of all real infinite differentiable functions A(z, y), (z, y) € E to be quasianalytical in y 
in this domain. Under the covering {w} of the semispace E we understand the set of all 
hemiballs w(£, p): 


w(é,p) = ((z,y):Ic E? +y? «p, y»0), EER”, p>. 


Lemma 1.2. For every function A(z,y) € (A) , which is not identically equal to zero, 
there exists a hemiball w € {w} and constants c > 0, a > 0 such that the inequality 


IMz,y)| > ey*, (z,y) € {w} 
holds. 


PROOF. Since the function A(z, y) is not identically equal to zero, and quasianalytical in 
y in E, the equality 


ool 
dy* y=0 
is impossible for all k, k = 0,1,2,.... So there exists a smallest number m > 0 and a 
point é € R” such that 
0* d(x, 0) n  9"X(£,0) 
—X e k=0,1,2,...,.m-1, zcR^, ue e (1.1) 


If m = 0,then the statement is obvious. P : 
Let m > 0. Every infinite differentiable function À(z, y) such that A(z,0) = 0 (see 
(Malgrange, 1966)) allows a representation A(z, y) = yh(z, y), where 


1 


haw) = [oa (1.2) 


[uU 


Combining (1.1) and (1.2) we have 


A(z, y) = y"! h(z, y), (1.3) 


where h(z, y) is an infinite differentiable function and h(£,0) #0. Let R(£, p) be a closed 
ball in R?*! of radius p > 0, centered at (£,0) such that h(z,y) £0, (z,y) € R(£,p). By 
continuity h(z,y) and À(£,0) # 0 such ball exists. 


Letc= min |A(z,y)|, and w = R(£,p) N E. By (1.3) we have 
(z,y)ER(E,p) 


lA(z, y)] = y" |h(z, y) 2 cy", (2,4) Ew. 


The Lemma is proved. 
By this Lemma it follows 


Theorem 1.1. Let M be a quasimonotonic operator with respect to the covering {w} 
of semispace E by hemiballs w(f,p). If Mà, = Mh, Ai(z,y) € {A} C {A}z, then 
A(z, y) = A2(z, y) in E. 
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PRooF. If Ai(z, y) Æ A2(z,y), then by Lemma 1.2 there exists a hemiball w € {w} such 

that |Ai(z,y) — A2(z,y)| 2 cy? > 0, (z, y) € o. In particular it follows from this that in 

the domain w either À, > Az or Àz > A; holds. Since the operator M is quasimonotonic, 

MA, # MA, which contradicts the condition of the theorem. The theorem is proved. 
Below we will use either Theorem 1.1, or Lemma 1.2 and its variations. 


1.2 INVERSE PROBLEMS FOR HYPERBOLIC EQUATIONS 


First we will illustrate the application of Theorem 1.1 to the investigation of the uniqueness 
of the solution of the inverse problem for the linear hyperbolic equation, and then we 
proceed to more general cases. 

Let R? be an Euclidean space of (z, y), z € R?, —oo « y < oo. Consider the problem: 
in the semispace y > 0 it is required to find a strictly positive function A(z, y) C C'(R?), 
even in y, if 


1. In the domain z € R?, y € R!, 0 € t < to(z,y), to > 0 there exists a unique twice 
differentiable solution u(z, y, t) of the Cauchy problem 
a 


= = Au + A(z, y)u + f(z, y,t), ulizo 


Qu 
at 0, 


at t=0 
where A is the Laplace operator in (z, y). 


2. In the domain y = 0, 0 € t € go(z) the function r(z,t) = u|,-o: 0 € t € go, 
go(z) > 0 is given. The inverse problem is reduced to the investigation of the 
operator equation MA = r(z,t), A € C'(E), E = ((z,y),y > 0). As covering {w} 
of the semispace y > 0 we take as above the set of all hemiballs w(£, p): w(£,p) = 
{(z,y): lz- EP +y? «p, y>O}, EER, p»0. 


Let o7 (£, p) = ((z,y) : |z -tP +y?<p*, y«0) 


Theorem 1.2. If a function f(z,y,t) is continuous and f > 0 , then operator M is 
quasimonotonic with respect to covering {w} of domain y > 0. 


PROOF. Let Ai(z, y) > 2(z,y) > 0, (x,y) € «wo = w(£o, po). Denote by ux(z, y, t), 
u2(z,y,t) the solutions of the Cauchy problem corresponding to the functions Aj(z, y) 
and As(z,y). Let À = Ay — 32, ü = u — us. In accordance with the Kirchhoff formula we 


have . 
: 1 fA(gü(gt-r) 1 f A(q)ux(q;t — r) 
EFE s nir S N (1.4) 


where 
2 


r? = J (zi - qi +(y—493)", dq = dqidgadgs. 


i=l 
Note that, if 


? 


= lf bgu(gt-r) 
u(z,y,t) = F(a,y,t) + dx i - dq 


rst 
and Fl... > 0, b > 0, then (in any case for sufficiently small t > 0) the inequality 
u(z,y,t) > 0 holds. This follows, for example, from the proof of the existence of the 
solution u(z, y,t) by the method of successive approximations. 
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By the condition of the problem and prepositions we have f(z, y, t) > 0, A1 — A; > 0, 
(z,y) € wo, AX > 0. If Ai(z, y) = A(z, —9) then M — à > 0, (z,y) € wo . Let t be 
such that the points q and (z, y) in (1.4) belong to &^ U wo. From this using the note 
above we conclude that by the strict positivity of the second term in the right part of 
(1.4) for sufficiently small t > 0 we have &(z, y,t) > 0, (z,y) € w Uw. If y = 0 in the 
last inequality, we obtain à(z,0,t) > 0. This means operator M is quasimonotonic. The 
theorem is proved. 

Let, as above, R? be a three-dimensional Euclidean space for (x,y), £ € R?, —oo < 


ð? 2 8? 
y < œ and A = + is the Laplace operator. 
Oy? { Ox ax? 


Denote by {a} a set ‘of infinite differentiable functions a(z, y) to be even and quasian- 
alytical in y. Let A(z, y) and p(z, y) be some functions belonging to set {a}, and u(z, y, t) 
a solution of the Cauchy problem, A > 0, p > 0, 


1 Oru 


; ie of = Au + F(z,y,t, n(z,y),u) + f(z,y,¢), 


»; (1.5) 
ulizo = at UMS 
t=0 


where F and f are fixed infinite differentiable functions. It is well known that a solution 
u(x, y, t) of such a problem for the more general case exists and is unique in some domain 
of (z, y, t). 

Here is considered the inverse problem: it is required to find one of the functions either 
A(z, y) or n(z, y), if in the domain |z| < r, 0 € t < go the function y(z,y) = ul,_o is 
known, where u(z, y, t) is the solution of the problem (1.5). Here the inverse problem is 
also reduced to the investigation of the equations MA = e(z, y), Mp = e(z, y). 


Theorem 1.3. Let function F(z, y,t, p, u) increase strictly monotonic in variables p, u 
in the domain p > 0, u > 0 and F(z,y,t,0,0) = 0, and the function f(z,y,t) > 0. Then 
the inverse problem has no more than one solution in the class (a). 


PROOF. For the proof of the theorem it suffices to verify the operators MÀ and My being 
quasimonotonic. At first we note, that for every A(z, y) € (a) and u(z, y) € (a) solution 
u(z, y, t) of problem (1.5) for sufficiently small t must satisfy the inequality u(z, y, t) > 0, 
i» Q. 

Really, in accordance to Sobolev's formula (see (Sobolev, 1963)) we have 


1 1 1 
=>- J [flodv + — J [F]edv + 7— J [ujAcdo. 


(Here we use the meanings of (Sobolev, 1963)) 
Since o and Ac have an order of increasing no more than 1/7 and o > 0, it is clear by the 
formula, that the sign of u is defined by the first term of the right part of this formula. 
In particular, if f(z, y, t) > 0, then u(z, y, t) > 0. 

Denote as wt and w^ open hemiballs of radius p > 0 centered at point € € R?, |£| < r, 


w*(£ p) = {(2,y),y > 0: |z — €^ 4 y? <p}, 


w (E, p) = ((z, yy < 0 : jæ — e +y? < p’}. 
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Let p(z, y) and p2(z, y) be two solutions of the problem (A(z, y) is fixed) and let pı (z, y) > 
po(z,y), (z, y) Ew = wt Uw. Denote by ui(z, y, t), ux(z, y, t) solutions of the Cauchy 
problem 


i . oe = Au; + F(z,y,t pi ui) +f, 121,2, 
tilio = a T = 0, 
corresponding to the functions pı(x,y) and u;(z, y). Let w = ui — ug. Then we have 
i : ow = Aw + F(z,y,t, py, w + U2) — F(z, y,t, ua, uz), 


Since pu (xz, y) > #2(z, y), (z, y) € w and as above u2(z, y, t) > 0, by F being monotonic the 
inequality F(z, y,t, ui, u3) > F(z,y,t, ua, ua) holds. So by using, for example, Sobolev’s 
formula, we obtain that (x,y) € ©, w(z,y,t) > 0 at least for small t. Setting y = 0 in 
the last inequality we have v(x, t) — y2(z,t) = w(z,0,t) > 0, this is in conflict with the 
problem condition. 

Let now p(z,y) be fixed and let A1(z, y), A2(z,y) be two different solutions of the 
inverse problem. As above, let w = uj — uz and Ai(z, y) > As(z, y), (z, y) € v. Now we 
have 


1 ^w 1 1 Y 8u 
M y ET = Aw + F(z,y,t, i, w + u2) = F(z,y,t, ua, U2) + (5; = x) n 
Ow 
wo = =| =0. 

S ôt t=0 


2 
By the condition of the theorem f(z,y,t) > 0, as is aS 


Therefore, by inequality 
Ai(z,y) > A2(z,y), (2, y) € w we have 


1 1 Or ur 
(5--x) aero (z,y) Ew, 

which as above is in conflict with w(z,0,t) > 0. By Theorem 1.1 and the obtained results 

the statement of the theorem follows. For reduction to other problems see (Lavrent'ev, 

1974). 


1.9 MULTIDIMENSIONAL INVERSE KINEMATIC PROBLEM 
OF SEISMICS 


Let R”+! be the (n + 1)— dimensional Euclidean space for (z, y),  € R^, —oo < y < oo, 
n — radius vector of the sphere N = {z, |z| = 1). 

Consider the problem: in semispace y > 0 it is required to find an infinite differentiable 
function A(z, y) € {A}, 0 < A < M < oo such that: 


1. The shortest line 4(a,b,n) C E exists of metrics ds? = \?(z,y)(dz? + dy?) joining 
points an, bn, —-oo < a < oo, —co < b < oo, n EN, and 
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2. In the domain |a| < oo, |b| < oo, n € Q the function 


T(a,b,n) = J A(z, y) dz? + dy? = MA 


+(a,b,n) 
is given. 
The problem of determining the function A(z, y) by r(a, b, n) is called the inverse kinematic 
problem of seismics, where A(z,y) is interpretated as the inverse of the velocity À = 
1/v(z, y) and r(a,b, n) is the time of moving of the perturbation by y(a, b, n). 
Let " 
atts) scm Uo) 


Seer hend] 


In R” we consider metric ds? = A2(z)dz? = A2(z) |dz|". Denote by 4? (a, b, n) the shortest 
line of metric dso joining points an, bn. Denote by 


1?(a, b,n) = J Jo ldz]. 


^! (a,b,n) 


Theorem 1.4. If r(a,b, n) 4 7°(a, b,n) for every a, b and n, then the inverse kinematic 
problem of seismic has no more than one quasianalitical solution in y: A(z,y) € {A}. 


PROOF. Show that under condition of the theorem the operator 


MA = J A(z, y)ydz? +dy?, A(z,y) E (A]g 


y(a,b,n)} 

is quasimonotonic with respect to the covering {w} of domain E = ((z,y) : y > 0). In the 
same way as Theorem 1.1 the statement will be proved. Let wo = w(£o, t) be a hemiball of 
radius £ > 0 centered at ĉo € R” belonging to the semispace, and let A1(z, y) > As(z, y), 
(z,y) € wo. Denote by +; the shortest line of metric ds? = A2(dz? + dy?), such that 
^i C wo. 

By the condition of the theorem the shortest line y; cannot belong to the hyperplane 
y — 0 totally, in contrast 


Ti = J o oae + dy? = Jz.) |dz| = TP. 
Ji Ui 


Using the inequality A1(z, y) > Az(z, y), (x,y) € wo and the last remark, we have 


n= EOD > [rley ar + ay? > [role faa ay? zi 
n n oe 


which means operator M is quasimonotonic. The theorem is proved. 
More general results are included in the papers by (Anikonov, 1971; 1975). 


REMARK . In the plane case, instead of the condition of Theorem 1.4 it may be taken as 
T(a,b) < r(a, b) + r(b, b), a < bı < b. By the fulfilment of this inequality every interior 
point on the shortest line *(a, 5b) of metrics ds? = A?(z, y)(dz? + dy?) cannot belong to 
the straight y = 0. 

For n > 2 variables a, b can be changed in the domain |a| < ao, |b] < bo, where ao > 0, 
bo > 0 are constant. 
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1.4 ON THE UNIQUENESS OF THE SOLUTION OF THE FREDHOLM 
AND VOLTERRA FIRST KIND INTEGRAL EQUATIONS 


As we noted above, inverse problems for differential equations reduce to the first kind 
operator equations, and often to the first kind integral equations. Let us consider multi- 
dimensional integral equations and prove the uniqueness theorems. 

Let R”+! be the Euclidean space for (z,y), £ € R^, -oo < y < coo, n 2 0 , let V be 
a compact metric set of v. Denote by w(f,t), € € R^, t > 0, (n + 1)— dimensional open 
set, with properties: 


1. w(ĉ,t) is contained in the semispace y > 0 and unicity {w} is the covering of this 
space. 


2. w(£,t) depends continuously of € and t, and limw(€, t) = (£,0). 
An example of the set w(£,t) may be a hemiball of radius t > 0 centered at point £ c R” 
w(é,t)={(2,y): I6" +y? «P, andy 0]. 


Let M(z,y,£, t, v) be a real continuous function on R?*! x R”+! x V possessing prop- 
erties 


1. M(z,y,€,t,v) > 0, (x,y) € w(€, 2), 
2. M(z,y,€,t,v) = 0, (z, y) € R^*' Va(£, t), 


where &(€,t) is closure of w(E,t). 
Let us consider this equation with respect to the function A(z, y) € {A}: 


MA = inf ‘| M(z,y,6,t,v)A(z, y)dzi ...dz,dy = w(£,t). (1.6) 
w(E,t) 


Note, that if function M does not depend on v, then the considered equation is the multidi- 


t 
mensional analog of Volterra equation. In particular, if n = 0 then MA = f M(y,t)A(y)dy. 
o 
t 
Equation y(t) = f M(y,t)A(y)dy is the first kind Volterra equation. 
0 


Theorem 1.5. If MA, = MA, à; € {A}, i = 1,2, then A(z, y) = A2(z,y), £ € R^, 
0<y<oo. 


PROOF. It is sufficiently to state operator M being quasimonotonic, defined by (1.6), 
AEC. 

As covering {w} of the semispace y > 0 we understand sets w(¢, t), satisfy conditions 
1 and 2. 

Let Ai(z,y) > A2(z,y), (z; y) € wo = w(bo, to), 4; € {A}, 1 = 1,2. Since the kernel 
M(z,y,é,t,v) is continuous and has properties 1 and 2, functions 


w;(£,t,v) = f! Mí(z,y,6,t,v)A(r,y)dzdy, i—1,2, 
w(E,t) 
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are continuous in the domain £ € R^, t > 0, v € V. Thus as V is a compact metric space 
there exist elements v; € V so that 


wi(£, t) T inf w(£, t, v) — wil, t, vi), i= 1,2. 


By the above we have 


w1(o, to) = J M6 esto vı) (z, y)dzdy > [| MG etos), y)dzdy 


IV 


J Mie, y, 05 to, v2) da(2, y)dedy = ws (6o, to) 


which means operator M is quasimonotonic. The theorem is proved. 
Suppose now that sets w(f,t) have additional properties: 


1. For every h > 0 set w(,t,h) = w(E,t)N{y >h, zc R^), t» A, is an (n 4 1)- 
dimensional open set, 


2. w(,t, h) depends continuously of £,t, h and limw(€, t, h) = (£, h). 


Define a class {A}. of A(z, y), given by the definition A(z, y) € {A}.o, there exists some 
sequence of numbers h;, 0 = ho € hy € ... € hm..., and some sequence of functions 
A(z,y), hi € y € Aigi, x € R^, analytical in y, which are infinite differentiable with 
respect to (x,y) such that A(z,y) = A;(z,y), x € R”, h; € y € hi41. In other words, class 
(Alo. consists of functions that are piecewise analytical in y and infinite differentiable 
with respect to (z, y). 


Theorem 1.6. If MA, = MX, A; € {A}oo, i = 1,2, then Ai(z, y) = A2(z,y), x € R^, 
y 20. 


PROOF. Let Aj(z, y) and A2(z,y) A; € {A}oo be two solutions of equation (1.6), such as 
MA, = MA;. By continuity of the kernel M(z, y, £, t, v) and conditions 1 and 2, functions 


wi(€, t, v) = f M(zx,y,£,t,v)à:(z, y)dzdy 
w(E,t) 
are continuous. 
Therefore v is a compact metric space, from this in particular it follows the existence 


of v; so that 
wi(€, t) = wil€, t, vi) = inf w;lé, t, v). 


Suppose that Aj(z,y) # As(z,y). From this and by the definition of class {A}. it 
follows the existence of the interval [a,b], 0 < a < b < oo so that the following conditions 
are true: 


1. A(z, y) = Ax(2, y), TE R”, y<a, 
2. à;(xz,y) are analytical in y when a € y < b, and (z, y) Z Ao(z,y), zv ER", a< y<b. 


Let w'(é,t) be the intersection of w(é,¢) with the semispace y > 0. By Lemma 1.2 and 
proposition Aj(z, y) Æ A2(z, y) there exists a set w} = w! (ĉo, to) = w(£o, to) N {y > a} so 
that 


Operator Equations and Inverse Problems 11 


1. Set w} belongs to the strip a < y < b, z € R^, 
2. |Ai(z, y) a A2(z, y)l 2 0, (z, y) € ug. 
Let Àj > Az, (z, y) € ud. Taken this into consideration and by Ay = \2, z € R^, y < a we 
have. 
w1(£o, to) = J M(z,y, ĉo, t, v1)Ardady > J M(z, Y, ĉo, t, v1) Agdady 2 w»(£o, to). 
w(£o,to) w(£oto) 


Thus the existence of the point is shown (£o, to) with wi(£o, to) > w2(ĉo, to). The last 
inequality is in conflict with the condition of the theorem. The theorem is proved. 


Example of nonuniqueness of a solution of the Volterra integral equation 
with a positive kernel. 


Let A(t) # 0, t > 0, A(0) = 0 be an infinite differentiable function with zero in every 
interval (0, a), œ > 0, as for example A(t) = exp (—1/t?) sin (1/t), (0) = 0. Show that 
there exists a kernel K(z,t) > 0, z £ 0, t Z 0 such that for every z > 0 the equality 


J K(z,t)A(t)dt = 0 


holds. 
Let E*(z), E- (x) be sets such that 


E*(z) = (t: A(t) > 0) n (0,2), 

E(x) = (t: Xt) < 0} n (0,2) 

Set 
| Hansda kreeg "sts 


{ -1, if t € E- (z), 


F()-1 o if tg E-(s), 


Introduce functions 
gt(e)= f xod, &()- f xoa. 
E+(z) E-(z) 
Define kernel K(z,t) by K(z,t) = —u*(z,t)g" (zx) — u^ (z,t)g*(z). It is obvious, that 
K(z,t) > 0, z #0. We have 
[ 60x04. f Kle dAd f K(, (08 = 9*(0)97(2)-9* (2) (a) = 0. 
g Et(z) E-(z) 


Thus the kernel K(z,t) is constructed such as for every z, z > 0 


J K(z,t)A(t)dt = 0. 
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Note that if \(é) is a quasianalytical function in the domain t > 0 then from the equality 
Jx K(z,t)À(t)dt 
0 

by Theorem 1.1 it follows that A(t) = 


On the uniqueness of solution of the Fredholm first kind integral equation 


Here we consider the multidimensional first kind integral equation with singularity on the 
diagonal or with this properties after repeated differentiation. Let R”be the n— dimen- 
sional Euclidean space for x = (21, 22,...,2n), and let Q, w be balls in R” with: 


Q={r:|z|< R}, w={z:]z|<r}, R>r>0. 


With £ we will denote the points of sphere annulus C = Vo, here the overbar, as usual, 
defines the closure of the set. Consider first the integral equation 


J 46.932) =u(é), dz = dzidz;...dz,, (1.7) 


where kernel M(z,&) has properties: 
1. Mí(z,£) is given for all z € R”, € € R”, continuous for € # z. 


2. For every q > 0 equality 
i lim , MG) —£)* = oo 


holds. 


As example of such kernels we can take M(z,£) = exp(e/ |x — €|) |z — €|°, where € > 0, 6 
can be any number. For € = 0 such kernels take place in inverse problems of potential. 
Let (A), be a set of functions A(z), z € ©, presented by A(x) = a(z)g(z) + b(z) where 
a(x) is some analytical function in closed ball ©, and g(x), b(z) are fixed measurable 
functions, so that function g(x) satisfies the inequalities 0 < Mo € g(x) € Mı. 
Note that if g(z) — 1, b(z) — 0 then set of function (A), includes the set of all 
analytical functions in ©. 


Theorem 1.7. Equation (1.7) has no more than one solution A(z) € (AV. 


PROOF. Suppose, that equation (1.7) has two solutions Aj(z) and A;(zx) A1(x) Æ As(x), 
A(z) € {A}a, : = 1,2. Denote as a(x) = a1(z) — a(x). We have f M(z,£)g(z)a(z)dz = 0. 
Show, that if a(z) # 0 then there exists at least one point & € C so that 
f M(z,£&o)g(z)a(z)dz # 0. By this the statement of the theorem is be established. Let 


s be a point on the sphere S = (z : |z| = r), B(s,t) is an n— dimensional open ball 
of radius t > 0 centered at point s € S . Denote by w(s,t) and C(s,t) the intersection 
B(s,t)Nw, B(s,t)nC. 

Show at first that there exists point so € S and numbers to and tı, to > t, , so that 
the inequalities hold 
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1. |a(x)| > A(r — |z|)?, £ € G(so, to) = do , where A > 0, a > 0 are some constants, 
2. M(z,£)g(z)a(z) > ko, x € & Woo, € € C (So, to), ko is some constant. 


Inequality 1 follows from Lemma 1.2, and estimate 2 follows from the continuity of func- 
tions M(z, £), and a(z), and g(z). 

Note that as a(x) = a1(z) — az(z) the inequality a(x) > A(r — |z|)", x € wo is true, in 
another case functions a1, a? may be changed. Then by the condition the kernel M(z, £) 
is continuous in domain é # x and when |z — é| — 0, M(z,£) — oo. Therefore for 
every s € S it exists intersection w(s, t) and C(s,t) so that M(z,£) > 0 when z € os, t), 
€ € C(s,t). By the above we include there exist intersections uo = w(so, to), Co = C (so, to) 
with 


1. M(z,£) > 0, x € wo, € € Co, 

2. a(x) > A(r — |z|)%, £ € wo, 

3. 3. M(z,&)g(xz)a(z) > ko, E € Co, z € Woo. 
Using estimates we have 


[M@,Oge)a(z\az = f M(x,8)g(2)a(c)de + f M(x, €)9(2)alx)de 


w\wo 


IV 


kovo + MoA | Mz, Or — lel) dx, 
wo 


where vp is volume of domain w\wo. Since lim M(z,€)|r —€* = oo, f M(x, é)(r — 
r— wo 


|x|)*dz — co , when € — so . Thus there exists a point ĉo with f M(z,é)g(x)a(x)dz > 0. 


The theorem is proved. 
Suppose now that kernel M(z,€) is indefinite differentiable in the domain € # z. 
Denote as D? a differential operator with respect to € = (61, 2,.--,€n) : 


Qf t£. On 
- Ot oth... od 
Let Mg(z,£) = D? M(z,£). Let kernels Mg(z, £) possess properties: 
1. Mg(z,€) satisfy condition 1, above, 


2. There exists a sequence of numbers {gm}, qm — oo so that 


lim M(z £z - €" = 00, Ifl o m. 
Consider the following equation with such kernel 


J 46.02) = «(0). (1.8) 


Theorem 1.8. Equation (1.8) has no more than one solution Az) € (A),. 
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PROOF. If A; = a;(r)g(z) + b(z), 2 = 1,2 are two nonequal solutions of (1.8), then the 
equalities 


J Ms(z,£e(z)g(z)dz =0, m=0,1,2,... (1.9) 


hold, where a(x) = ai(z) — a;(z), Mg(z,é) = D^ M(z,£&), |B| = m. As in the proof of 
Theorem 1.1 estimates state 


1. Msg(z,£) > 0, £ E€ wp, Ce Cy, 
2. a(x) > A(r — |r|), r € o, 
3. Me(z,£)g(z)a(z) > km, € € CP, £ € o No. 


From this and the inequality above follows 


nro > kmvo + AMo J M,(z, €)(r — |z|)". 


If the number m is sufficiently large, then f Ma(z,&)(r — |z|)?dz, where € — so tends 
wo 
to infinity. Thus there exists a point £o such that f Mo(z,£o)a(z)g(z)dz = 0, |B| = m, 


which is in conflict with (1.9). The theorem is proved. 


The example of nonuniqueness of the solution. 


Let the kernel M(z,£) = 1, where |£| > r, |z| < r/2 and let A(z) # 0 be an infinite 
differentiable function in circle w = {z : |z| € r} so that 


1. A(z) 20, z € {x: |z} > r/2), 


2 A(z)dz —- 0. 


iz «7/2 


It is obvious that f M(z,£)A(z)dz = 0. 


1.5 ON THE UNIQUENESS OF THE SOLUTION OF INTEGRAL 
EQUATIONS OF THE FIRST KIND WITH ENTIRE KERNEL 


In this section we consider the multidimensional integral equation of the first kind 


w(x) = J F0», zED, (1.10) 
R” 


where D is a domain in the real Euclidean space R”, n > 1, k(z, y) is a complex-valued 
kernel in R” x R”, represented in the form 


k(z,y) = [ (pa) e o7 dpag. (1.11) 
B 
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Here an integrand function (p,q) is complex-valued, continuous in R^ x R^, finite with 
compact support B. 

It is known (Ronkin, 1974), that the kernel k(z, y) is continue on an entire analytic 
function of exponential type. In the present section we give conditions to domain B and 
function p(z, y), guaranteeing the uniqueness of solution A(y) of equation (1.10) in the 
class of continuous complex-valued finite functions. 

Let a(t) > 0, o(0) = 0 be a continuous function in R!, e > 0, 6 > 0 be fixed numbers 
and let space R"(p) be defined in the following way 


R"(p) -(a:(pa) «€R'). 
Consider set 


w(p) = fa $37 =g) + < alpn), dq» o} ; 


The set w(p) is an open hemiball in R”(p) of radius a(p,) centered at 
Po = (Pı, P2,---,Pn-1,0). Further we propose that domain B and function y(p,q) = 


pı(p, q) + ip2(p, q) have properties: 
1. For any p, 0 < p; < £, i = 1,2,...,n, intersection B(p) = BN R(p) is non-empty open 
set in R”(p) and B(p) C w(p), 


~ n-i 
2. In a domain B = BN (CD Dp +< J the following inequalities hold: 
i=1 


y1(p,q) 20, v:(p,3) 20, «(p q) + vx(p, a) > 0. 


In the case where n = 1 as example of a domain B satisfying condition 1 may be the 
domain with 


B= {(pi,q1):0< qı <a(pi), 0<p <1}, 
where a(t) > 0, a(0) =0, 0 € t € 1 is a continuous function. An example of domain B, 
not satisfying condition 1 is the square 


B= {(pi,q):0<m <1, 0«q«1) 
Theorem 1.9. Let conditions 1 and 2 hold. Then equation (1.10) has no more than one 
continuous finite solution A(y). 


PROOF. It is clear from (1.10) and (1.11) that for any continuous finite solution A(y) 

of equation (1.10), function w(z) is an entire analytic function. Therefore, if w(z) — 0, 

z € D, then w(x) = 0, z € R”. Let w(x) 20, z € R”. Show, that A(y) = 0, y € R^. 
We have 


| ks, y) dy = f |[ eem A(y)dy = 0. 
R” R^ 
Changing the order of integration we obtain the equality 


if [et etm / eii dpdq = 0. (1.12) 
B R^ 
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Let 
ila) = f e™Ay)dy. (1.13) 
R^ 
By such designation equality (1.12) acquires the form 


J olp, a)e*X(a)dpdg = 0. 
B 


Passing on to an iterated integral we come to the relation 


J | J aion erp (1.14) 
Rr |B(y) 


Here B(p), as above, is the intersection of the domain B by the plane R” (p) of variables 
q of dimension n, passing through p and orthogonal to space R” of variables p. Let 


X= f eo. (oie. 


B(») 


In this case we rewrite equality (1.14) as 


J A(p) e" dp =0. 
R^? 


Thus the Fourier transform of continuous function \(p) is equal to zero. Therefore \(p) = 
0, p € R” and we have 
J erlad = 0. (1.15) 
B(p) 

Because under the hypothesis of the theorem the sought solution A(y) is finite continuous 
function, then it follows from (1.13) and by the Paley-Wiener theorem (Ronkin, 1974) 
that A(q) = A1(q) + iA2(q) is an entire analytic function. 

Therefore, real functions A(q), 1 = 1,2, are entire analytic functions too. 

Determining the imaginary and real parts in (1.15) we obtain equation system 


S (à — p2ă2)dq = 0, 
B(p) 


f (pire 92À1)dq = 0, P 
Bip) 
Show that \1(q) = 0, À;(q) = 0, q € R^. From here by virtue of (1.13) equality A(y) = 0, 
y € R” follows. We first prove the statement concerning the analytic functions. 
Let po = (p, p9,..., p9 ,,0) be a fixed point, r > 0 be a fixed number and w be an 
open hemiball of radius r > 0 centered at po: 


n-i 
w= [iden «20. 
i=l 


Denote as {w} the set of all open hemiballs centered at (p1, p2,- - - , Pn-1,0) and such, that 
wW <W. 
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If f(q) is a real entire analytic function, then there exists an open hemiball 2 € {w} 
such, that for all g € 2 one of two inequalities holds: either f(q) > 0, or f(q) « 0 

Suppose, that this is not so. Then in every hemiball w € {w} there exists a point q” 
such, that f(g”) = 0. In particular equality follows from this 


f(a)l, co =0, lal <r. 


Because f(g^) = 0 and f(qg) = 0, where gf = (q?, q2,..., q2 ,,0) then in every hemiball 
there exists a point q^ such, that 


ô 
e =0, P Ewe {uw}. 
Os q=9" 
ð 
This leads as above to equality a. =0, |g| «r. 
n |g,.—0 
Hence, for any m, m > 0, : 
"nu 
a T 0, lel <r, 
Oqm qn=0 


which contradicts the relation f(q) # 0. The statement is proved. 

If we suppose now that, although one of entire functions M(q ), ilg ), satisfying equa- 
tion system (1.16), is not identically equal to zero, then by virtue of the statement proved 
above, there exists € > 0, € < € and a point p such, that for any 


q € w(p y= {acd Ee pi — qi)? +4? < alpn), n> 0}, Ücpi;«é, i=l 2e 


one of four relation holds 


1. Ài(a) > 0, Ã2(4) > 0, A (q) + Aa(q) > 0, 
2. (q) € 0, (q) € 0, i (g) + Às(q) < 0, 
3. Ai(q) 2(4) > 0, Ar(q) — Às(a) < 0, 
4. \i(q) 2 0, A2(q) € 0, ig) — Aa(q) > 0. 


Under the hypotheses of the theorem at & « 6 for all (p,q) 


n=l 
(pg) E€ Bn foa togmtq« Jt 
1 
the following relations take place 


v1(5,4) 20, v2(p4)2 0, «i(p,q) + p2(p,4) > 0, 


and 
B(p) = BOR (p) <u(p), 0«pi <£. 


That is why, together with inequalities 1-4, at least one of equalities (1.16) cannot be 
performed for all p € R”. The theorem is proved. 
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1.6 EXISTENCE AND UNIQUENESS OF A SOLUTION TO 
AN INVERSE PROBLEM FOR A PARABOLIC EQUATION 


Problems of finding the coefficients of differential equations from information about their 
solutions are called inverse problems for differential equations. As a rule, inverse problems 
are nonlinear. Their study, especially in the multidimensional case, is often connected with 
significant mathematical difficulties. We will give a method for studying solvability ques- 
tions for certain nonlinear inverse problems for differential equations, using a parabolic 
equation as an example illustrating this method. 

With the aid of the Fourier transform we are able to reduce the inverse problem under 
consideration to a boundary value problem for a nonlinear integro-differential equation 
which is fully acceptable for investigation. Methods of potential theory are applicable 
to the boundary value problem thus obtained; under the appropriate restrictions on the 
data, this leads to determination of solvability, uniqueness, and stability conditions for 
this inverse problem. 

In a domain Q = (—00,o0) x D, D C R”, OD — Ts € C, we consider the equation of 
parabolic type 

p(z)u, — Au = 0. (1.17) 


Inverse problem. Find function u(z,t), p(x) > 0 satisfying (1.17) such that 


ulp = es, t), P=To x (—00, oo), (1.18) 
ulizo = w(t), zE€D, uo(z)],.., = lio; (1.19) 
V.ucL'(Q, ueL?*(—-oo,00: W(D)n C(D)), (1.20) 


u, € L*(—o0,00): C(D), p(x) € C(D). 


We introduce some notation. G(z, y) is Green's function of the Dirichlet problem for 
the domain D, ü(z,£) is the Fourier transform of a function u(z,t) with respect to t and 


(D) = max f IG(s, v)| dv. 
D 


First we consider problem (1.17) and (1.18). 


Lemma 1.3. Let $(£,s) € W;"(To) NC(To); à = 0, |E] > R, and s € To. Then, for any 
real-valued function p(x) such that 


p(z) €C(D), |p| < 1/Ru(D), 
there is a unique solution to problems (1.17) and (1.18) such that Vu € L?(Q), u € 


zm : W2(D) n C(D)), us € L*((—oo,oo) : C(D)) and à(z,£) = 0, |£| > R, 
rcli. 


A proof of Lemma 1.3 is based on the study of the boundary value problem 
Av-iép(z) —0, vl -$(£), z€D, 


and is carried out as in (Bubnov, 1984). 
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EXAMPLE. Let A(£) and B(£) be continuous functions, 0 € € € Ro, Ro > 0, let p bea 
positive constant, and let u(z,t), z,t € R! be the function 


ued- foe (fE) fm (E eo) 40+ moo (E e] 


that is entire analytic in t and satisfies the heat equation 


Qu O*u 


at? da” 
According to the Paley-Wiener theorem, 
a(z,£)=0, |z|» Ro. 


Our study of the inverse problem is based on the following lemmas. Let 


ALS sup J plan, 
zé ôn 
To 
A=min|- f ie f e2Zarvae], wo = iD) vuol R. 
D PA. On ? d 
« 0 


Lemma 1.4. Let the conditions of Lemma 1.3 hold, and let 
. .. [ .0G 
Ay = min - j it J fiat » 0, 
EKR To 


Az > wo + R(8w041)!?, 8R?°A > wo. 


Then the inverse problem (1.17)-(1.20) is equivalent to the boundary value problem for 
the nonlinear integro-differential equation 


i£Auo(z,£) —— 0, vl, = vo(s,£) (1.21) 


SG T 
KI«R 


in the class of functions v € C(D)N W2(D), v(z,£) 20, |z| > R, z € D; 


min 


= J d > Ro > 0, 


KKR 


1 
Ro = 3 wo + Ap + [lwo — A2} = 8RwoAi| "| . 


Taking the nonlinearity into account and using methods of potential theory, we prove the 
following result. 


20 Chapter 1 
Lemma 1.5. Let the conditions of Lemma 1.4 be satisfied. Then there is 6(R, A1, A2) > 0 


such that for any uo with 0 < Aug < 6 there exists a unique solution v(z,€) to problem 
(1.21) that belongs to the class 


B= T e€ C(D)NW3(D), v(z,&)=0, l> R, zeD, 


bege aA 


EINE 
a 2 S | d 


From the lemmas formulated above we get our main result. 


Theorem 1.10. Let the following conditions be satisfied: 
é(6s) e W"(To) AC(T), $=, lt] » R, 


lizo = uo(z)lo , 0 < Aug € C(D), 


: [0G E 
min l- J i€ J 5 dod | = Aa > 0, 


lélk<R To 


Ag > wo + R (8w9A,)'/? , 8R!A, > Wo. 


Then there is a 6(R, Ai, A2) > 0 such that for any uo with 0 < uo < 6 the inverse problem 
(1.17)-(1.20) has a unique solution(p, u); moreover, 


im Auo(z, €) 
Me) = — pies DAE 
EKR 
u(z,t) = EGET 


where v(z,£) is a solution to problem (1.21). 


REMARK T. his theorem can be briefly formulated as follows: if in addition to the con- 
ditions of the first boundary value problem for a parabolic equation we require that a 
solution be regular in t, —oo < t < oo, then, under suitable restrictions on the data of 
the boundary value problem, one can also uniquely find the function p(x) occurring in 
the equation, i.e. one can solve the inverse problem. From the method it is clear that 
this fact permits a significant generalization for other boundary value problems and other 
evolution equations including equation with variable coefficients. 
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1.7 ON UNIQUE SOLVABILITY OF AN INVERSE PROBLEM FOR 
A PARABOLIC EQUATION 


We formulate conditions for the unique solvability of the Cauchy problem for a semilinear 
integro-differential equation, and offer a method for its investigation (Anikonov and Belov, 
1989). Such problems are reductions of inverse problems for parabolic equations in cases 
where the Fourier transform of their solutions with respect to chosen variables exists and 
some other conditions are satisfied (see (Anikonov, 1986)). Analogous inverse problems 
for parabolic and hyperbolic equations were investigated by different methods in (Bubnov, 
1987a,b). 
Let (z1,..., 24-1) be a point in the Euclidean space R?^!, 


Wrest) = {(t, z) |ti t€ to, z€R'"!)astrip in R^, 

G[,,5] = (62.2) (tz) E€ Mep 2€ R'} a strip in R^, 

G° = {(z,2)|2eR™, 2€[-aal), 

Giu] = {(t,2,2) | (tr) E€ Mp) zE [-a,a]} and a = const > 0. 


Let us examine the following equation in Gyo,rj: 


Ou(t, x, z) @u(t, zx, z) Ou(t, x, z 
ACA = Multaa) € ECEE eat eh 4 o(t,2,2), (1.22) 
where 
n-1 ^u n—1 


L(u) = + T Ero. 


aij —— 
Gc. Or0z; 1L 


üj- 


We assume that a;;, b; and c are sufficiently smooth functions of t given on [0, T] and, in 
addition, c(t) € 0 and 


n-1 
per < > aj(t)&€j, VEER, te€([0, T], p= const » 0. (1.23) 
ij-1 


The function q(t,z, z) is given on Gyo,rj. 
The function a(t, x) (0a/0z = 0) has yet to be found. We assume that the following 
condition is satisfied: 


Uleo = (hz) (t,x) € Tir] (1.24) 


and that the Fourier transform of u(t, z,z) with respect to z exists: 


D : 
W(t,z,y) = on J u(t,z,y)e "dz, 

ed (1.25) 
u(t,z,y) = | W(t, z,y) e^*dy. 


Let us introduce a Cauchy condition for equation (1.22): 


u(0,2,2) = uo(z,z), (z,z) € R”. (1.26) 
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With the help of conditions (1.24) and (1.25), problem (1.22) and (1.26) can be reduced 
to a Cauchy problem for a integro-differential equation, not containing a(t, x). Indeed, by 
applying the Fourier transform to equation (1.22), setting z = 0 in (1.22) under conditions 
(1.24) and (1.25), we obtain the equations 


iud = L(W) — y?W +iyaW +Q, 


dp 2 oo ; oo 
X - Le) K. A Warsia f AW4A + q(t, 2, 0), 


from which we eliminate function a and obtain the problem 


o = LW) -pW + 4 — fos f swal +Q, 
f. AWdA 


WY(0, z, y) = Wo(z, y), (z,y) ER". 


(1.27) 


—00 


Here 


1 7 —izy 
Q(t,z, y) > x J q(t, z,y)e dz 


is the Fourier transform of q and 


* F 
Wolz, y) = z- f uols, y)e™ dz, 


O(t,z) = Po — L(y(t, z)) — g(t, z, 0). 


We examine the simplest case, when Q and Wo are real-valued functions with compact 
supports in y, belonging to a finite segment [—a, o]. Then (1.27) reduces to the problem 


d = 4) -W+ P das fias] +Q, (1.28) 
t f \Wdd A 
W(0,z,y) = Wo(z,y, z€R"', y€[-o,o]. (1.29) 


Equation (1.28) is a nonlinear integro-differential. It consists of partial derivatives 
as well as integrals of the solution with respect to the parameter, which also appears in 
the coefficients and the initial data of the problem. These are interesting mathematical 
objects which have not yet been studied. We examine the solvability of problems (1.28) 
and (1.29) and some properties of its solutions. 

Below we denote by C*'(Gf, ,.;) the space of functions that have continuous derivatives 
in Gt, ta] in the space variables of order up to and including k, and a derivative of order 
l in the variable t. 

As to the input data of problem (1.28) and (1.29), we presume the following: 


are bounded in absolute value (1.30) 


$,q, Wo € (Gig); all kth derivatives in z and y 
by constants ck, — k — 0,..., 4, respectively in Gs 
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Q(ti,z,y 20, We(zy)y20, teloT] ze R"!, y€í-o,a] (1.31) 
Im (2) = f Wi(s,)dA 2520, zeR™. (1.32) 


Let us split problem (1.28) and (1.29) according to the method of the weak approxi- 
mation (Yanenko, 1971; Belov, 1986): 


18W, 
2 8t 


1 
=L(W,), nr<t<(nt+ 3 )r, (1.33) 


: ; 1 
r~ —y?W, + e c fa + j xwa +Q, (n4 5 \r<t<(n+1)r, (1.34) 
f. AW,dd 


-Q 


W,(0, x, y) = Wo(z, y), (z, y) EG”. (1.35) 
Here n = 0,..., N — 1 and N, = 4t, 0< t € T. 
We start with the Cauchy problem 


1 Ov(t,z,y)  —— yv(t,z,y) | 


ari mes os [estan] — y^v(t,z, y) + Q(t,z, y), (1.36) 
2 t f. Av(t, 2, y)dÀ 


-a 


v(to,z,y) = volz, y), (z,y) € G°. (1.37) 


By (1.30’), (1.31) and (1.32) we denote conditions (1.30), (1.31), and (1.32) with Wo 
replaced by vo. 


Lemma 1.6. Let ® < k, k = const, and let conditions (1.30’)- (1.32’) be satisfied. 
Then there exists vo, 0 < vo € T, such that problem (1.36) and (1.37) is uniquely solvable 
in the class C*! (Gg, 1.4.1) for any v € (0, vo] and any to € [0,T — v]. vo depends only on 
the constants 6, k and a. 


PROOF. Let us linearize problem (1.36) and (1.37) by making a shift 0 > 0 with respect 
to t of some terms in (1.36): 


2 ot f v*(t — 0,2,y)dÀ $a 
E y^v*(t,z, y) $ Q(t, x, y), ss 
v4| = vo(z, y), (z, y) e G*. 


tSto 

It can be shown that there ons vo depending on 6, k, and a such that for all v € (0, vo] 
and to € (0, T — v] the solution v^ of problem (1.38) exists in the interval [to, to + v] for 
any fixed 0 c (0, 6} and uniformly in 6 c (0, 6] and (t,£) € Ife 4. 


[rot -8,2,y)Ad PERS (1.39) 


—a 


Estimate (1.39) allows us to obtain uniform boundedness (in 0) of the family of solu- 
tions (v?) and their first derivatives with respect to z and y. Hence, as 0 — 0, uniformly 
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in every compactum of Gite to dup the functions v* converge to v being solution of the 


problem (1.36) and (1.37) in C*"(Gft, 14,3). As a result 


J Av(t—6,2,y)33 2 p> 0, (t,£) € Hp toss}: 
This allows one to prove the uniqueness of the solution in the class specified above. 
The lemma is proved. 


REMARK 1.1. In consequence of condition (1.317), the inequality v(t, z,y)y > 0, (t, 2, y) € 
Gt, toy) IS realized. 


REMARK 1.2. Let v be the solution of problems (1.36), (1.37) in the interval [to, to + v], 
the existence of which is guaranteed by Lemma 1.6. As follows from remark 1.1, we can 


apply the generalized mean value theorem (Fikhtengol'tz, 1966) to the integrals f A?vdàÀ 
and f Xvdà: 


J M v(t,z, A)dÀ = m(t, z) J Av(t,2,A)dd, k=1,2. (1.40) 


-a -a 


It is clear that m, are smooth functions that satisfy —a < mı (t,£) < a and 0 < m(t, z) < 
a’, By multiplying (1.36) by y, integrating the result over [—a,a] and then applying 


(1.40), we find that 
Iy(t,z) = J Av(t, z, A)dÀ 


-Q 


is a solution of an ordinary differential equation 


dI, (t, 
x z) = m(t, r)®(t, z) + (mi(t, z) ^ m(t, x))L (t, x) + Ig, 
where z is a parameter. Let us consider the equation with constant coefficients 
dj(t 
te = —ka - o’ i(t). (1.41) 


Its solution j(t), satisfying the initial condition j(0) = 6, is monotonic decreasing and 
strictly positive in an interval [0,t,]. The value 4, depends on 6, k and a. Since m6 > 
—ak and m? — m; > —o? we have 
Iv(t,2) > (0 2 (t), to StS ta (1.42) 
provided that 
I (to, 2) = J Avo(z, A)dA > j (to). 


—« 


Theorem 1.11. Let || € k = const, and let conditions (1.30)-(1.32) be satisfied. 
Then there exists t, € (0, T] such that problem (1.28) and (1.29) are uniquely solvable in 
C? (Gio, ]). The value of t, depends on the constants ô, k, and a. 
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PROOF. By applying the maximum principle to equation (1.33) (see conditions (1.22), 
Lemma 1.6, and the inequalities (1.42) of remark 1.2), we can easily show that a solution 
W,(t,z,y) of problems (1.33)-(1.35) exists in the interval [0,t,] for any fixed 7 € (0, ro] 
and, uniformly in 7 € (0, ro], 


Iw,(t,2) > j(t) > j(ts), to<t<t., (t,£) € lot) (1.43) 


where j(t) is a solution of (1.41) and j(0) = 6. Estimate (1.43), the maximum principle 
for equation (1.33), and condition (1.30) allow us to prove the estimates, uniform in 
T € (0, To], 
x sup |D‘w, + >D sup [D.D^w, 
l4 Cio. c. I&|«2 Gio... 
where D, = 0/dt, D* = dl /dch... ck ay'n, k = (ky,..., kn) is a multi-index, and 
|k| = kı +... + kn. Because of (1.44) as in (Yanenko, 1971) and (Belov, 1986), we can 
prove that, uniformly for each compactum of Gi}, W- converges to W € C?1( (9...) 
along with its derivatives of first and second order with respect to the space variables z, y, 
and is a solution of problem (1.28) and (1.29). The uniqueness of the solution is proved in 
a standard way by proving that the difference of two possible solutions of problem (1.28) 
and (1.29) from the class C^! (Gt, , ;) is identically zero. The theorem is proved. 


<c, (1.44) 


REMARK 1.3. If ® = 0, then we can take k = 0, and in this case the solution j(t) of 
equation (1.41) with the initial condition j(0) = ô is strictly positive in the whole interval 
[0, T]. This guarantees the solvability of problem (1.28) and (1.29) in this interval. 


1.8 FORMULAS IN MULTIDIMENSIONAL INVERSE PROBLEMS FOR 
EVOLUTION EQUATIONS 


Let D be a domain in the real Euclidean space R?*"*! of variables (z,z,t), 
T = (21,...,2m), Z = (Z1,---,2n), t € Rl. We consider the inverse problem of the simul- 
taneous determination of two complex-valued functions u(x, z, t) and A(z,t) (OA/Oz, = 0, 
2=1,...,n) in D which satisfy 


u = Au + Lu + Az, t)Bu. (1.45) 
ulo = uwo(z,z) ulo = f(z,t). (1.46) 


Here A, B, and L are linear differential operators such that A and B act with respect to 
zx and are independent of (z, t) while L acts with respect to z and is independent of (z, t), 
the functions uo(z,z) and f(z,t) are given, and (z, z,t) € D. 

As examples of operators A, B, and L we may consider linear differential operator 
with constant coefficients Aa, Ba, and La such that 


A=) AD, B= Y B.D, L-.L,Dj. 
In particular, with 


m o? 8? 
A=) a ba n=], 
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arises the inverse problem for parabolic equation, namely, to find functions u(z, z, t) and 
A(z, t) satisfying 

Ou 
Ano Au + À(z,t)Bu, ulo = uo(z,z),  u|.g- f(z, t). 

For Bu = ¢(z,t), where y(z, t) is some function, the inverse problem (1.45) and (1.46) 
consists in finding a source function J of the form J = A(z, t)p(z, t). 

We study the case when the initial state u|, , = uo(z, z) can be factored in the form 
of a product a(z)b(z), and establish explicit formulas for the solutions u(z, z, t), A(z, t) of 
(1.45) and (1.46), which introduces a constructive element in the theory of such problems. 
In what follows we assume that A, B, and L act in a well-defined way on the functions 
under consideration, and the integral and other transformations used are also well defined. 


Lemma 1.7. Suppose that the function f(z,t) such that Bf # 0, that y(z,t) is a 
solution of the equation 


Ot = Ly, 
and that qo(t) = y(0,t), (x, z,t) ED. Then the functions 
of dyo 
SAS) vo- I 
f(x, t)p(z,t) e dt 
u(z,z,t) = ——————, A(z,t)- 
( ) polt) (nt) Bf eo 


satisfy equation (1.45) in D and the conditions 
u|,.o = f(z,t), ulizo = a(z)b(z), 
where a(x) = f(z,0) and b(z) = v(z,0)/vo(0). 


Let L be a differential operator with constant coefficients and c(w) its symbol; that is 
L(e'«*) - c(w) e, w € R*. 


Lemma 1.8. Let f(z,t) and q(w), (z,t) € D, w € R” be complex- valued functions 
such that the expressions 


of = c(w)t+iwz 
a= AS, BY, il e q(w)dw 


are meaningful and Bf #0, f e°l)q(w)dw 4 0. Then the functions 
R^ 


Flayt) f ertonqo)da 
u(z, 2z,t) = RT 
J e q(u)do 
R^ 
à 
(a = AJ) ET -f f eg (w)e(w)dw 
Ma, t) = ——— 2 _ 


Bf J eM g(w)dw 
R^ 
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satisfy equation (1.45) and the conditions u|, , = f(z,t) and ul,_) = a(z)b(z), where 


J 40999 
a(z) = f(z,0) b(z)-E* 


J q(w)dw 


R^ 
The following assertion holds. 


Theorem 1.12. Suppose that L — tz i 


(1.45) and (1.46), where a(z) = b(z,0) ind b(z) is a continuous bounded complex-valued 
function in R” with b(0) = 1. If the conditions in Lemma 1.7 hold, then the solution 
u(x, z,t), A(x, t) of the inverse problem (1.45) and (1.46) can be represented in the form 


z € R^, and uo(z, z) = a(z)b(z), in problem 


f(2,t) f b(y) exp(—ly — 21? /4t)dy 


u(z,z,t) = R” 
J exe 1t)y 
n 
X2.) = (E - AJ) COSE ILOG (-£) ay 
- ^w 5 gays | alee (-5) al 
x 


(ores J by) exp (-£) a) 


COROLLARY . Suppose that Bu = ¢(z,t) in (1.45) and that the conditions of the 
Theorem 1.12 hold. Then the source function J(z,z,t) = v(z,t)A(z,t) and the solution 
u(z, z, t) of (1.45) can be found from the data (1.46) by means of the formulas 


ol!) = rers J b(y) exp(—ly — 21? /4t)dy, 


F(z, thelz,t) 


u(z, z,t) = ws (t) = v(0, t) 
dyo 
CV. 
ot dt 

J(z,z,t) = TO e(z, t) 


From Lemma 1.8 we deduce the following assertion. 


Theorem 1.13. Suppose that uo(z, z) = a(x)b(z) in (1.45) and (1.46), where z € R^, 
a(z) — f(z,0), b(z) has a Fourie transform and b(0) — 1, and that the conditions in 
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Lemma 1.8 hold. Then the solution u(z,z,t), \(z,t) of the inverse problem (1.45) and 
(1.46) is given by 
f(s,t) f ee) 

R” 


u(z,z,t) = 


Hu 


f e*t a(w)dw 


R^ 


of c(w)t c(w)t 
(2 — Ar) fe ( "tey fs "t a(w)c(u)do 


A(z,t) = 


X 


Bf | e*"*q(u)ào 
R^ 


where 


ao) = f (y) eds. 


R” 
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Inverse Problems for Kinetic Equations 


2.1 KINETIC EQUATIONS 


Kinetic equations characterize the continuous nature of motion of an object, they are 
the fundamental equations of the natural sciences. These equations are widely used for 
qualitatively describing physical, chemical, biological and other processes at the micro- 
scopic level. In view of the great role kinetic equations have, as a rule, integro-differential 
equations can be used. These are often nonlinear relative to one or several distribution 
functions F(z, p,t) describing the dynamic state of an object with respect to the spatial 
position z, the momentum p and the time t. The Liouville, Boltzmann, Vlasov equations, 
a chain of the Bogolyubov equations, the Fokker-Planck equations, various equations of 
radiation transport, quantum equations for the density matrix and the Wigner function 
(see Alekseev, 1983; Bogolyubov, 1981; Case and Zweifel, 1967; Ferzeiger and Kaper, 
1972; Haken, 1987; Klimontovich, 1982; Kompaneets, 19777; Lifshitz and Pitaevskij, 1979; 
Prigogine, 1980) are examples of kinetic equations. 
Below F(z, p,t) is a distribution function. 


1. Liouville equation: 


A OF OH OF 2) 


F 
3 Hess {F,H} = Sm Op; 2; 


2. Boltzmann equation: 


OF 
RP -Pj = St F. 


j=l 


3. Transport equation: 


F 
ad 4 grad, F + oF = Jis, Q, V)FAN +q. 
4. Transport equation in tomography: 


oa E big F ps nr t, Er, Q, Q) FAE'dQ' +q. 
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5. Kinetic equation in problems of integral geometry: 
Qgrad, Fü--oF =0, InF =% > grad, YN = —c. 


OF OF OF 
f La 
Iz opta sing + K(z,y9) 57 (z,y, v) 


K - first curvature. 


6. Fokker-Planck equation: 


A +{F,H} = J [W(p+49,9)F(z,p+49,t) — W(p,a)F(z, q,t)] dq. 
R” 


7. Quantum equation: 
OF i 1 
2 =. oa — S hne.t 
gp US MISERE aan! I (s ; t) 
R2n 
—-$ (s + ; hy,t)| e(»-») F(x, p t)dp'dy + St F +q. 


Kinetic equations with regard to quantum effects have the form 


oF i 1 
Sr +F, H} +aF = ayn J je (z- 59.1) 


~ (s + : hy, 3] eiv(r-P') F(z, p t)dp'dy + St F + q(z, p, t). (2.1) 


Here F(z, p, t) is the quantum distribution function, H is the Hamilton function, ® is the 
averaged potential, a(z, p, t) represents the absorption, A is the Planck constant, St F is 
the collision integral, q(z, p, t) is a function that characterizes the sources. 

Given the appropriate information at the boundary, the main direct problems for 
equations of type (2.1) describe the evolution of an initial distribution. As a rule, all 
the coefficients in equation (2.1) are considered to be known. Analysis of physical prob- 
lems connected with kinetic equations (for example problem of tomography) shows that 
in addition to the distribution function F(z,p,t) certain coefficients contained in these 
equations should be ascribed to the quantities sought. 

The other words some coefficients are unknown in these equations. The problem of 
determining some functions involved in the kinetic equations, e.g. a or 6, apart from 
the distribution function F'(z, p,t) is conventionally called the inverse problem for kinetic 
equations (see Anikonov, 1975; Anikonov, 1978; 1984; 1985; 1986; 1987; 1989; Anikonov 
and Amirov, 1983; Anikonov and Bubnov, 1988; 1989; Lavrent’ev and Anikonov, 1967; 
Lavrent'ev et al., 1971). In these problems the main information is the trace of the 
distribution function on a certain (z, p, t)-manifold, e.g. z, = 0 or |z| = 1, etc. Physically, 
inverse problems for kinetic equations are connected with determining interaction forces, 
absorptivities, sources, scattering indicatrices and other physical quantities. Some inverse 
problems of this kind may be formulated as control problems (Anikonov and Bubnov, 
1989). 
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We repeat that our basic information for inverse problems is the function Fo(s, p, t) = 
F(z,p,t)lr, s ET, p€ D CR”, tı <t < tz, where T is a manifold in R” of variable z, 
for example 

T-í(z:|zr| 21) or (z:z,-0). 
The other words we must know the trace of the distribution function F(z, p,t) on T. 

The studies concerning predominantly multidimensional inverse problems for kinetic 
equations suggested and developed by the author (Anikonov, 1978; 1984; 1985; 1986; 1987; 
1989; 1991; Anikonov and Amirov, 1983; Anikonov and Bubnov, 1988; 1989; Anikonov 
and Pestov, 1990) are reported in this work without attempting to survey the whole field. 


2.2 AN EXAMPLE OF AN INVERSE PROBLEM FOR KINETIC 
EQUATION 


We consider the inverse problems of determining two functions (F,q(z,t)) or (F, A(z,t)), 
0q/Op = 0, OA/Op = 0 in the region |p| < oo, |z| < a, |t — to] < b, a > 0, b » 0, 


OF OF OF 
Ot + P + Mea = q(z, t), (2.2) 
F|z=0 = Fo(p,t), ipl <œ, |t—to| < b. (2.3) 


Equation (2.2) is the one-dimensional kinetic equation (Vlasov equation for q = 0) 
A(z, t) = —06/0z, and 9 is the potential. 


Theorem 2.1. Suppose that Fo(p,t) is an analytic function in the region |t — to| « 6, 
|p| < oo while for any t, Fo is the entire function of exponential type in p which belongs 
to the space Lz. Then there exist such numbers a, > 0, bà > 0 that inverse problems 
(2.2) and (2.3) to determine (F,q) have an analytical solution in the region |x| < ay, 
|t — to| < bi, |p| < co. 


Theorem 2.2. Let the conditions of Theorem 2.1 be met and 


OF 

=| #0. 

Op j 
Then there exist such numbers a, > 0, b > 0 that the inverse problems (2.2), (2.3) to 
determine (F, A) have an analytic solution in the region |z| < ai, |t — to] < b1, |p| < oo. 


2.3 ONE-DIMENSIONAL INVERSE PROBLEMS 


In one-dimensional inverse problems explicit inversion formulae can often be obtained 
for the solution. The following example illustrates this situation. Let us consider the 
one-dimensional kinetic equation (2.1) with 


@=0, StF=0, g=0, H= 5p + &(2). 


In the case it runs as ar oF 
— — é ‘ tas = 
az? $'(z) Op T a(z)F =0. (2.4) 
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Given the information Fo(p) = F|.-o, p € R!, in equation (2.4), we shall seek the 


functions F(z,p) and a(x) (Ga/Op = 0). Let us assume that the function Fo(p) > 0 is 
differentiable, and introduce the functions 


fle) =F (yar), s) =5 (y2), 


ld 7 g(z)dz 
= — — > 
oR) = 5, dp]vp-z P^ 0, 


where 


joys BR. qu By 


Theorem 2.3. If the function Fo(p) is continuously differentiable, and Fo(p) > 0, 
(x) > 0, x > 0, 6(0) = 0, then the following formulas hold: 


a(x) = 2b (&(z)) $'(z), 


Pen) =e rra) + f ra] r20, pcR!. 


REMARK 2.1. In addition, assume that the functions (x) and F(z, p) satisfy the relation 


oo 


$"(z) = 4rc J F(z,p)dp. 


—O00 


Then according to the formulas of Theorem 2.3 under appropriate restrictions on Fo(p) 
we obtain the equation for the potential $(z) 


$"(z) = P(8), 
where function P is known and depend only on the initial information Fo(p) by the formula 
T 1 021 1 
P(z) =4no f exp ICCEDESPCCEECEDE, dp. 
S = 


In view of this, given $(0), 6(0) and ’(0) the function (z) may be considered as 
the sought function along with a and F. Here we omit the formulation of results. 
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2.4 MULTIDIMENSIONAL INVERSE PROBLEMS 


This section deal with the multidimensional version of the inverse problem considered in 
section 2.1. In the multidimensional case the unknown coefficient also turns out to be 
found by an explicit formula used as a solution of the Cauchy problem for an equation 
of the evolution type with respect to a particular spatial variable. Here we also restrict 
ourselves to one example. Let H — lp in equation (2.1), a and F be the sought functions, 
and the operator St be given. We define the operator A by the formula 


n-1 4 
Ar = E OF oa J [5] e80- F(z yl t)dp'dy — St F — q. 
R?” 


at 2 ga” (27) 


Let the function 
Folz’, p,t) = Fl, z' = (21,25,..., 24-1) 


be given, and 


Op. 
Note that the operator A does not involve differentiation with respect to z,, it is inde- 
pendent of a and known, if ® and St are given. 
Theorem 2.4. Let a twice differentiable function F » 0 be the solution to the Cauchy 
problem for the evolution equation with respect to the variable z,: 


1 


= IE 


0 


-[r AF| (zp,)dz, F\z,=0 = Fo(z', p, t). 


Then the formula 


a= (ars =) [F 


is valid. 
REMARK 2.2. Under the condition 


Ab = 4r0 J F(z,p,t)dp 
R^ 


and the boundary condition for $(z,t), the latter may be considered as the sought func- 
tion, along with a and F, similar to the case above. In classes of entire functions the 
solvability of the Cauchy problem similar to that of Theorem 2.4 may be studied on the 
basis of Banach space scales (see Lavrent'ev et al. , 1971). 


Some results in the analytic case. 


In the Euclidean space R?^*!, n > 2 we consider the kinetic equations 


OF OF OF 06 
we (Z Iz t Op; i) S 


Below it will be convenient for us to represent the collision integral St F as follows: St F = 
aF + b where a and b are functions of variables (zx, p,t) which are possibly generated by 
the distribution F(z, p, t). 


(2.5) 
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Let 


(p = arctan (A | q = VP + Pi, 2=2n-1, Y= fn, 


n=1 
OF OF 06 
apum 3 sens Op; Oz; 2). 


j=l 


With these variables (2.5) can be written as follows: 


oF + ae += or + Zils a cos y + — op sin 
pi Sg PETS dnt a PE gy P 
10F [06 9o 
Ice f gi = A 2. 
230 (z cos — 3 sing) + LF aF +b (2.6) 


Below it is assumed that the functions ®, a and b in (2.6) do not depend on the variable 


p, i.e. 
09 = Qa ðb 


9g " Op ap 
We formulate the following inverse problems: 


Problem I. Find the functions F, a and b if the functions Fo = F|y=o and 9 are known. 


Problem IT. Find the functions F, 9 and b if the functions 69 = ®|,29 and Fo = F|,-o 
are known. 


Problem III. Find the functions F, 9 and a if the functions Fo = Fly-0, ®o = Ply=o 
and b are known. 


All functions considered below are assumed to be analytic in the neighbourhood of 
the point (0, qo), qo > 0 of the space R?^*!, and periodic in q with period 2r. 


Theorem 2.5. Suppose that for each inverse problems I-III the conditions of the same 
number below is satisfied. 


1. Foly=o ee Foly=r # 0, 


OF, OF, 
A o D 0 
4 e y=0 E v pEr 7 > 
OF, OF, 
3. ok Bo ke 0. 
P E y=0 * e ? y=n á 


Then the uniqueness of the solution of these problems is guaranteed in the class of 
analytic functions. 


Theorem 2.6. Suppose that 


oo 
Fy- Y uf cos ko + v? sin kg, 
k=0 
where |u®| < Ae-?*, || < Ae-P*, k = 0,1,2..., A > 0, p > 0 are constants and 


conditions 1-3 of Theorem 2.5 are satisfied. Then there exist analytic solutions of inverse 
problems I-III. 
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Consider the inverse problem for kinetic equation in the multidimensional case: 


OF OF OF 


[Ar "is o P rt POE 


find F(z,y,q, p,t), o(z, p,t), if function F|,-o = Fo(z,q, p.t) is known (a, f and St are 
known also, St F is Boltzman the collision integral). 


Lemma 2.1. Function F(z,y,q,p,t) is the solution of problem 


eal [r+ Ea- rans rb amm F|yzo = Fo. 


Theorem 2.7. In the analytic case of the inverse problem of finding F > 0, ø has no 
more than one solution. 


Let us consider the uniqueness of determining the potential 6(z, t) from equation (2.1) 
using the functions a,(p,t), b(p,t), k = 1,2, p € R”, a € t € B, which are defined by 
the relations 


as(p, t) = Fl. b.(p, t) Ss UE az; LL 


j=l 


(2.7) 


zr? 
r-r. 


where zf, x9 are fixed points of the domain D C R”. Assume that c,(p,t) = St ax, and 
that there exist Fourier transforms 9, à, bk, ¢ of the functions 9, ap, bk, c, with respect 
to (x, p) respectively. It turns out that the presence of the quantum term in (2.1) allows 
us to formulate the following result. 


Theorem 2.8. Assume that the functions 


moa z "f (84/00 + b — êr eci"dy, k=1,2 
R^ k 


are defined correctly. Then the formula 


ô e ) E A, e o) = Az e G1) 
h , 


sin2(z2 — 23, z) 
holds. 


REMARK 2.3. The additional information (2.7) permits posing the problem of simulta- 
neously determining three functions F, a, ® from (2.1) using Theorem 2.4 and 2.8. 


2.5 A UNIQUENESS THEOREM FOR THE SOLUTION OF 
AN INVERSE PROBLEM FOR A KINETIC EQUATION 


Let Q be a domain in the real Euclidean space R2"*!, n > 1, of the variables (z,p,t) 
where z € D and D is a domain in R” with a smooth boundary while the variables p c R^ 
and t € R are subject to the inequalities |p; — p?| < a; and |t — to| < b, where a; > 0, p?, 
i=1,...,n, b » 0, and to are fixed numbers. 
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In Q we consider the kinetic equation 


oW ðW ðH OWOH 
r; Opi Opi ðr) P 2.8 
ðt + 2 (= Op; Op; a) A(z, p, t)W (2.8) 


with on the right side A(z, p, t), satisfying the equation 


n gÀ " 
j=l Op; Oz; 

We pose the inverse problem: find functions W(z, p, t) and A(z, p,t) in Q if the Hamil- 
tonian H(z,p) € C?(Q) is given and the trace of the solution W(z, p,t) of (2.8) on the 
boundary T of Q is known, i.e. W|r = Wo(z, p, t) > 0, (z, p,t) € T, where Wo is a known 
function, Wo > 0. 


Theorem 2.9. If in the domain Q the matrices (0? H/Op;Op;) and (—0*H/0z;,0z;), 
i,j = 1,...,n, are positive definite, then the inverse problem has no more than one 


solution F(z,p,t) € c7(Q), A(z,p,t) € C*(Q). 


We sketch the proof. Let F = In W, then we have a linear inverse problem: 


OF OF OH OFQOH 
or (ann ouo] = 
F |r= Fo =In Wo. 
Suppose that F(Z, p, t), A(z, p,t) is a solution of the inverse problem such that F |r= 0. 
We shall show that F(z,p,t) = 0 and A(z,p,t) = 0 in Q. By the hypothesis of the 
problem the function A(z, p, t) satisfies the equation 
n 07 
2/3» yr; — * 


Hence 


z OF " {OF OH OFOH "OF QA "n 0 QA 
2:5 - m; (5 Es E 2.22) gs ED (rz). 


j= z ô 1-1 j=l 


It turns out that 


»m OH OF OF 0H OF OF : = 0 [OF(8F0H OF OH 
2 42, \OpjOp; 02,02;  ÓzjÓr; Op; Op; 2 6 Op; (3x; \ Ox; Op; Opi Aa; 
E e our HE 
255 Ox; |Op; \Ox; Op; Opi Ox; 2 55 Oz; Op; Ox; Op; 
1 4 3 (ðH OF OF ð (OF OF 
- iÈ a (annar) 12S La E in) us 


_ 3 (POF), 8 (AF OF\)_ ð (AA 
Oz; Ot Op; Op; \ Ot Oz; UE Oz; Op; } © 
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Integrating both sides of (2.9) over Q and considering the condition F |r= 0, we obtain 


f z H OF OF _ 9H OF OF 
Op; Op; Ox; Ox; Ox ;Óz ; Op; Op; 


2.) dzdjdt = 0. (2.10) 


Q i,j=1 


By the hypothesis of the theorem the matrices (0?H/Óp;Op;) and (—0?H/0z,0z;) are 
positive definite. Therefore it follows from (2.10) that OF /0z; = 0 and OF/Op; = 0, 
i=1,...,n which with F |r= 0 give the equality F(z, p,t) = 0, (z, p,t) € Q. From (2.8) 
it then follows that A(z, p, t) = 0. 


2.6 THE GENERAL UNIQUENESS THEOREM 


We consider the simultaneous determination of three functions (F, H, A) from (2.1) where 
F is the distribution function, H is the Hamilton function, and À = St F +q. The 
uniqueness of the solution to the inverse problem can be proved under special restriction 
concerning both the data and the sought functions (F, H, A). Let D C R” be a domain of 
variables p with a smooth boundary; Q = D x Di, z € D, and 0Q be the boundary of Q. 
The data for the inverse problems are the functions Fo, A, B, C, Ho defined as follows: 


Fo(s,p,t)=F |r, seT-ÓD, peR", actxf, 
(2.11) 
A(z,p) =F liza, B(np)-Fles C(2,p) = |, Ho = H hog. 


t=a 


Suppose that the following conditions are satisfied. 


1. The gos Hamilton function 2m is analytic in p € R” and H(z,p) € C?(DxR"7). 


aT hE € 0, 5 i9: EDI 2» o|£^, VE € R^, with constant a > 0. 


ij-l 


2. 
P "t az 


3. |D^H(z,p| € (N + |p|)", lal € 2, N > 0, m > 0 are constants, D^ being a 
differentiation operator with respect to variables (z, p). 


4. The known potential 6(z,t) belongs to the class C?(R" x [a,b]), pe R^, a €t «x f. 


n 


IDIS] < (Mi + |pl)y™, Iyl <2, Ny > 0, m; > 0 are constants. 


n OS ! 
6. P» zz; <0, VE c R”. 


7. The given function obeys the conditions 


A € C*(Q), pr d -&&; > 0, Y xa. &6 «0, VEZ O. 


$j-l ij-21 


Qo 


. The sought quantum distribution function F(z,p,t) € C?(D x R” x [a, b]) together 
with its first and second derivatives, rapidly vanishes by 


Jim |D*^Fi|p" 20, Vm»0, lo| <2. 
lp] co 
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9. The unknown function q(z, p, t) is twice and continuously differentiable and satisfies 
the equation 


10. The operator 


is supposed to be known, defined on the set of distribution functions and to map 
this set onto itself. 


11. For any two distribution functions Fi(z,p,t) and Fo(z,p,t), Fi # Fo, the following 
inequality holds 


IG lgrad, (Fi — Fo)? + (F, - Fa) (MF; — MF,)\ dzdpdt > 0. 


Let us make some remarks concerning the conditions stated. Conditions 2, 6 and 

7 imply the convexity of the functions H(z,p), A(z,p) and ®(z,p) with respect to the 
corresponding variables. In the case n — 1 these functions may be as follows 

l p 


=5P, peR', Ace, 


pl < = 
P J^ 
$(z,t) = —,(t)z? + B(t)z + c(t), 4» 0. 


Conditions 3, 5 and 8 characterize an increase and a decrease, respectively, in functions 
H(z, p), $(z, p) and F(z, p,t) at infinity. Condition 8 is usually considered to be satisfied. 
Sometimes the sources q(z, p, t) are independent of p or z. These cases are embraced by 
condition 9. Items 10 and 11 state the restrictions concerning the aprioriunknown operator 
St. Condition 11 is connected with the monotonicity of the operator 


T 2 a" 


which, unlike operator St is considered as given. 
The inverse problem is: find the functions (F, H, à), given (2.11) for equation (2.1) 
and provided that conditions 1-11 are satisfied. 


Theorem 2.10. The inverse problem has no more than one solution. 


PRoor. We first prove that the Hamiltonian H(z,p), x € D, p € R" is uniquely deter- 
mined. 
Setting t — 0 in (2.1), we have 


Clap) UH) = coss. | loeo AG papdy + StF + go(2,p)- — (212) 
R?^ 


Here 


qo = q(z,p,0), [®]o = le (s - ; by, 0) -$ (s 4 5 hy, 0)] ] 
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If H, and H; are two solutions of (2.12) satisfying H. jlag = Ho, we consider the difference 
y = Hı — H4. We have 


(o, A} = St A — St A + qor — dq =o, lag = 0. 
The right-hand side St; A — St; A + qo1 — qo2 = Ao(z, p) can be shown to satisfy 


a By 
ürjÜp; | 


jo 


Indeed, by conditions 9 and 10 


E oa {(St2 A- St; A) T (qo = do2)] = = MA- MA - 0. 
j 


Thus, 
(eA)- X. E =0, elo (213) 
, — ^0 "M3 cM . 
ico Ox; Op; 


Using (2.13) one can verify directly that 


n 2 
EZ a= 5 RERO E sini 
+ o 2 (ees 


Recalling the condition ylag = 0 and making use of the essential fact that the domain Q 
is a Cartesian product D x Di, x € D, p € Dy, we then obtain by integration 


| ^ (ss Op Oy aA 2e te) asap = 0 


À cà \Ôp:ðp; Ox; 02; — OziÓz; Op; Op; 
From this and condition 7 we obtain Op/ Oz; for j = 1,...,n, which in turn leads to 
the relation y(z,p) = 0, (z,p) € Q. Thus, H1 = Hz in Q. Since by condition 1 the 
Hamiltonian is analytic, it has a unique continuation from the domain Q to D x R^, on 
which it will henceforth be taken as known. 

Let (F1, 41) and (f2,.2) be two solutions of the inverse problem satisfying all the 
hypotheses of the theorem. Writing F = F, — F; and q = qı — q2, from (2.1) and (2.12) 
we obtain 


= pr t{hA}= Come [9] ev (77? F(z, p’, t)dp'dy + St; F, + Stz Fz +q, (2.14) 


Fltaa = 0, Fl = 0, Flzer = 0. 
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Differentiating (2.14) with respect to p;, multiplying by QF/Oz;, and summing over J, 


we get the equality 
"OF ð {OF 
FH 
© dx; Op; (s TAS ) 


- oe ap (On aa) pele P) F(x, p', H)dp'dy + St; Fi + Sto Fy tq. (2.15) 


The following relations can be verified directly: 


1 


sea CAE aa La ae 
2 i ot Oz; Op; 2 Ox; \ Ot Op; 2 Op; \ Ot Oz; E Ox; 0p;Ot" 
m x Rajal » 9H OF OF OH OF OF 
1 Ox; àp; i Opp; Oz; ax; ^ O2,0z; Opi Op; 
14 0 |OF(0F0H OFOH\| 1 4 89 |OF (OFOH  OFOH 
2 42, Op; (3x; \Ox; Op; pi Oz; 2 55 az; Op; \Ox; Op; p; Oa; 
Ee a a TA 
Ge 2 Ox; \ Op; Ox; Op; 2 ja Op; Ox; Ox; Op; 


Since by hypothesis 


aq o 
=0 and M= 
jan 0zj0p; i » No 


we have the equality 
z ôF ð 
z 0 rF 


We now establish that 


; n . 
/ / / 3 hon Op; Oar, i [b] ev?) Pdp'dzdt 


a j=l 


b 
1 2 ES 
22. J J J Y: 2. |F[ v;lë]dydzat, (2.19) 
2h j-1 Oz; 
a DR I= 
where i 
p P per SR t iyp! 3, 
F asy | Tere dp’. 
We have 


OF Oi 1 1 "V 
- E = iy(p'—p) / ! 
— On; dp; come | lo (s- 5 hy,t) $ (s t 3 hy,t)| e x F(z, p',t)dp'dy 
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z OF i PES ci 
= x Ji J [®]y; e" €^? F(z, p’, t)dp'dy 
j J R2?” 


j=1 


-È A J (dlv; F(z, y,t) e "dy 


whence by integration we obtain 


n AGE (ry J J (lysF e "rápdpázat 


=; 3 f 3 E n OPE p y t)[O]y;dydzdt, 


a DR” j=1 


where 


= 1 P 
F (et) = ros | Fte V? dp. 
R” 


We prove next that 


7 E b 
1 n OF* å 1 8 
JI Oz; Wydene ml | [X3 |E] yjl@ldydedt. — (2.20) 
Let 
F= sj wa | Fit dp = P(z,y,t) + iQ(z,y,t). 


R” 
Since the function F is real-valued, P and Q are even and odd in the variable y, respec- 
tively. We have 


j f f DE F[9]y;dydzdt 


a D R^ 31 


ay ES (Pe) si(a- PIE) tensis 


a DR? 
Since the function [6] is odd in y, the same is true of the function 


leis (3 ups ac). 


for all j, 1 € j € n, and we have 


As -Q — pêg 22) Ieluv =0. 


This establishes (2.20). 
We now show that 


Sue -fe (s- 599) - (2+ 5 hy,t)] 20. (2.21) 


j=l 
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Indeed, 
að (c—fhy,t) O0e(z4 E 


us. [9] = iu eege. às; 
We fix x, y and t, and adds the function 


n [O6(z—lzyt) 06(zlzyt 
so- £s ez Ee) morbum) 220 


Since g(0) = 0, we have g(z) = g’(6z)z,0 < 0 < 1, and so 


n 09? (z - 3 Ozy,t) 04? (z + ]0zy,t) 355 
Yi OzjÓz, ÓzjÓz, dd 222) 


Setting z = h we obtain from (2.22) 


09? (x — $ bhy, 09? (x +} Ohy,t 
=—h » prins t mn 


By condition 6 


n [89^ (r7 $t) aet (2 $9he] a 
k,jz1 Oar, Oat, y;Uk & V. 


Therefore, g(h) 2 0, which corresponds to (2.21). 

Since F = Fy — Fy, we have F|,-, = 0, Fh = 0, Fir = 0. 

Integrating (2.15) over z, p, and t, recalling (2.16)-(2.20) and using the specific prop- 
erties of the divergence terms and the Cartesian product representation D x R” x [a,b] 
of the domain, we obtain the fundamental identity 


dpdzdt 


y y [.£H ara, OH OF oF 
2 "Ozjür;Op;Óp;  Opjp,Oz,Ór; 


a DR” & 
1 ; 2 ð 
+z J f [Fl 2 3; (Plysdudedt 
a DR? j= 


b 
+ | J [ (Fi - F) (MF; - MF;) dzdpàt = 0. (2.23) 
2 DR” 
Since 


> 0, 


2 @H OFF, Y aFY z OH OF OF 
Opj;Op; Oz; -: Oz; 82,02; Op; Op; 


$j-1 t=] 44-21 


and (as proved above) ae 


j=l 


ES 2 0 (2.21) implies that 
Oz; 


j J J B lgrad, (F; — F))|? + (F, — Fe) (MF, - MF;)} àzdpdt <0, 
a D R” 
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which for F, # F, contradicts condition 11. Therefore we must have Fy = F. From (2.1) 
it follows in addition that \; = A2. The theorem is proved. 

Some conditions of the theorem are essential. We give an example of the nonuniqueness 
for the case when H(z, p) is specified. 

Assume that Q = {(z,p):z? +p? <1}, ie. Q is not a Cartesian product, and 
F(z,p) 2 1 — e1-6*t5), H = 4 (2? + p°). We have Flag = 0, (F, H} =0, F £0. 


REMARK 2.4. Under the condition 


A® = 4ro / F(z, p,t)dp 
R^ 


or with the information of type (2.5) available, it is possible to pose the question of 
simultaneously determining four functions (F, H, A, ©). 


2.7 THE EFFECT OF THE 'REDUNDANT' EQUATION 


In this section we restrict ourselves to the linear equation of transport because this effect 
may be present in more general kinetic equations including nonlinear ones. Let us consider 
the radiation transport equation 


10F — : ; r, 1 
x Taree P= Rn) )FdQ + =q, 


where Z = (x,y,z) € DC R?,t € Ri, |Q] = 1, D is a domain, B is a unit sphere centered 
at the origin of coordinates. Expanding the solution F(z,f2,t), the indicatrix f(z,Q), t) 
and the source function q(z, Q, t) in terms of spherical harmonics gives 


21 4 1(1- | m [)! 
4x. (I | m |)! 


1/2 
Yim(Q) = | | (—1)t/2t-m pil (cos 0) e", 


: t ; 1 1 
Q = (cos ysin 9, sinysin 8,cos#), 0€ < 2z, E «0x 2 


We can obtain an infinite system of equations for the expansion coefficients vij (z, t) 
of the solution F(z, p, t), see (Case and Zweifel, 1967). This system of equation runs as 


1 OYIm 1 1/2 199m 1/2 1/59 11m 
Lm *gi*1em (I 1 = m) TEE e (= m) (Le m) Em 


1/3 .a 
E (s. -i 5) [omae m - ^us 
-(I— m 4-2) (I — m + Dj) 
E 9 Q9 1/2 1/2 


(Lm 4 1) P m 4 29 met j 
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= Vino(chi—1)+—Qim [mS (2.24) 


Here fi(z,t) are the expansion coefficients of the scattering indicatrix f(z,t, QQ’), Qim 
are the expansion coefficients of the source function qg(z, Q, t). 

With respect to the inverse problem, the system of equations (2.24) has a very inter- 
esting property: it remains a closed system of evolution equations with respect to variable 
z while neglecting an infinite set of relations corresponding to the subscripts | =| m |. 
This property depends only on the term NV,F and allows us to determine, apart from the 
solution F(z, Q, t), various coefficients of the transport equation, if the proper boundary 
information is available, e.g. Fo = F|,-o. Furthermore, in the case of finite expansions the 
elimination of the sections and the scattering indicatrix from (2.24) leads to a sequence 
of linear problems. We shall restrict our consideration only to these systems of equations. 

Let 


Pim = Aim + iBim, Qim = Pim + iqim, 0< l< N, 0O<m< l, 
VN-4jm = 0, Qn4jm = 0, J 2 1, Bro = 0. 


Theorem 2.11. The closed system of equations (2.24) with respect to the real functions 


Qim, Bim Which does not contain the unknown coefficients (cf; — 1), l = 0,1,2,..., N, is 
as follows: 
Banm _ 1 ( Nm i. anm- | ÜBNm-i 
ðz  2\N—-l+m Ox Oy 
E eee * (Senmer | 98m 
2\N+1l+m Ox Oy j 
OBNn = 1 ( N+n i OBNn-1 ER Oann-1 
dz 2\N-14n Oz Oy 
1 ( N-n j5 OBNa+1 " OaNnn+1 
2\N+l4n Ox Oy : 
m —0,1,...,N, n=1,2,...,N 
Oo im au 
au = 


ðz ^ U mpl my?" 


2l + 1 Qai 2l +1 1/2 1/2901, 


v 


H 1/2 1/2 O01 1,m-1 OBi-1,m-1 
+5 (I+ m) (l+m-—1) (Sg +a 


1 Oo1 1m OB. 
—L(1— my? (1— m — 1)1/2 | mH OP emt 
^ (1 — m) ^*(I — m — 1) ( ER + By 


alg. 1/277 1/2 Oetisim—1 , biti m-1 
zC m42)^'(I—m-41) o EE MEE REN 


1 
-lq m4 124 m 4 2) (eto | Bins 
2 Ox Oy 


Qim E +100, 2141 


P Mic ENARRARE Ert en en ee 1/2 
(L— m) A(1 + my? v pu + (21 4- 1) 


t v ôt 


ða, 
Oz 
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| QU" -1)'? (aues 4 Hin) 


2 Ox Oy 
4 v2 aiya l-1 4 Oii 
2 Ox Oy 
| Qi 1)" (21 + 2)? (0a ii Bua 
2 Oz Oy l 
Of ia _ Bu 
Pua = (an) nya. 
21410 21+ 1 OBi+1,n 
E E A ne pe - (L£1- n) (I m 4- 1a bnn 
1 1/2 _ 1)1/2 OBi-1n-1 , 9O1-1,n-1 
+, (I+n) (Ln — 1) ("ze NO" E 
1 = 1/277 _ 4 9/2 OBt41n—1 Altin- 
^ (LI — n 1) ^* (I n- 2) (Mpa + E UD 


Oy 


OB. isa + 00141541 
Ox Oy 


OB, 20-41 1984 
qu + (21 + 1) Bz 


m — n) 1/2 OB ing O01 141 
;(- (l-n +1) (Sue y ES 


-5 (Ln 4 1) (L4 n 4 2) 72 ( 
Bin 
+ (ay aye lo TD 
| (20) (2l - 1? (88.114 4 Ooy~11-1 
2 Ox Oy 


v2 (Pee n 2 


2 Ox Oy 


NUR 1) 701 4-2)? (Bui " 90141141 
2 Oz Oy f 
l21,...,N. 


Assume that ain, n are the solutions to the above system of equations, then the sought 
coefficients o(cf; — 1), L = 0,1,2,..., N, may be calculated using ain, fi, by the formula: 


m dau 1 1/3904 
e(cfi — 1)ou = a AFI [eren E ug 
| (QI) "(21 = 1): Qoi i3 + OBi a 

2 Ox Oy 
į v2 ayi,- 2 OB. ia 
2 Ox Oy 


(21 + 1) (21 + 2)? [ðo | Bia 
rt or YJ Oe 
2 Oz Oy 
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As an example, we give the following result. 


Theorem 2.12. Let the coefficients aĝ (x,y, t), Bf(z,y,t) in the expansion of the given 
function Fo(z,y,t,2) = F|.-o be non-zero in a certain neighbourhood of the origin 
of (z,y,t)-coordinates. Then there exists a neighbourhood of the origin of (z,y,z,t)- 
coordinates such that the inverse problem of simultaneously determining the function Yim 
and o(cf; — 1), | m |< l has no more than one solution. 


We should emphasize once again that the effect of 'redundant' equations may be 
present in more general kinetic equations as well, since this effect depends only on the 
term vC) grad F. 


2.8 PROBLEM OF SEPARATION 


A solution of kinetic equations appears to have another interesting property: under a 
certain restriction on the source using the kinetic equation solution trace it is possible 
to represent the distribution function as a sum of two terms, one produced by the initial 
state and the other generated by the action of the source. 

We shall consider a kinetic equation of the form 


Of agrad, P o SUP + f, teR', eR”, zcR', (2.25) 


where St F is the collision integral, and f(z,v,t) is a certain function. The restriction 
concerning the right-hand side of (2.25) is as follows: the function q(z,v,t) = StF + f 
is independent of the absolute value | v | of the velocity v. In the linear theory this 
restriction is definitely satisfied for isotropic scattering or for Of /0 | ? |= 0. In nonlinear 
cases, for Of /ð | v |= 0, the function q is independent of | d | if the kernels of the operator 
St are independent of | v |. The general solution of equation (2.25) can be written as 


1 
F(z,9,t) = Fy(z — vt, v) + fats — àt(1 — £), v, t£)d£. 


Let 6 = af, | Q |= 1. Since q is not assumed to depend on a, then the following 
equality is valid: 


1 
F(z,of1,t) = Fy(z — oft, aN) +t J g(x — aMt(1 — £), 0, t£)d£. (2.26) 
0 
Let 2 and z be fixed, and for Q = No and xz = zo 
1 
a(a,t) 2 F, bot) Fo, c(a,t) = J qd£. 
0 


From (2.26) we obtain the equality 
a(a,t) = b(a,t) + c(a,t) (2.27) 


Here the function b is generated by the initial state Fo, and c by the function q with 
z = zo, Q = Qo. 
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Theorem 2.13. If the functions a, b, c are analytic in the neighbourhood of the origin 
of coordinates, then the function a in equality (2.27) uniquely determines the functions b 
and c. 


Pnoor. Bearing in mind possible generalizations, it is worthwhile performing the proof. 
It is clear from (2.26) that the function b depends on the variable combination (—at, o), 
while the function c depends on another combination, (—at,t), and thus we may write 
down the relation 
a(a,t) = b(—ot,o) + tc(—at, t). 

Here it is clear that in expanding the functions b and c in the Taylor series, the powers a 
of are greater than or equal to the powers of t for the function b, while, on the other hand, 
the powers of t are strictly greater than the powers of o for the function c. Therefore, b 
and c can be uniquely separated out of the expansion of function a, i.e. it is possible to 
split the distribution function. The theorem is proved. 


REMARK 2.5. On splitting the distribution function the inverse problems of determining 
the function q can be solved by the familiar techniques assuming the appropriate infor- 
mation is available. These problems will belong to the integral geometry and consist of 
an integration along lines. 

Of course, other results achieved in separating a solution are also of interest including 
those concerning more general kinetic equations. Note also that the quality of being 
analytic for the functions in Theorem 2.13 is essential: there are proper examples for the 
opposite case. 


2.9 DIFFERENTIAL AND INTEGRO-DIFFERENTIAL IDENTITIES 


In studying the uniqueness and existence of solutions of multidimensional inverse problems 
for kinetic equations the integro-differential relations are of substantial use (Anikonov, 
1985; 1986; Anikonov and Amirov, 1983). We shall give some examples of them here and 
some applications in Chapter 4. 


Theorem 2.14. If F(z,p,t) is a twice differentiable function and 


zd OF OF 
LF => "s FH BF- 
{ h x. ðr; Op;’ 


j=l 
then the identity 


OF BEP, e OF OLF - @H OF OF 9H OF OF 
LBF = ee, 
m 1 Op; “On; * 8s; pj t (ara Op; Op;  Op,Op; Oz; 2) 


is valid. 


The quantum equations of type (2.1) turn out to be similar relations which generalize the 
previously known ones. Here we restrict ourselves only to two equalities. Let F} and F, 
be two infinitely differentiable solutions of (2.1) with a — 0, which rapidly vanish at 0 
with respect to the variable p, such that F| = F; for t = a, t = 8, z € OD. 

Let 


F=F,-P, MR =Y 


j=l 


t F; = 
a Bp kK=1,2, 
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1 1 
[6] = lo (s Ms 5 wt) -6 (z+ TO) , 
2 2 
where F(z, y,t) is the Fourier transform of the function F(z, p,t) with respect to p. 
Theorem 2.15. The equalities 


J J [esa FP Fors Meas ff f 
D R” 


a DR” a 


ir 


a 


F 


! 32; (6]dydadt = 0, 
ji i 


lgrad, F? + F(MF, — MF) 


n 


oo h?" 
2 


BE psenzrogr| dzdpdt = 0 
m=0 (2m "E 1)! Jaj=m+1,|B|=m41 j 


hold. 


COROLLARY I. f M is a monotone operator and the given function 9 is convex with 
respect to the variable z, then F 0, zc D» pe R",a<t <b. 
The following identity containing the curvature tensor holds for the operator H in inte- 
grating over geodesic lines; 
A 

HF = Ẹ V;F 
where 
ð 
V; is the operator of covariant differentiation with respect to the metric tensor gij, T5, 
are the Christoffel symbols of the tensor gi; (see (Pestov, 1986)). 


h i pk 
Vi=V;- 3 


Theorem 2.16. The identity 
h |? "OCEANS Gi Apes 
lv F — Rat t* Vi FV'F + V; (VIF -HF) 
v h ACA E S h Gus 
=Ñ; [vir (VIF - evir) +2 Vi (FV'HF) -2FQHF 


h e el Th hk os ih 
holds, where R;;,; is the curvature tensor, Vi = Pao V'zg"V;,V'-g"V;,Q-V'Vi 
and F is a twice differentiable function. 


Theorem 2.17. If functions F(z,y,z), Q(z,y,z), d(z,y), c(z,y) belong to C?(D), 


D € R? and 
MF = (= - cos 6 + (E — ) sin 6, 
OF OF 
'F = | — — inĝ — | =— -d . 
M'F ( e) sing (& ) cosa 
Then the identity 


no Daw 00 | (OF ? (aF 2 
M'F3-MF - MF5-M'F = —- \(-4 «(2-3 
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ae ca OF DH 
"oy Ox Oz Ox \ Oy Oz Oz V Oy Ox 


holds. 


The use of this kind of identities for geometrical problems is exemplified in (Anikonov 
and Pestov, 1990). 


2.10 SOLUTION EXISTENCE PROBLEMS 


The results concerning the existence and uniqueness of the solution in three inverse prob- 
lems for the transport equation are given in this section. The general technique applied 
here is developed in (Anikonov, 1986; 1987; Anikonov and Bubnov, 1988; 1989) and in 
works cited there. By way of illustration, we restrict ourselves to a linear problem to 
determine a single coefficient. 

The boundary value problem for the transport equation 


ĉu + (w, Vzu) + o(z,w)u = e, J $5 w judu + h(t, 2,w)f(z,w), (2.28) 
a 


ulr_ = e, (2.29) 


is considered, where t € R, z € D C R^, D is a strictly convex bounded domain with 
boundary y € C?, w € Q, Q is a sphere of unit radius in R” and T. = {(z,t,w) : t € R!, 
z € y, (w,n) < 0, n being the unit vector of the outer normal to y}. 

Three inverse problems are posed: 


1. Find functions u(z,t,w) and o(z,w). 
2. Find functions u(z,t,w) and c,;(z,«). 


3. Find functions u(z,t,w) and f(z,w), in (2.28) if conditions (2.28) and (2.29) are sup- 
plemented by data for t = 0 and the stipulation that the solution u(z,t,w) behave 
regularly with respect to the variable t on (—oo, oo), i.e. by the following conditions: 


u m uo, uo h= p lizo; J- = {(x,w) Ir € Y (n,w) < 0}, (2.30) 
max fot LA D? [la + [Veal] dà «oo, p22, (2.31) 


where ù is the Fourier transform of the function u(z,t,w) with respect to the variable 
t. 


That the inverse problems are well-posed is guaranteed by the following assertion. 


Theorem 2.18. Let c, o, € C'(D) x C(Q), g € C and let the condition 


fita] 


2 


. 6d? (i + løslo 


50 Chapter 2 
2 


Ef) = max. f a 13 D" [If + Vei] da 


+ (ra. ag IV zc. llo 


"Inr 
R 


be satisfied. Then, for any functions hf and q such that 


le; llo 
0 


J |Vg| dw 
2 


max f (1+ 11” [fif + [Veh s[] aa < os, p22 
there exists a unique solution u(z,t,w) of problems (2.28), (2.29), and 
max f (1+ 1A D” [lèl Iva] dà < cB hf), p22, 


—O00 


|®(z,w) lp = max | #(z,w) |, d= diam D. 


Boundary value problems 1-3 posed above turn out to reduce to an investigation of the 
boundary value problems for integrodifferential equations (see (Anikonov and Bubnov, 
1988)). This reduction is formulated in the following lemmas. 


Lemma 2.2. If v(z,,w) is a solution of a boundary value problem 


—idv + (w, Vev) +v e £ J val = s, | godu + fh, (2.32) 
ue a Q 
hase. Hil Ses: eR” (2.33) 
of the class - 
max f a 1A D [lv P +V] dà «oo, p22, 


then functions u(z,t,w) and o(z,w) sought in inverse problem 1 are computed from the 


formulas 
oo 


o(z,w) = Rh t B J Avd\, u= / ve Mq), 
o o 


—00 


Fo = 0,(z,w) ET — h(0, z, cw) f (z,w) — (w, V suo). 
a 


Lemma 2.3. Ifv(z,A,w) is a solution of a boundary value problem 


i 


J guodw’ 
a 


—iÀv + (w, Viv) + ov = |F — Í AvdA E + fh, (2.34) 
-co n 


2020, (2.35) 


vh- -$, 


J guodu” 
n 
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of the class 


foime fecit ioter em zn 


then functions u(z, t,w) and o(z,w) sought in inverse problem 2 are computed from the 
formulas 


TR m J Noddy: u= / veg), 
f guodu! / E 
2 


(w, Vv) + ouo — f(x, w)h(0, z,w) 


Fy = 
J guodu” 
Q 


Lemma 2.4. If v(z,\,w) is a solution of a boundary value problem 


—idv + (w, Vv) + ov 


= 0 | guoda/ + wees) X i AvdA + (w, Vato) + oto — o, J guodw'|, (2.36) 
k(0,z J J 


v|.-60, {h(0,z,w)| > 6 >0, (2.37) 
of the class 
foina J aD [Lo P+ 1 ex em p22}, 


then functions u(z,t,w) and o(z,w) sought in inverse problem 3 are computed from the 
formulas 


1 LUE 
f(z,w) = (0, 2,0) le Vzuo) + ouo — s. [ stai Ei | ; 
u(z,t,w) = f ex 


The solvability of problems (2.32)-(2.37) is established in the following theorems where 
the condition of the smallness of the data has a somewhat cumbersome form. 


Theorem 2.19. Leto, € C'(D)xC(Q), g € C'(D)xC(0)xC(Q), uo € C*(D)xC(Q), 
| uo 12 6, h(0,2,5)f € CD) x C(Q), E(Q, hf) « oo, and p 2 2 and let 


LA 


ane | [Fol Ista, IV-£ilo 


Bp 3) Pp 3) V uolo 


ne (Ine Fol? + lo IE (2P2R] | 4(12)?d? R? 
Fsllo 


WeFolo | tot 


falas 
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+ lieso 


2 
a 


Ry = max Jos |)” [IP + 1V1] dd + max / (1+ 1 pe [ig + 


Uu 


ZI 
R 


~00 


p22. 
Then there exists a unique solution of problems (2.32), (2.33) and 


max J (1+ 1A D? [l6? + vo] a < cB(¢, hf). 


=00 


Theorem 2.20. Leto € C}(D)xC(Q), g € C(D)xC(Q0)xC(Q), uo € C?(D)x C(Q), 
Sguodu'| 2.6 » 0, h(0,2,w)f € CHD) x C(Q), and E(Q, f) < oo. Let Q and hf be 
(9) 


real function and let 
2 


; , OVER 
: 6?(2p — 3) 


nr 
2 


12d? (ier Mo f 1g | de! 
2 


0 


2 
eme [Ics Iv fisia] HM "4 | Vag | dw | 

R 0 0 

2(12)2d? R? j ; 
u 
< 
oH fiala w| |f Tema] + oa ETE <1, 
0 


p22. 


Then there exists a unique solution of problems (2.34), (2.35) and 


max J (+ LA D?" [iô]? +1V.6)"] da < cB(¢, hf). 


Theorem 2.21. Let o,o, € C(D) x C(Q), g € C (D) x C(Q) x C(Q), uo € C?(D) x 
C(Q), h(0,z,w) #0, E(Q, hf) « oo and p> 2. Let 
1 

fisia + 
2 0 


os 22 1/21? 
f 413 Dy» [| dA 
h(0, z,w) 2p — 3 


6d? |lelo + leslo 


+ les llo 


Jig lau" 
2 


VAZIA 
Q 


0 


* (ia. + I Vzoc;l 
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Then there exists a unique solution of problem (2.36), (2.37) and 


max f (1 A p [lol - Iv.of] dà < cE( hf), p22 


100 


The proof of these theorem is based on the study of nonlinear integral equations equivalent 
to the corresponding boundary value problem. Note that the conditions of Theorems 2.18- 
2.21 are satisfied with a sufficiently small domain D or with input data sufficiently small 
in norm. 


2.11 AN INVERSE PROBLEM OF MATHEMATICAL BIOLOGY 


At present time the theory of differential equations is used for the quantitative and qual- 
itative description of biological processes, for example, the morphogenesis or evolution 
of population (see (Prigogine, 1980; Murray, 1977; Svirizhev and Logofet, 1978; Pilant 
and Rundell, 1986; 3Volterra, 1931; Kolmogorov, 1985; Kolmogorov et al. , 1985; Lecture 
Notes in Biomathematics, Springer Verlag, Berlin (1980); Keyfitz, 1968; Rundell, 1989)). 
Here, often evolution equations of the form 


ĝu = Au + F(z,t,u), (2.38) 
arize, where the vector u(x,t) characterizes quantitatively the considered biological sub- 
stance depending on time £ and a classification variable z, and the operator A and the 
vector F determine the propagation and interactions of elements of the substance, respec- 
tively. Pertaining to type (2.38) are the Lotka- Volterra systems of equations in the theory 
of evolution of populations, the Turing system of equations describing chemical processes 
in cells, and the Kolmogorov-Petrovskii-Piskunov equation in the theory of propagation 
of genes. 

For equations of type (2.38) usually initial boundary-value problems are considered 
and questions of a qualitative aspect with the present biological interpretation. 

Presently, inverse problems of mathematical biology are discussed, consisting in defin- 
ing not only solutions of the considered equations , but also certain coefficients entering 
into these equations (see (Rundell, 1989)). At the present time new statements are made 
and systematic studies are undertaken of inverse problems of mathematical biology, where, 
besides the solution, the sources, as well as the coefficients of diffusion, transport, absorp- 
tion, birth, death, etc. are described. Here, as in the inverse problems of mathematical 
physics, the information must be the biological measurements and constraints connected 
with the initial boundary value data and qualitative behaviour of the solutions. 

In this section an inverse problem of determining the absorption coefficient is studied 
in the equation which describes the evolution of populations, taking into account memory: 
find two functions u(z, v, t), A(z,t) in the domain t € R, v € R, z > 0 such that the 
function A(z,£) does not depend on v. It can be a generalized function of the variable, 
and m 

x t a, — Ae J A(z,t — r)u(z,v,T)dT = 0, (2.39) 


-00 


u |z=z,= a(v, t), (2.40) 
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where (x) > 0, z Æ zo, (zo) = 0, a(v,t) are given function, and zo is a fixed number. 

The biological interpretation is the following: t is time, z is the physiological growth of 
the specimen, v is the velocity of its change, u(z, v, t) is the number of specimens, A(z, t) 
is a coefficient characterizing birth and death, the presence of the convolution in (2.38) 
characterizes the consideration of memory (Kolmogorov, 1985), ®(x) is the potential. 

Without accounting for physiological growth and for ®(z) = 0, v = 1 the inverse 
problems for the equation u; + uz + Au, considered in (Rundell, 1989). Let us also observe 
that in (Rundell, 1989) the desired function \(z,t) can depend on the solution u, for 
example through the formula À = F(z,t, f udv), where F is some function. 

R 


EXAMPLE. Let the given function a(v, t) in the inverse problems (2.39), (2.40) be repre- 
sented in the form 


i -1 
a(v, t) = -e |-; (r v-- an) | ; 
fe?) J g (50? - gn?) dn ie ü oss G ;v) 


where f(y) > 0, y > 0, g(y) > 0 are differentiable functions. Then the functions u(z, v, t), 
A(z, t) determined by the formulas 


1 
f(v? +20(z)) + J g ($0? - ly +8) dn 


x exp E (e + 9(z)) JL (s v? -5n soja) |, 


A(z,t) = -38 ts(9(2)77 ze ) 
satisfy (2.39), (2.40). Observe that the Fourier dr of the functions u, À have the 


form 
ZR je ORE G v? Z EXE jj]. 


^ 


A = ~iw — 29(®(2)) Seu? 


u(z,v,t) = 


Let us consider at first the inverse problems (2.39), (2.40) in the spectral statement: 
find functions w(z,v,w), A(z,w) in the domain z > 0, v € R, w € R such that 


Ow ðw 06 . 
al Bo De + (A(z, w) + iw) =0, (2.41) 
v [24 G(v,w). (2.42) 


Theorem 2.22. Let the function G(v,w) be greater than zero and twice continuously 
differentiable with respect to the variable v, à? (0,0) = 0. Then 


w(x, v,w) = A(v? + 2(z),w) exp fe G v? — PU + TOR ar ; (2.43) 
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À(z,u) = 2B((2),u) 5 — iw, (2.44) 


A(y,w) = y á( /y,w)à(— vyw), y2 0, (2.45) 


B(y,w) => n F In ae (2.46) 


aed 


PnRoor. For the proof of the theorem it is necessary to establish that the functions w, À 
defined by (2.43), (2.44) satisfy (2.41) and (2.42). Let us note that since the given function 
à, à, (0,w) = 0 is twice continuously differentiable with respect to v, the functions A(y,w), 
B(y,w) for each fixed w are at least once differentiable with respect to y. We have the 
relations 


Ow i 7 1 1 
2 2A; (v? + 26(z), v) exp j B (; ye z” + 9(2),w) d 


+A(v? + 20(2),w) exp {/ B E v — PU t a(z), w) a 


—U 


[n v? -5n ? 4 (zx), w) dn] 2È ; (2.47) 


—vV 


= a (v? + 26(z),w) exp [j B G v? — 7 + TOR J 
+A (v? + 26(z),w) exp {/ B G v? — Pv t a(z), w) 7 


—U 


JEG v? -5n 24 (2), w) an} 


+2A (v? + 26(z),u) B (6,0) exp j B G v? — Ly + OR an) 2 (2.48) 


Multiplying (2.47) by v, (2.48) by (— 00/ Oz) and adding the expressions obtained, we 


arrive at the equation 


að f /1 1 
ar mir -2B(8, w) 5 Al? + 2®,w) exp t B (5 v - "id t a(z), o) «| , 


=v 
which can be rewritten in the form 


Ow _ 9w06 


which in fact completes the proof of the first part of the theorem. 
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Now let us show that w |;-4,— G(v,w). By hypothesis ®(z9) = 0. Therefore, from 
(2.42) we have 


v |e=2,= A(v?,w) exp [j B G y= T ^u) 7 : (2.49) 


At first let v 2 0. From (2.49) we obtain the equation 


mm putari ) 
w |; A(v wwf] aG Q9 dg). 


Setting lo - lp — y we find the relation 


w |.-2,7 A(v?,w) exp we] (a. s m dy 


ify. Neu. âl, w) 2 
Aaaa B= E à(- 25, w) Vy — 


then by Abel's formula of inversion of an integral equation 


w |[;-4,— yà(v,w)à(—v,w) exp I» em = ü(v,w). 


Since 


w 


Consequently, for v > 0 the equation w |z=z,= â(v,w) is established. 
Let us prove the equality w |,-2,= 4(v,w) for negative v. From (2.49) we have the 
relation 


m 


op 


1 


w(zo, —v,w) = A(v?,w) exp 4 —2 P B (5 v? — jT») dn?, v20. 
0 


2 


Repeating the previous argument, we obtain the required equation 
w(ro,—v,w) = à(—v,w), v0. 


The theorem is proved. 

Let us return to the inverse problem (2.39), (2.40). Besides the restrictions presented 
above on the variable v for the existence of the Fourier transforms of the functions w, and 
X restrictions are necessary with respect to the variable w of the information â(v, w), 


LF ; 
a(v,w) = z- COD 


From the representation of the given function à in terms of the function A(y,w), B(y,w) 


in the form " 
a 2 1, 1l, 
à — A(v DESEE -3" an} (2.50) 
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and formulas (2.44), (2.45) it is clear that there is a unique correspondence between 
the functions à > 0 and the functions A(y,w), and B(y,w), y > 0, characterizing the 
even and odd parts of the functions @ with respect to v. We will, in correspondence 
with this remark, assume that the information à is represented in the form (2.50), where 
A(y,w), and B{y,w) satisfy the following conditions: A, B are continuously differentiable 
with respect to y, the function A along with its derivative with respect to y is rapidly 
decreasing in w, and the function B along with its derivative with respect to y is bounded. 
Under these restrictions we have formulas for the functions u(z, v, t), A(z, t) in the inverse 
problems (2.39), (2.40): 


oo 


u= J je + 20(z),w) exp [j B G v? — 57 + OR 2 elt dus, 


—oo —v 


oo 


A(z,t) = J [rocno -i tà 


—O00 


CHAPTER 3 


Geometry of Convex Surfaces in the Large and Inverse 
Problems of Scattering Theory 


3.1 GEOMETRICAL QUESTION OF SCATTERING THEORY 


In applications we often meet situations, where a source of a physical field or the totality 
of a physical field is a certain manifold. Physical fields can be electromagnetic, acoustic, 
seismic, thermal, gravitational, etc. It appears natural to name such manifold-emitting 
sources manifolds. Earthquake sources, stars are examples of emitting objects. The 
capability of emission may have various causes: it may be an internal property of the 
emitting object, or the manifold becomes emitting after external action, or both cases are 
taking place at the same time. The main problem in this case consists in the reconstruction 
of the emitting manifold by the field, which it creates. This chapter is devoted to these 
questions. Our main attention is devoted to the inverse problem of scattering theory, that 
is the reconstruction of the scatterer by the scattering information. Nevertheless, many 
results, we believe, may be useful in other situations. 

The physical formulation of a variety of inverse problems consists of the following: 
there is a certain body-scatterer in space, which is irradiated by an electromagnetic, 
acoustic or another field. 

The body scatters the field: the form of the scatterer can be found by scattering in- 
formation. In physical problems a body often is irradiated from different directions at 
one or several frequencies in such a way that extensive information on the scatterer is 
obtained from scattering. The amplitude of the scattered field is a physically authentic 
information. It is very essential from a mathematical viewpoint, how and where is the 
field measured, is the source monochromatic or not, what is the nature of the direct wave, 
etc. Since scattering problems are very complex, asymptotic methods are widely used at 
investigations and applications. Frequencies, and the distance from the scatterer of the 
measurement point are usually parameters, by which the asymptotic expansion is pro- 
duced. In this case geometrical characteristics of the scatterer such as the square of the 
orthogonal projection on a plane, Gaussian curvature of boundary, characteristic func- 
tions, support and other functions arise in the asymptotic expansion as coefficients. This 
geometrical information is very important, because in certain cases it permits effectively 
the determination of the form of the scattering body, i.e. to solve the inverse scattering 
problem. In this connection it is especially essential that methods of geometry in the 
whole and methods of integral geometry may be used with success. They lead to new, 
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and sometimes exhaustive results. It needs to be underlined once again that geometrical 
methods operate asymptotic. Therefore theorems of existence, uniqueness and especially 
solution stability of inverse scattering problem are important here. 

Consider a first example, showing the geometrical characteristics in a scattering prob- 
lem. Many problems of wave propagation lead to the boundary problem for Helmholtz 
equation. So if a plane wave uo(z,p, A) = exp(iA < z,p >) (A being the frequency), 
spreading in the direction p, |p| = 1, falls onto scatterer K , bordered by a closed surface 
B, then for the scattering field u(z, p, 4) we have the external problem 


Au+\u=0, zc€R"VK, 


T(u(z, p, A) + exp(il <2,p »)ls = 0, 


zoo 


En Qu . ni 
u =0 (\x|)" yn zizi xd iuz} = o (|z) yn 


where T is the boundary operator of Dirichlet, Neumann or impedance condition. At 
m = 3 an exact solution of the direct problem (i.e. finding the scattering field u(z, p, A) 
from known boundary B of scatterer K ) is given by the Kirchhoff formula 


1 ð eiAle-yl OQu(y,p A) ePlz-ul 
, „À = eal , ,A a re RUE oe , 
u(z P ) 4r J {uy p in, |z — y] Ony |z — yl Sy 
(ny being outer normal in point y € B). 
In the case of a more general equation the solution u(z,p,A) is computed by the 
formula 


1 ô u(y, p, À 
vini)» f funo Glena - P Bua 8 
B y v 


where G(z, y, A) is the fundamental solution. In this case solution u(z,p, à) may be a 
vector, G may be a matrix. In particular, in the case of the dynamical equation of 
elasticity theory with constant coefficients elements G; of matrix G are given by relations 


1 eikzlz-yl 1 a? eiilz-yl _ eikalz-yl 
Ge =a thi a Ss a a ee 
2zu " |e —y| 2rpM Ox,0x; iz — yl 


A? A? 
£ , B= adit a and p being Lame constants. 
at 2p H 
The inverse scattering problem consists of determining the boundary of a simply con- 
nected scatterer K by information on the scattered field u(x, p, à), to be known to a 
certain set M of the variables z, p, À, i.e. by function u(z, p, à)|m. In applications usually 
a known function 


where k? = 


Alp A= lim: Iz e ^ Flu(lz| v, p, A) 


is considered |v| = 1, |p| = 1, A > 0, which is called the scattering amplitude or the 
intensity of the scattering field |A| in the direction v of observation. By the Kirchhoff 
formula we have 


= 1l ~iA<v,y> |: 3 
A= ds J ià < ny, v > u(y, p, A) 95, ^ P A)| ds,. 
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Consider the Dirichlet problem, and the scatterer smooth and strictly convex. Apply- 
ing Kirchhoff's approach (see (Taylor, 1981)), for the scattering amplitude we obtain the 
following expression 


" 
A(v,p,À) = eal posme (€ ny,v » +|< ny, p »p] ds,. 
B 


In particular, the scattering amplitude in the forward direction v — p, is proportional 
to the area F'(p) of the shadow of scatterer K ; 


id ir 
A(p, p, A) = on J < p, ny > ds, = 5; (P. 
«pny»20 
In case of v = —p (backscattering) the formula 


iA ; 
A(—p, p, à) = E J eT»? < p ny > dsy = p(p, À) 
<Pmny> <0 


holds. 

Notice that if function p(p, ) is known for all directions p and frequencies \, then 
computing p(p,A) + p*(—p,) (* denotes conjugation) and applying the Ostrogradskii 
formula we obtain the relation (Levis, 1969) at y(y)=1,y EK; 4Y(y)) 0 y € K, 


oo oo oo 


„PAP A) ete pA d= f f [ue9a- (A, p), 


—oo —oo —00 


i.e. (A, p) is the Fourier transformation of the characteristic function y(y) of scatterer 
K. In practice measurements of p(p, À) are produced in a limited range of frequencies 
land directions. 

Another way to solve the scattering problem uses geometrical optics. In this case the 
asymptotic expansion of the scattering amplitude at A — oo contains a summand, with 
a Gaussian curvature of border of the smooth strictly convex scatterer as a main term of 
expansion. So in (Majda and Tayler, 1977a) it has been proved that for v # p, 


fim Ap Jf = y(r p) KG), 


where y is reflection coefficient, K (7) the Gaussian curvature of the border of the scatterer 
in point ğ with normal n = (v — p)/ (|v — pl). 

It is possible to come to geometrical characteristics of the scatterer according to (Ma- 
jda, 1976a) in the T way: consider the mixed problem for hyperbolic equation 


c Qu 
2 = oer az; s [estne =| >» TEK, 
i j= 
ðu 
t|... 
where ajj(r) are the elements of a smooth positive symmetrical matrix with aj; = ôi; 
for |z| > p (domain K is contained in a ball |z| < p). Let u(x,t) be a solution of this 
problem. Let 


ult=o = filz), 


= fa(z), Tu|-x = 0, 


K*(s, p)- im, tu, [(s + t)p.]- 
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Functions K+(s,p) characterize the behaviour of the solution along ray z = (t + s)p at 
large positive and negative time. Consider the scattering operator S, which is defined as 
mapping S: K~(s,p) — K+(s,p) and permit representation (Majda and Tayler, 1977a) 


K*(s,p) = f fsi — 5,p,v)K~ (5, v)dšdv, 


where kernel S(s, p, v), s € R!, |p| = |n| = 1 is a distribution in D'(R! x 5? x S?). It 
turns out that for kernel S(s, p, v) the asymptotical expansion 


S(s, p, p) = 2r (= Ko") 8 (s + 2h(p) + i&(s + 2h(p)) 


+ (smooth members), 


holds for p from a certain open subset (dense everywhere) of unit sphere S?. In this 
formula y; are points on OK, in which < y;,p >= h(p), h(p) is the support function of 
convex hull of surface 0K; K(y;) is the Gaussian curvature in point y;. Moreover, it turns 
out that 

A(p) = max sing supp S(s, —p, p) = —2h(p), 


i.e. we can obtain the support mA of the convex hull of the scatterer from information 
about scattering operator. 

At last, we provide the result, connected with the convex hull of the scatterer. Let 
$(z) be a finite function of real variable z € R™ and let f(z), z € C™ be determined by 


the equality 
p xe : 
mc $ -iezz»gq . 
fe) = C) [seas 


The scattering problem leads to a similar integral, while, as noted by, for example, 
(Levis, 1969), the support of function $(z) plays the role of the scatterer. Therefore, 
defining (x) becomes an important problem. By the Plancherel-Polya theorem the 
support function of a convex hull of support of the function ®(z) is equal to 


h(p) = sup lim gin lf(e+iRp)]. 
zERm Roo 

In real situations and when investigating inverse problems frequency A, as a rule, is 
fixed and there is only a limited range of directions, for which the radiation intensity is 
measured. Therefore it appears expedient to pick out principal geometrical information, 
connected with the radiation intensity, and to put this outlined component into the base 
of the investigation technique of inverse problems to recover the emitting manifold. 

By formula (3.1) analysing of various investigation methods of problems, in particular 
scattering, leads to the necessity of study of problems connected with the representation 
of the amplitude of the field in the far zone in the form 


A(n, A,t) = ps eO) o), y, n, 1)ds,, (3.2) 
Bi(n,t) 


where B(n,t) is the illuminated part of manifold B(t) in the direction of n, f and j 
are certain functions, t is time, À is the frequency, and n is the direction of reception of 
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emission. In the case of a regular convex surface B(t) formula (3.2) may be rewritten in 
the following way 
A(n, à, t) = J ef v(m Gd, y(p, t), n, t) Ri Rady. (3.3) 
<n,p>>0 


Here y(p, t) is the focal radius of surface B(t) in the point with normal p, |p| = 1, Ri (p,t), 
R2(p,t) are main curvature radii of surface B(t) as functions of normal p and time t. By 
solution of the Minkowski problem it is possible to find, that y — Mz, where M is an 
operator, giving z = R,(p,t)R2(p,t) corresponding to point y(p,t) of surface B. That is 
why (3.3) can be rewritten as 


A(n,A,t) = J eif Mn) 30. Mz, n, t)z(p, t)dwp. (3.4) 
<n,p>>0 
As it was mentioned above, the really measured information is |A|. That is why, and by 
the above in the stationary case when 
AA)e f e) 950, Mz, n)z(p)dw 
<n,p>>0 


from which we separate the main geometrical information, connected in with the scattering 
intensity. Namely, for |A| we use the representation of the solution of the inverse problem 


Al = af |< n, p >| Ri Fadu, + bf x(< n, p >) Ri Radu, + Q. (3.5) 


Here w is a unit sphere in R’; f|« n,p >| Ri R;du, is the area of orthogonal pro- 
w 


jection of the concave surface to plane, orthogonal to n; x is the Heaviside function; 
f x(« n, p >)RıRıdw is the area of the illuminated part of the concave surface in the 


direction of n; a and b are constants; Q is the rest. Relation (3.5) is taken as the basis 
for further investigations. In this case if P^! is the inverse operator 


Pz= afi< n, p >| 2(p)dw, + JEG n, p >)z(p)dwp, 
then from (3.5) we obtain an equation of the second kind with respect to the sought 
function z(p) = R Rz, which has the form 
2(p) = A(p) + Tz, 


here 


Á(p =P JA], Tz = P'Qz, 


and to which standard methods of the exploration technique may be applied, founded, 
for instance, in the theory of fixed points. In this connection principal attention will be 
given to operator P: 


Pz= «fic n,p »|z(p)du, + b f x(< n, p >)z(p)dwp. 
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3.2 INTEGRAL EQUATION OF THE FIRST KIND 


In this section integral equations of the first kind are considered which are the result 
of inverse problems of scattering theory. In particular, geometrical aspects to define the 
convex surface by the functional from its orthogonal projections (shadows) and illuminated 
parts. 

Induce notation is used in this section. Let n,p be points of the unit sphere w = 
z € R2, |z| =1. For function f(n), given at w, let f(y,n) be denoted as its average value 
on the circumference y = const , when point n is a pole of the spherical coordinate system 


Jen) = z- f fæ, 7). 


Here g,t are spherical coordinates of the point p € w with respect to pole n. As f* are 
the even and odd parts of the function f(n), that is f* — 4 [f(n) + f(—n)], A is the 
Laplace-Beltrami operator on w, <> means a scalar product, w(n) = p Ew, <n, p>>0 
is a hemisphere, depending on n € w, p(n, p) = a |« n, p >|+ bx(« n, p >), x is Heaviside 
function, a,b are constants. 

Consider the first kind equation with respect to measures p on w. 


f(n) = [ P(n,p)u(dup), n Ew. (3.6) 


Theorem 3.1. Ifa £ 0, b 40,a+6 4 0 then equation (3.6) has not more than one 
solution p. 


PROOF. Let p be a solution of equation (3.6). Show, that it is unique. Since at any p € w 
the function |< n, p >| is even in the variable n, then by (3.6) we obtain equations 


b b 
ft(n) =a / I< n.p >| u*(durp) + 5 u* (dwp) + ; i u*(dep), (8:7) 


Fin) = 6 | n (dup), (3.8) 
w(n) 
where l'(n) 2 p € w, «n, p >= 0 is circumference, and pt, u^ are the even and odd part 
of measure p, respectively. That is ut = 4 [u(Q) + u(—Q)] for any set of Q on w, —Q is 
the set diametrically opposite to Q. Show that equation (3.7) identifies uniquely with the 
even part ut of measure, and equation (3.8) identifies with the odd part u^ of measure 
n. 
Consider integral 
J pt(dwp), nEw. 
T(n) 


Let u = f + faing be an Lebesgue expansion of measure u to the amount of being 
absolutely continuous and singular with respect to Lebesgue measure m on w. Then 


PICS 
P(n) 
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Support of the singular measure has zero Lebesgue measure, therefore u*(T(n)) = 0 
on the set of zero m-measure. Integrating both parts by sphere w and noting that 


/ |< n, p »| do, = 27, Tis (F(n))du, = 0, 
we obtain the equality 


f f* (n)du, = 2x(a + b) J ut (dwp). 


w 


With this we can rewrite the equation as 
+ b 4 
fi(n) =a f I< np >| ut (dwp) + hng (T(n). (3.9) 
From hereon for any function f(n) is designated by fi(n): 
-lp i + 
hin) = f*)- aaah | f*(n)dos. 


Furthermore, by proof of (Aleksandrov, 1937a,b; 1938a,b) and with the assumption that 
Hsing (F(n)) = 0 almost everywhere at w, we obtain that equation (3.3) uniquely identifies 
with the even part of the measure u. For the proof of the uniqueness of the solution of 
equation (3.8), obviously, it is sufficiently to show, that from the equality 


J W~ (dwp) =0, nEw, 
w(n) 


follows u^ = 0. Draw two arbitrary plane through a straight line, passing through the 
origin and belonging to a plane. Let these be planes P(ni) and P(n;) 


P(n)-z€R?, |z| 21, <2,n;>=0, i=1,2. 


These planes divide sphere w to four parts. One of these part (no matter which) we denote 
as Q. From the equality 


|W (dey) f (du) =0 


w(n1) w(nz) 


because of the oddness of u^ (y^ (—Q) = —u-(Q)) it follows that 
f c» - 0. 
Q 


We restrict those Q, which lie in the half-space z3 € 0. Let w~ be a hemisphere which 
belongs to this half-space. Reflect it reciprocally unique to plane z3 = —1 projecting from 
the origin. In this mapping region Q at sphere w defines a rectilinear strip T(Q) at plane 
23 = —1, and changing Q one can obtain any strip. The integral 


J J(z1, x2) (dzi, dz2), 
T(Q) 
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corresponds to the integral 


RC 
Q 
where J is the Jacobian of mapping by projecting. That is why from the equality 


J 1 (dep) =0 
Q 


follows the equality 
pe (Jdzi,dz;) = 0 


a<<zrv><p 


for any vector v = (vi, v2) and for any o and £. 

In this way, the single-valued solution of the equation brings about the uniqueness of 
the reconstruction of the measure by its meaning at half-spaces. A positive answer to this 
question is given in an article by (Khachaturov, 1954). The theorem is proved. 


REMARK I. f a Æ 0, b = 0, or a + b = 0, then from equation (3.6) only the even part u* 
of measure can be determined uniquely, and by b # 0, a = 0 the odd part. 


Let in equation (3.6) measure u be absolutely continuous, and z(p) be density of 
measure u. Consider the following equation with respect to z(p): 


f(n) = af |< n, p »|z(p)dup + b | z(pjdwp = Az. (3.10) 
By Theorem 3.1 it follows, that equation (3.10) has no more than one solution z(p) € 
Li(w). 
Theorem 3.2. The eigenvalues A, of equation (3.10) are 


2n (a +d), r — 0, 


2k — 3|! 
zs 4ra(—1)* 2E FO r 22k, k=1,2,..., 
= N 
Prof- BEIDE, r-2bRl k-01... 


The eigenspace, corresponding to eigenvalue \,, consists of all spherical harmonics of 
order r. 


PROOF. Take a set of spherical harmonics Y,(p) of order r in equation (3.10) and make 
sure that it satisfies the equality AY = A+, Y, where 


AY,(n) = [ T(n, p)Y;(p)dwp- 
Let r = 0, then Yo = const and 


AYo = af |< n, p >| Yodoy +b J Yodwp = 27r (a + b)Ys. 
w w(n) 
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If r = 2k, k = 1,2..., then we have (see (Blaschke, 1916)): 


b 
a f |< n,p >| Ya(p)dwp + 5 f'Ya(pdes = «fI mp 7lYa(p)do, 


([2k — 3|)!! 


= tral- Ea OMT 


By r=2k+1,k =0,1... we obtain 
AYan =a / «n.p >1Yanss(p)dwp +b f Yansi(p)dup. 
w(n) 


Since Y2x41(p) is an odd function at w, the first addend here is equal to zero. For cal- 
culating the second we accept point n as a pole of the spherical coordinate system y, t, 
and 


AY = b f Yous(p)dey 
w(n) 


aai 7 i 
EU s [oe [P (cos 7) [ogee cos mT + 63,4, sin mr| sin ydy 


m=0 


> ord [Pager in 


2k — 1|)! 
= 2rb(— 1) AM Gens. 


where pim are the associate Legendre polynomials, P, = pp, Here we seized on the 
equality of (Gradshtein and Ryzhik, 1980) 


[2k — 1|)! 


| tantoa = =(- DA t2)! 


Show now that a9, 41 = Yor41. For the average of the function Y2,41 on the circumference 
y = const we have 


Qn 

n 1 

Yaai(n) = z; | Yanl r)dr 
0 


Qn 
= = X / PE (cos y) D cos MT + b, ,, sin mr] dr 
= P2rx41(cos 1) 0241 
Supposing in this equality g = 0 and considering that 
Yos (y, n)lg-o = Yorgi) Pox4i(1) =1, 


we obtain a2,,, = Yogi. 
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In this way the spherical harmonics Y, are eigenfunctions of operator A. As the spher- 
ical harmonics form a basis in L2(w), the equation does not have any other eigenfunction. 
The theorem is proved. 


Theorem 3.3. If in equation (3.10) function z(p) belongs to class C*(w), then function 
f(n) belongs to class C**!(u), k = 0,1.... Equation (3.10) is analytical on sphere w 
solution z(p) if and only if function f(n) is analytical at w. 


PROOF. Continue function f(n) given on w, to all values z, z # 0 supposing 


f(z)- )cefKk z,p >|z(p)dwp + b J z(p)dwp, z € R?. 


(x,p)>0 
From à 
US J piz(p)dwp + E piz (pdsp, ¿i= 1,2,3, 
Oa; |e 
(z,p)>0 T(z) 


where l'(z) = (p € w,< z, p >} is the circumference at w, it follows, that f(r) belongs 
to class C! (R? V 0), if function z(p) is continuous at w. Let z(p) € C'(w). Using (see 
(Blaschke, 1916; 1930)) 


u(p)dsp= f < Ar(piz(),p) — 2p:2(p)p, 2 > dwp (3.11) 
T(z) w(n) 
following from Green’s formula it is possible to write 
l 
J piz(p)dsp = l J < Ai(piz(p), p) — 2piz(p)p, z > dwp. 
T(z) (z,p)>0 


Here A, is the first differential Beltrami parameter. Then 


S us d d 
ümór; ^ de sj J Pons - E piz(p)dsp 
F(z) lre ) 
bz; 
rt J < Ai(piz(p), p) — 2p:z(p)p, z > dwp 
{z,p)>0 
+ iP “sf < Ai(piz(p)), p > pjdsp 
oo 
Ee J < Ai(piz(p), p) — 2piz(p)p, £ > dwp. 
7 (z,p)>0 


From here it follows that f(x) belongs to class C?(R*\0). Consistently applying formula 
(3.11) to integral by I'(z) we obtain for z(p) € C*(w) function f(z) € C**1 (RT X0) and, 
hence f(n) € Cw). 

For the proof of the analyticity we use the following fact (see (Sobolev, 1974)): let on 


sphere w a summable function A(p) be given and let h(p) = bu Y,(p) be its expansion in 


spherical harmonics. Then the function h(p) is analytic on Par w if and only if 


IY«(p| < Cexp(-vk), k&=0,1,2..., 
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and C > 0, v > 0 are constants. 
Suppose that f(n) is an analytic function on w and 


foem. Yale) = "ES S AOP pn sla, 
k=0 


is its expansion in spherical harmonics Y;(n). 
From Theorem 3.2 it follows that it is possible to write the solution of equation (3.10) 


as a series: i i (2k +2)! 
+ 
ET es gy cadi uS 
ARS aeta) Ta 2 D (2k 3i 
1 &, (2+2)! 
tara 24 D a apne”) 
Because of the fact mentioned d and f(n) being analytic |Y;(p)] € Cexp(—vk), 
k —0,1,2..., with certain positive constant C and v. Therefore 


Ya.(n) 


2k 4-2)! "m 
|n ie prem se ens 
(2k + 2)!! —(2k+ Io, 
e ipse] «e 


for a certain constant Cj > 0 and n > 0. 

In this way, terms of expansion of the function z(n) by spherical harmonics decrease 
exponentially. Hence z(n) is analytic on w. If z(n) is an analytic function, f(n) is analytic 
as well. The theorem is proved. 

Consider an equation of the first kind with respect to the odd function v(p), p € w, 
on sphere w: 


ae x nom (3.12) 


w(n) 


Theorem 3.4. If h(n) is an odd function of class C3(w) at w, then there exists a unique 
continuous solution of equation (3.12) and the formula of inversion holds: 


Vi me ct J SAND ads (3.13) 


an <n,p> 
(n,p)»0 P 


PROOF. Suppose at first, that function h(n) is an odd analytic function at w. Expand it 
in a series in spherical harmonics: 


h(n) = » Yos (n). (3.14) 


Since spherical harmonics Y?!4, are eigenfunctions, and the corresponding eigenvalues are 
(- y (BE = 1p" 
(2k + 2)!! 


represented as series 


(see the proof of Theorem 3.2), the solution of equation (3.12) can be 


vp = SD oe panel) (3.15) 
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and, as the terms of this serie decrease exponentially, it is an analytic function. 
Set expansion (3.14) in integral 


1 Ah(p) 
E" n 
(n,p)>0 


use equality 
AY (p) = —(2k + 1)(2k + 2) Yos (p), 


integrate term by term, pass on to spherical coordinates y, 7 with respect to pole n and 
use equalities of (Gradshtein and Ryzhik, 1980): 


Yorti(y.7) = Pai (cos y)Yari (n), 


1 
Pya(t), , 4 QE B 
/ pod (“oa E652 


We obtain 


AY Y 
ER tees A HM 
27 i PT p Qn : «np» p 


LS AY; (p) 
Qn = i «np» ap 
k=O in, p)>0 


a alie Yok+ı (p) 
= -zz LOk Drs " n dun 
nP, 


1 j TY, (y, T) sin 
= —— Y + 1)(2k + 2) far | EHE T 
2m i-o COS y 


5 
oe P. 
- Vk iar 2) [ Dny, (n) 
k=0 0 COS ^y 


(2k)! 


= Sak + 1)(2k + DD i s Yen) 


k=0 


Ea 2k 4-2)! 
- YCO ripreso). 


Compare the obtained series and expansion (3.15), and conclude that the solution of 


equation (3.12) is given by formula (3.13). 
It is possible to represent solution v(n) in the form 


v(n) = h(n) — f tg yh"(y,n)dy, (3.16) 


where À" = Gh, n). To obtain (3. 16) start with (3.13) in spherical coordinates y, 7 
with respect to pele n. Using an expression for the Laplace-Beltrami operator in spherical 


70 Chapter 3 


coordinates, we have 


1 Í Ah(p) | 


2m «np» 


v(n) 


P 
(n,p)>0 


= 2v 
trc Ahip(y,7)) . 
= E / SET sin yd 
$ 1 Qn 
= - / tg 7 H J d dy 
2h Oh a 8k 
- -jedi I eR s ; 


Dj 
- Jes (y,n)d4 — feton- z J E 
0 


= À(0,n) - À( Sn) — Je yh" (y, n)dy 
d 


1 $ tgy [Oh(s2r) | 0h(,0) 
J wa | ðr 9s ^h 


Since 


2r) OMD) e io gg 


Or 


and because of the oddness of function h(n) average h (3,0) = 0, finally we obtain 


v(n) = h(n) — f tg sh (9, ndn. 


If function h(n) belongs to class C?(w), then approximating it in norm of space C?(w) by 
a sequence of analytic functions, we establish the correctness of the theorem. 


Theorem 3.5. If in equation (3.10) a Z 0, b # 0, a + b £ 0 and function f(n) belongs 
to class C*(w), k > 4, then the formula of inversion holds 


"M. (A*2)fi(p)x mp >| 
z(n) = "825i J ANT dwp 
(n.p)? >0 t=0 
ud p A as. (3.17) 


4nb 5 <n,p> 


PROOF. As in the proof for the even and odd parts of the function z(p) in Theorem 3.1, 
we obtain equations 


fip) = f 1< n.p >12*(P) dup, (3.18) 
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f'(n)-b Te 1(p)dwp. (3.19) 
w(n) 
It is known that equation (3.18) is uniquely solvable and its solution is given by (Blaschke, 
1916; Pogorelov, 1976) 


+n) = 18 J (A *2)h(px mpl, 
ab s e V<n,p >? -t p 


8r? ôt 
(np)? >t t=0 


2)Al (yn 2 dy 


= (A+ 2)flp Mu J fla (3.20) 


cos y 


In this way, the even part of the function is determined uniquely and constructively. 
The odd part of z^ (n) is determined uniquely from equation (3.19) by 


z (n)— 


e (n) i tg yf "dy. (3.21) 


EJ ee ES xi 


Summing (3.20) and (3.21) by equality z(n) = z+ + 2^ we obtain (3.17). 


REMARK I. f a = 0 or a + b = 0, then from equation (3.10) only the odd part z^ (n) of 
the unknown function z(n) is determined uniquely, and if b = 0, then only the even part 
z* (n) can be found uniquely. 


Studying the inversion formula of equation (3.10) leads to the following result. 


Theorem 3.6. Ifa +0, b X 0, a+b + 0 and function f(n) € C*(w), k > 4, then solution 
z(n) of equation (3.10) belongs to at least class C*-*(w) and the following estimate holds: 


lze S M IHf(le« (3.22) 


where constant M depends only on a and b. 


PROOF. By the previous theorem it is possible to write the solution of equation (3.10) in 
the form 


7 (A 2A nm) 
on) = (AF 2) Al—) + z— [a 
tah ln i a | yf "da. 


Since function fi(n) is even on w, equality (A + 2)fi (7, n)]| = 0 is fulfilled. Us- 
=F 


ing this equality and the existence of a continuous derivative (A + 2) fily, n) " we can 
represent it as follows: 


[a 2G 2)| = (1-5) f[iasni( (1-5)+5 DI dt. 
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In the same way, because of the oddness of function f^ (n) to w, f", QM = 0 and 
-$ 


arok G oee e 


From the last formula it follows that z(n) belongs to class C^-*(u). 
Prove the estimate. Using again formula (3.23), we have 


x 


1 1 PE -y 
lele S g A+ DA MNc + 1- / 21793. 


COS 4 
d E (A +2)fi) lza Qn |b] |o Mhe 
pe 
1 : T 
tsa | G- Ji dy sup |f^"|. 
hew 
Estimate every addend in the right-hand part of this inequality: 
1 2 2b — 
A+2 € — — 
ICA  2)f()leu € Jal IA flou + la| M loot Tn jal a+] HCOE 
~ it 
sap (A +2)fi] |< qq la bi jsup [/" |, 
sup B < sup |Z" : 
yn yn 
(te )t dcs dnd = - 2G 
2nb ) \2 TS EN C9 EO J cosy Xe 


G is the Catalan constant. From these reduced inequalities the estimate follows: 


ı (13,1, Ml 
an (a A lèl T la| la + 3) IF) lew + ——_ 
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The theorem is proved. 

On the basis of the inversion formula of equation (3.10) by a standard manner more 
general integral equations of the first kind may be examined, including nonlinear ones. 
Consider, for instance, 


fin) = [ P(n, p): (phos + Tz, 


where P(n,p) = a|« n,p »| - bx(« n,p >), and T be a certain operator, using the 
inversion formula, we have 
z(n) = f(n) + Bz, (3.24) 

where f = P-! f, B = P-'T. Relation (3.24) is an equation of the second kind and, for 
instance, by corresponding restriction to T operator B may be contracting, which allows 
us to use methods of the theory of fixed points. 

Another application of the inversion formula can be connected to the problem of 
integral geometry. Consider, for example, the following problem: to find function u(n) € 
C™(w), n € w, if function y(n) is known: 


y(n) = 2a f u(p)ds + b J (A + 2)u(p)dup. (3.25) 
Pn) w(n) 


In the case a # 0, b = 0, we obtain the problem about the recovery of u(p) by its integral 
along large a circle sphere. Minkowski showed through the use of expansion in spherical 
harmonics, that the even part of u(p) is uniquely determined for b = 0. Later Funk 
obtained a solution of this problem in a closed form, using the inversion formula of the 
Abel integral equation (see (Funk, 1916)). 


Theorem 3.7. Let a £0, b 4 0,a4- 6 Z 0, v(n) € C”(w) and f ne(n)do = 0. Then 


there exists a solution u(p) of equation (3.25), u(p) € C? (w) and inversion formula holds 
1 
u(p) = q | *) {(1- < np »)In(1- <n,p>)} dwp, 


where 


i) em J (A + 2)¢i(p) |< n.» >| 


= dw 
V< n, p >? —1 5 


t=0 


~ 8x? dt 
(n,p)?>¢ 


d f Av (p) 
4nb d «np» 


dwp. 


PROOF. . With the use of Green's formula we obtain 


2 J u(p)ds = T |< n, p »1(A + 2)u(p)duy. (3.26) 
Pin) w(n) 


Therefore equation (3.25) with respect to function u(p) in accordance to (3.25) and (3.26) 
acquires the form 


e(n) = f {al<n,p >| + bx(< n,p >)} (A+ 2)u(p)dwp. (3.27) 
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The solution of equation (3.27) exists and is unique, if 


[ ^3. =0. (3.28) 


w 


Indeed, let 


= Y Y,(n), Y¥,(n) = — fe (p)Pk(< n, p >)dwp 


w 


be a series of expansion in spherical harmonics of the left-hand part of equation (3.27). 
Because of (3.28) for the first harmonic we find Y,(n) = 0. Then function 


oo Az} 
P= popra Oh 


where eigenvalues \, of equation (3.10) are solutions of equation (3.27). If y(n) is infinitely 
differentiable on w, then function A' is bounded to w, for all values of l of operator A. 
That is why, and from the estimate of the derivative of spherical harmonics (Sobolev, 
1974) 


|A°Y,(p)| € Cr le] lo| 2 m 


Lal)? 
C > 0 is constant, infinite differentiability of solution u(p) follows. 

Consider equation (3.27) as an integral equation of the first kind with respect to 
(A + 2)u(p). According to the inversion formula of equation (3.10) we find 


m = 1d (A + 2)¢r(p) |< n, p >| 
(A*2wm)-90)- - gay [ UT duy 
(n,p)? >t t=0 


_ J Ae" (p) 4 
r ips x 


Applying the Weingarten formula (Blaschke, 1930; Pogorelov, 1973), from (3.29), we 
obtain the wanted result: 


u(p) = = f 9(9)(0— < n,p >)In(1— < n, p >)}dwsp: 


3.3 UNIQUENESS 


Here some results are given on the univalent definiteness of a closed convex surface, 
if the functionals of its orthogonal projections and illuminated parts are known. The 
investigation method consists of the following: at first the problem has to be reduced to 
the integral equation (3.1), and then the results are used to solvs the problem of Minkowski 
and Christoffel. All surfaces are, as a rule, supposed to be strictly convex and sufficiently 
smooth. 
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Consider the function 
p(n) = 2aF(n) + bS(n), 
where F(n) is the area of orthogonal projection of the closed convex surface B on the 
plane P(n) = {xz € R?,« z,n >= 0}, S(n) is the area of the illuminated part of the 
surface in the direction n, a and 6 are constants. 
F(n) and S(n), can be represented as 


F(n) = > n n,p >|G(B, dwp), S(n)- J G(B, dwp), 


a w(n) 


where G(B, -) is a surface function. Using these representations it is possible write function 
(n) as 
eins af |< n,p >|G(B, dwp) + b J G(B, dwp). (3.30) 
w w(n) 
A closed convex surface is defined uniquely, apart from a translation, by setting the surface 


function (Aleksandrov, 1937a,b; 1938a,b). Therefore and because of the uniqueness of the 
solution of equation (3.30) (see Theorem 3.1) the following theorem takes place. 


Theorem 3.8. A closed convex surface is defined uniquely except for a translation by 
function e(n), n € v, a Z0, bz 0, a - b z 0. 


Let G,(B,-) be the first function of curvature of a closed convex surface B, let L(n) be 
the length of a boundary orthogonal projection of the surface B on the plane P(n), M(n) 
is the average integral curvature of the illuminated part of the surface in the direction n. 
Functions L(n) and M (n) are determined by 


1 1 
Lo - | I< mp»|O(B,duy), M(n) = 5 f 609,429. 
w w(n) 


Consider the function 


p(n) = 2aL(n) + 25M (n), 


where a and b are constants, a # 0, b # 0, a+b # 0. Through the use of the representation 
for L(n) and M(n) (3.30) can be written in the form 


p(n) = af |< n, p >|Gi(B, dwp) + b " G(B, dwp). (3.31) 
4 w(n) 


Expression (3.31) is an equation of the first kind with respect to the function of curvature 
G,(B,-). The solution of the given equation is unique using Theorem 3.1, but since a 
closed convex surface is defined uniquely, except for a translation, by setting first the 
function curvature, the following theorem is equitable. 


Theorem 3.9. A closed convex surface B is defined uniquely apart from a translation 
by function p(n), n € w. 


REMARK I. f surface B is centrally-symmetric, then Theorems 3.8 and 3.9 are equitable 
at b — 0, that is, a convex centrally-symmetric surface B is uniquely defined by the area 
and length of its orthogonal projection on a plane in all directions. These are well-known 
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facts (see (Blaschke, 1916)). A noncentrally symmetric convex surface is not uniquely 
defined by the area of its orthogonal projection. So, a piecewise analytic convex surface 


of revolution, which form has parametric representation 


A T 
z = ——cosq, =0, ae Ji — 0, 
Và P, y n p y 


dia d "sese 
is, when being a sphere, a surface of constant luminance, that is, it has a constant area of 
projection on any plane (see (Blaschke, 1916)). A convex noncentrally-symmetric surface 
is also not defined uniquely by the length boundaries of its orthogonal projection. For 
example, a convex surface, having support function H(z) = h(z) + |x|, where A(z) is an 
odd function, is a surface of constant coverage, i.e., it has a constant length of projection 
on any plane. 


Theorems 3.8 and 3.9 establish the uniqueness of reconstruction of the convex surface 
B, but do not give a constructive method to determine its form. When surface B is 
sufficiently smooth, the problem of determining its form by y(n) or (n) is reduced to 
the inversion of equations (3.30), (3.31) and to solve the equation in partial derivatives of 
the second order for the support function of the surface. 

Let Ri(n), Ro(n) be the principal radii of curvature of a closed strictly convex, and at 
least twice continuously differentiable surface B as function of unit vector n normal to B. 
In such cases the surface function G(B,-) has a continuous density z(n) = R,(n)R2(n) 
and ¢(n) is written in as 


v(n) = aza] |< n, p >| Ri(p) Ra(p)dwp + b J R, (p) Ra(p)duy. (3.32) 


w(n) 


Consider (3.32) as an equation of the first kind with respect to the product of the principal 
radii of curvature R,(p)R2(p). Using Theorem 3.5 leads to following result: 


Theorem 3.10. Let a 4 0,6 Z 0, a+b z 0, let B be a closed convex regular surface. 
Then: 


l d A42 «np» 
Ry(n)Ra(n) = — Bni di Cep ala, 
<n p>? >t Y »P nf 
Ae (p) a 


i LE m 


Hy(n) Hix(n) 
Hafn) Hz(n) 


H2(n) Has(n) 


Ha(n) Haa (3.33) 


Ry(n)Ra(n) = | 


pas H,3(n) 
Han) Hy(n) 


2 
Here H;;(n) are the derivatives Fade 


of the support function H(z), calculated at z = 


new. 
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From the theorem it follows, that the recovering of a regular surface B by the function 
y(n) is reduced to the solution of equation (3.33) for its support function. Formula (3.33) 
is well-known (Bakel'man et al. , 1973). 

At the conditions to surface B as above, the first function of curvature G1(B,-) has a 
density, which is equal to the sum of the principal radii of curvature. That is why it is 
possible to write function v (n) in the form 


Wa) = a f |< n,p Pl) + Pop), + 8 J o fais. 


w(n) 


Theorem 3.11. Let a £ 0, b 4x 0, a+b #0, and B be a closed strictly convex regular 
surface. Then: 


id) eem J (A + 2il) | mp >| 4 


1 
8r? dt V<n,p >? -t 


<n,p>?>t "E 
J Av (p) 4 
rJ «np zu x 


H(n) = Z f [Ri(p) a9) {(1- < mp »)hn(17 < n,» >)}. 


(3.34) 


Using (3.34) the question recovering regular surface B by the function v(n) is solved 
constructively in explicit form (Materon, 1975). 

Consider now the problem of determining a closed convex surface and a function, 
given at its surface, on the condition that the functionals from its orthogonal projection 
and illuminated portion are given. Similar problems appear also in inverse problems of 
scattering. 

Let B be a closed strictly convex smooth surface in R? and let u(n) be a continuous 
function, given at B as a function of the unit normal vector of this surface. Consider 
functions 


ep) = a f I n,p Iu) Ri) Ro(p)dos b f u(p)Ri(p)Ra(p)dwy, — (3.85) 


w(n) 


exp) = az f |< n.p >|u(p) [Ra(P) + Ra(p)] dos b. f wp) (Fa (p) + Fa(p) dwp, (3.36) 
w(n) 


where Rı, Rz are the principal radii of curvature, a;, b; are constants, a; Æ 0, b; # 0, 
ai + b; # 0, i = 1,2. The geometrical sense of the introduced function at u(p) = 1 are 
described above. 


Theorem 3.12. A closed strictly convex, twice continuously differentiable surface B 
and a continuous function u(n), u(n) > 0, n € w, are uniquely defined by functions vj (n), 
v2(n), n € w (surface B can be restored up to a translation). 
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PROOF. Because of the uniqueness of the solution of equations (3.35) and (3.36) functions 
$,(n) = u(n)Ri R2, 9;(n) = u(n) [Ri + Ry] are uniquely found by functions y(n) and 
Y2(n), respectively. Considing that u(n) > 0 at w we obtain equality 


R,(n)Ro(n) 


Ri(n) + Rin) - an); 


where ®(n) = ©,/® is a given function. It is known that this function defines surface B 
uniquely up to a translation (Pogorelov, 1973). That is why the principal radii of curvature 
R; and Rz are uniquely established. With equality 6; = u F4 Rz we find function u(n). 
The theorem is proved. 

In scattering or emission problem one has to determine both the scatterer B(t), the 
form of which can vary in the course of time t, as well as the function, characterizing the 
density of the source on it. Taking this into account consider equation 


(n, t) = J P(n, p)u(p)ds, (3.37) 
B(t) 


where P(n,p) = a|« n,p >| + óx(«n,p»),a Z 0, b Z 0, a -- b 0, ds, is an element of 
area B(t). 


Theorem 3.13. Let surface B(t) be closed at any t > 0, strictly convex, twice continu- 
ously differentiable, and let the following equality hold: 


d 
ac (p t))| = e [Fi(p,0) + Ra(p,0)], 
t=0 
where Ri(p, t), Ra(p,t) are the principal radii of curvature of surface B(t), t > 0, a £ 0 is 
a constant. Then the continuous function u(p) > 0 and the form of surface B(t) for any 
t > 0 are uniquely defined by function 9(n,t), n € w, t > 0. 


PROOF. Consider relation ds, = Ry(p,t)R2(p,t)dwp equitable for smooth surface, and 
rewrite equation (3.37) as 


O(n,t) = | P(n, pulp) Ra (p, t)Ra(p, deo. (3.38) 


Differentiating equation (3.38) with respect to t and using the conditions of the theorem, 
we obtain 


dé 
dt 


= a f P(n,p)u(p) (Fi (p, 0) + Ra(p, 0)] dwp- (3.39) 
t=0 w 
Function Wi(n,t) = u(n)Ri(n,t)Ro(n,t) from equation (3.38) is uniquely defined by 
function (n, t), but from equation (3.39) function W2(n) = u(n) [Ei (n, 0) + R2(n, 0)] is 
uniquely defined by function ®{(n,0). From this surface B(0) is uniquely defined exept 
for a translation by the function 


Ry (n, 0)R2(n, 0) 


Wine R(n, 0) + Ra(n, 0)’ 


Then the principal radii of curvature Ri(n,0) and R;(n,0) of surface B(0) are uniquely 
found. That is why from W2(n) = u(n) [Ri(n,0) + Re(n,0)] we find u(n). Knowing u(n) 
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from W,(n,t) = u(n)R,(n,t)R2(n,t) we define Ri(n,t)R2(n,t). Using the uniqueness of 
the solution of the Minkowski problem by function W, /u we uniquely define surface B(t) 
for any t > 0, apart from a translation. The theorem is proved. 

Let the area and length of a boundary of an orthogonal projection of convex surface 
B on plain in all directions be given. The possibility to define surface B by the functions 
F(n), n € w, L(n), n € w is given by the following theorem. 


Theorem 3.14. Let B, and B, be closed convex analytic surfaces in R?. If Fi(n) = 
F,(n), Li(n) = L2(n) for any n € w then surfaces Bı and B; are equal. 


ProoF. Area F(n) and length L(n) of the boundary of a closed convex, at least twice 
continuously differentiable surface B are defined by 


F(n) = È f [Rilp)Ro(p)]* le mp >| dup, (3.40) 


Ln) = = f IRP) + Ra(p)]* I< sp >l dup, (3.41) 
where Hj and HR; are the principal radii of curvature as functions of the unit normal 
vector. Let H(x) be a support function of the surface B. It is well known that E, and 
R; are eigenvalues of matrix |H;;|| of the second partial derivatives of the function H(z). 
That is why we can write 

AH = —(R + R3), 


from which we obtain that 
AH* = —(R, + R3)*. (3.42) 


here A is the Laplace operator. Hence, if function L(n) is given, then from equation (3.41) 
(Rı + R2)* is uniquely defined, and from the equation above the even part of support 
function H(z) is uniquely defined (Christoffel theorem). 

In order to prove the theorem, it is sufficient to establish that from setting F(n) (or 
from setting [Ry(n)R2(n)]* ) the unique definiteness of H- follows up to an arbitrary 
linear function (addition of a linear function to H means a new choice of the origin of the 
coordinate system). 

Since Rı Rz is equal to the sum of main minors of matrix |H.;||, then filling Hi; as a 
sum of even and odd parts and picking out parameters u and v on sphere w by specific 
manner (Buseman, 1958), for odd part, H~, of the support function we obtain equation 


HH. 7 Hz,’ = (uv), (3.43) 


where #(u, v) is known function. 

Our problem is the definition of H^ up to a linear function and multiplication by (-1) 
(this is equivalent to rolling of surface B) from equation (3.43). 

Show, that there exists a point (uo, vo) such that (uo, vo) > 0 (it is assumed that 
surface B does not have a center of symmetry and therefore H^ is non-linear). 

Since H- is an odd function, there exists a closed curve | = (u(t), v(t)) such that 
in domain D, the bounded curve l is simply connected and H^ = 0 on l. Therefore, 
if everywhere on w function p(u,v) € 0, then in a domain D we would obtain surface 
H- = H^ (u,v) of nonpositive curvature, equal zero at l, which is impossible. 
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Let (uo, vo) be a point, in which % > 0. Let K denote a circle of certain radius centered 
in point (uo, vo), in which  » 0. Because of the continuity of (u,v) such a circle exists. 

If now H, and H; are support functions of surfaces Bı and Bz and convex in the circle 
in one side (it is always possible to roll one of the surfaces), then adding, for instance, a 
certain linear function to Hj, one can obtain H, = H» on a certain closed curve lj, lying 
inside circle K. Since H} = H} everywhere on w, then Hy = Hy on h. 

That is why, taking into account that H7 and Hy satisfy the same equation (3.43) 
and furthermore Hy = Hj at lj, then Hj = Hj everywhere in the domain, bounded by 
lı (Bakel’man, 1965). And by supposing B; and B; analytic everywhere on the sphere 
Hy = Hj. The theorem is proved. 


REMARK . (Campi, 1986) giives un example showing that the condition of analyticity in 
this theorem is essential. 

In 1926 (Matzumara, 1926) proved, that a closed convex surface, of which area and 
length of boundaries of its orthogonal projections on any plane are constant, is a sphere. 


In connection with this result it is naturally to expect that if isoperimetric defect 
L? — 4x F(n) is small for each isoperimetric projection B(n) of surface B, then surface 
B differs a little from a sphere. Let R and r be radii of described and inscribed balls of 
convex surface. 


Theorem 3.15. If L? - 4sF(n) X &?,n € o, then R-r < E, 
PROOF. We again use representations 
1 
F(n)= -3 [I< np>lG(B dup), L(n)= f H(p)dsp 
u P(n) 


for the area and the length of the boundary of an orthogonal projection. 
Let 


1 1 
s=- / F(n)e, M =z- / L(n)du,.. 


Numbers 5 and M are the equal area and average integral curvature of surface B, respec- 
tively, and connected by isoperimetric inequality (Blaschke, 1916) 


M? — 458 > Av? (R — ry*. 
Using this inequality and Cauchy-Bunyakowsky inequality, we have 


(R- ry 


lA 


2 
ii |M? - 4x5] = E zi ( Í Hà, -4 J F(n)dw, 


1 1 1 
| cs PESE 2 c 
E Fak (das f dum 1 [erras 
E 
4r? 


2 
J [L?(n) — 4x F(n)| du, < aser "iB 


< , 
T2 


wW 


i € : 2 Ss 
that is ER — r < —, as stated. Notice that no regularity condition of surface B has been 


T 
used in this proof. The theorem is proved. 
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3.4 EXISTENCE 

This section is devoted to the theorems of existence of closed convex surfaces with a given 
functional of its orthogonal projection and illuminated parts. 

Theorem 3.16. Let F(n) be a continuous, strictly positive, and even function on unit 
sphere w such that — |z| F (à ) is conditionally positive defined on R?. Then there 


exists a unique closed convex centrally-symmetric surface B (except for a translation), for 
which function F(n) is the area of its orthogonal projection on plane < z,n >= 0, n € w. 
If function F(n) is k times continuously differentiable (k > 4, analytic and 


a f asyriplemprla| co (344) 
dt J<n,p >? -t 
<n,p>2>t t=0 


then surface B is at least k — 4 times continuously differentiable. 


PROOF. Let us prove the existence. By theorem (Matheron, 1975) there exists a unique 
symmetrical measure yz at w such that function F(n) permits representation in the form 


F(n) => fI» »lu (dos). (345) 


Since measure y is symmetrical on w, and function F(n) is continuous and strictly positive, 
then the following conditions are fulfilled: 


f ru (den) = 0, [I< np ola (dwp) > 0, nEw. 


w 


Hence, as it is proved in (Aleksandrov, 1937a,b; 1938a,b), there exists a unique closed 
convex surface B, the surface function of which is given by measure p. Surface B is 
centrally-symmetric, because measure p is symmetric. But then, by representation (3.45), 
F(n) is the area of the orthogonal projection of this surface B on planes < z,n >= 0. 

Pass on to the proof of the regularity of surface B. Let F(n) € C*(w). Consider 
function 


NE X (A 2)F(p)|« n, p >| 
z(n) = "aii J Enee = dwp) <0. 
<n,p>?>t t=0 


This function is the unique solution of 


= |< n, p >| 2(p)dwp 
2 w 

and, by Theorem 3.6, it belongs to class C*-*(w) (and is analytic if F(n) is analytic). 
Then with necessity function z(n) is the density of measure p, but since measure p is a 
surface function, z(n) = K ^! (n), where K (n) is the Gaussian curvature of surface B. This 
Gaussian curvature K (n) of surface B belongs to class C*~4(w) and, by condition (3.44), 
is strictly positive on w. Therefore from the result of the regular solution of the Minkowski 
problem (Pogorelov, 1973; 1978) it follows: if k — 4, then K (n) is continuous, and the 
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surface is smooth and of class C; if k = 5, then K(n) € C'(w) and this inplicates that 
the surfacebelongs, at least, to class C?, if k = 6, then K(n) is (k — 4) times continuously 
differentiable, hence surface B is at least (k —3) times continuously differentiable. Finally, 
if F(n) is analytic, then K (n) is a positive analytic function and surface B is analytic. 
The theorem is proved. 


Theorem 3.17. Assume a regular (k times continuously differentiable, k < 7, or ana- 
lytic) strictly positive even function L(n) on sphere w, satisfying condition 


J(n) «0, J(n)-J;(n) <0, new, 


where 

d (A + 2)L(p) |< n.p >I 

— SS dup] , 
V<n,p >? -t 


€«np»?»t t=0 
and index s denotes differentiation in point n along the arc of a large circle. Then 
there exists unique (up to translation) regular (belonging to class C*-*(w), 0 < a < 1, 
or accordingly analytic) closed convex centrally-symmetric surface B, which length of 
boundary of the orthogonal projection on plane < z,n >= 0 is equal to L(n), n € w. 


PRoor. By its evenness and regularity function L(n) admits representation in the form 


Ln) = 5 f Ie nip »lp(p)ào, (3.46) 


w 


with a certain even, not necessarily positive, function p(p). Equation (3.46) is uniquely 
solvable and function p(p) is determined by 


u^ ilid (4+ 2)L(p)|« n, p >| 
(n) = - 4r? dt i. V<n,p >? —t dum 


<n,p>? >t t=0 


By condition of Theorem 3.17 we have 


p(n) 20, p(n)— pj, (n) 2 0, 


and furthermore 
J np(n)du, = 0. 
since p(n) is an even function on w. 

Thus, function p(n) satisfies all condition of the existence theorem of a closed convex 
surface with given sum p(n) of the principal radii of curvature (Pogorelov, 1973, 1978). 
According to this theorem there exists a unique (except for a translation) closed centrally- 
symmetric surface B with a given sum of the principal radii of curvature. 

If L(n) € C*(w), k > 7, then function p(n) € C^-*(w), but then the surface of p(n) 
belongs to class C*-**(w), 0 < o < 1 (see (Pogorelov, 1973)). If L(n) is analytic, then 
p(n) is analytic and, hence, B is analytic (Pogorelov, 1973, 1978). The theorem is proved. 


In this theorem tougher condition are set to function L(n) than to function F(n) in 
Theorem 3.16. This circumstance is the result of condition in the problem of the existence 
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of a convex surface with a given sum of the principal radii of curvature, which are more 
limited with respect to the Minkowski problem. 

Proof of next two theorems are conducted in a similar manner with use of the inversion 
formula of the integral equation (3.10), the existence theorem of a convex surface with the 
main function of curvature radius given and results of the regular solution of problems of 
Minkowski and Christoffe 


Theorem 3.19. Let a regular (k times continuously differentiable, k > 4, or analytic) 
function p(n) on unit sphere w in R? be given, satisfying conditions 


1. 
J 00e. =0, 
2. 
id (A + 2)&(p) |< n.p >| l Ag" (p) 
———— B—M————————d —— J — ——— duy > 0, 
87? dt Í V< n, p >? -t B 47? <n,p> “p 
<n,p>? >t t=0 w(n) 
necu, 
where 


exp) = “o*(p) = mec J p(p)dwp, 


a and b are constant, a # 0, b # 0, a+b #0. Then there exists a regular, at least 
(k — 3) times continuously differentiable (or respectively analytic) closed convex surface 
B, for which y(n) is equal to the sum of the area of its orthogonal projection on plane 
< z,n >= 0 multiplied by 2a, and the area of its illuminated part in the direction n 
multiplied by b. Surface B is determined uniquely, a part from a translation. 


Theorem 3.20. Let regular (k times continuously differentiable, k > 7, or analytic) 
function (n) on unit sphere w in R? be given, satisfying conditions 


1. 
[$6092 =0, 
2. 
p(n) $0, pln)- p(n) <0, nEw, 
where 
l d (A + 2) (p) |< n.p >| 1 Av" (p) 
Er LA T APPI P dies wee J IV P) diss 
p(n) 8r? dt / V<n,p >? -t ac 47? «np» P 
<n,p>? >t t=0 w(n) 


tb) = 00) - rore | 99d, 


a and b are constant, a # 0, b X 0, a+b 4 0. Index s denotes differentiation in point 
n on the arc of any large circle. Then there exists a unique (up to translation ) regular 
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(belonging to class C*??, 0 < a < 1 or respectively analytic) closed convex surface B, 
for which function y(n) is equal to the sum of the length of the boundary of its orthogonal 
projection on plane < z,n >= 0, multiplied by 2a, and the average integral curvature of 
its illuminated part in the direction n, multiplied by 2b. Surface B is determined uniquely, 
except for a translation. 


3.5 STABILITY 


Theorems of stability definition of a convex surface by functionals of projections and 
illuminated parts are of interest from two points of view. From the standpoint of ge- 
ometry there are very importanestimates of the proximity of surfaces, which are close 
to the functionals determining them. On the other hand some of the functionals appear 
in investigations of inverse problem of scattering, used in some approximation. In this 
connection theorems of stability allow us to confirm the stability definition of a scatterer 
by data of scattering. 

At first we provide an estimate of the stability definition ofa regular convex surface 


by function y(n) = 2aF(n) + bS(n). 
Theorem 3.21. If two closed, strictly convex regular surfaces Bı and B, have functions 
yi(n) = 2aFi(n) + bSi(n), yo(n) = 2aFo(n) + bS;(n), a Z 0, b #0, a+b #0 k times 
continuously differentiable, k > 4, and the inequality 

lex(n) — vi(n n)lc«u € Se, e>0 


is fulfilled, then at a certain location on surfaces B and B; the difference of their support 
functions H,(n), H(n) satisfies the inequality 


\|H2(n) ia Hi(n)leu, < Ces, 
where C > 0 is a certain constant. 


PROOF. Let z;(n) denote the product of the principal radii of curvature of surface Bj, 
j — 1,2. Then functions y;(n) can be written in the form 


ej(n) =a J |< n,p >| zi(p)dwy +b f z;(p)dwp. (3.47) 
w(n) 


Using inversion formula (3.17), we have 


z 
2 


ap) = gl +2980) [EE Cae 


COS y 


1 ~-i 
tz (n)— A fur (y,n)dy, 


where 


$;(p) = “9*(p) - E5 / p(p)dwp. 
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Estimate the difference z; — z1. By Theorem 3.6 


lz2(n) — a(n)llou, S Cla, b) lox(n) — exeo » (3.48) 


where C (a, b) is constant, depending on a and b. 
Let Q be any set of non-zero measure at w. Since surfaces By and B, are smooth, 
then it is possible to write their surface functions GU)(Q) in the form 


GQ) = | 25(n)dun. 
Therefore 


[G(Q) - GQ)| < ft) - 2 (0) dwn 


< mes Q max |za(n) — ai (n)| < 4r llza(n) 7 aln)llow)- 


Hence, taking into account estimates (3.48) and |ly2(n) — e1(n)le«u, € € we obtain 


IG? (Q) m G(Q)| < C46, where C, = 47C (a,b). Using now the Volkov theorem on the 
stability of the solution of the Minkowski problem (Volkov, 1963), we conclude, that at a 
certain location on surfaces B, and B; the difference of their support functions H, and 


Hy satisfies the inequality 
|H2(n) — Hi(n)| € Ce”, 


with a certain constant C > 0. The theorem is proved. 

To obtain the estimate of the stability definition of a regular convex surface by function 
p(n) = 2aL(n) + 25M (n) an estimate of the stability of the solution of the Christoffel 
problem is used: if sums p(n) and p2(n) of the principal radii of curvature of closed 
convex surfaces B, and B, differ by not more than € then at a certain location their 
support functions H, and H, differ by not more than (In2 — 1) €. From this and the 


4 
theorem of stability of the solution of the integral equation it follows 


O(n) =a f |< n.p >I [Ra(p)  Ra(p)] dup +b f [Ri(p) + Ra(p)] dwp 


u w(n) 
ata #0,b40,a0+6b0. 


Theorem 3.22. If two closed, strictly convex regular surfaces B, and B; have functions 
V;(n) = 2aL;(n) + 2b5M;(n) k times continuously differentiable, k > 4, and inequality 


lis(n) - (n) Se €>0 


holds, then at a certain location on surfaces B, and B, the difference of their support 
functions Hi, H; satisfies inequality 


1 
[Fa(n) — H Gan, <C (12-7) & 


where C > 0 is a certain constant. 


CHAPTER 4 


Integral Geometry 


In this chapter the problems of integral geometry are considered. The methods and results 
concerning this type problems are presented. Existence and uniqueness theorems, as well 
as some inversion formulas are given. Some results are proved. 

The problems of determing unknown objects (functions, differential forms, tensor 
fields) their integral known on manifolds (curves, surfaces) are called the problems of 
the integral geometry. 

The problems of the integral geometry are connected with differential equations, in- 
verse problems, group representations and one can find them in the geophysics, astronomy, 
or medicine. 

The first results here were received by (Funk, 1916; Radon, 1917). They founded the 
integral geometry and discovered the inversion formulas. 

Later, new problems of integral geometry were set, many theorems of the uniqueness 
and existence were proved, new inversion formulas were discovered and new applications 
were given. 

The problems of integral geometry are studied in the works of (Blaschke, 1916; John, 
1955; Gelfand et al., 1966; 1969; 1980; Gelfand and Goncharov, 1987; Kostelyanec and 
Reshethyak, 1954; Helgason, 1980; Semyanstyi, 1960; 1961; 1966; Kirillov, 1961; Solman, 
1976; Lavrent'ev and Romanov, 1966; Lavrent'ev et al., 1970; Lavrent'ev and Bukhgeim, 
1973; Lavrent'ev et al., 1986; Romanov, 1967; 1978; 1987; Lavrent'ev and Anikonov, 
1967; Anikonov, 1969a; 1969b; 1978a; 1978b; 1982; 1983; Anikonov and Romanov, 1979; 
Anikonov and Pestov, 1990; Anikonov and Stepanov, 1991; Anikonov and Shneiberg, 
1991; Bukhgeim, 1972; 1979; 1983; Mukhometov, 1977; Blagoveshchenskii, 1986; Pestov, 
1985; Pestov and Sharafutdinov, 1987; Sharafutdinov, 1989; Goncharov, 1988; Palamodov 
and Denisjuk, 1988). 

Applications and information are given in the papers by (Tikhonov et al. , 1987; Nat- 
terer, 1986). One general problem of the integral geometry is to find the function 


O^ f 


oo n m > 
LEVER): pide aar e| c OX ded 


if we know the function 
f= f Hedu, ver’, 
B(y) 


where B(y) C R” is a smooth manifold, y € R”, and y is a measure on B(y) (Gelfand 
et al., 1966). The Radon transform is 
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flep)= f fe» y= (op), 


(w,z)=p 
Biy)={z:(w,z)=p}, wE R”, l|-21 pcR, zcR* 


In this chapter we shall discuss some methods of the integral geometry, especially in the 
curve case and some applications. At first we give some review. 


4.1 INVERSION FORMULAS 


4.1.1. The Radon formulas 


The inversion formulas in problems of the integral geometry play a very important role 
in the theory and applications. 


Theorem 4.1. (Radon, 1917; Gelfand et al., 1966; Helgason, 1980). Let a function 
f(w,p) be the Radon transform 


flwp)= | fm. 


(z,w)=p 
Then the function f(w,p), w € 2, Q = {w, |w| = 1}, p € R, satisfies the conditions: 
1. f(w,p) € C*(Q x R), 


2. . 
ai Ym > 0 


> 
Op = , a 2 0, 


sup(1 + |p|”) 
peR 


^ 


3. J f(w,p)p*dp is a homogeneous polynomial of the degree k depending on w = 
R 


(w1,..., Wn) and the inversion Radon formula 


—4T (n-1)/2 n "re Ki 
ta) = AO TD a [ j(w,(w,2))aw 
T(5 ) Q 


holds, where A"7* is the Laplace operator of degree (n — 1)/2. 


Other similar formulas are contained in the papers by (Semyanistyi, 1960; 1961; 1966; 
Gelfand et al. , 1980). 

If n = 2, z1 = T, t2 =y, £? +y? <1 
Ken- [fed 0<p<2, -1<p<1, 


z cos y+y sin y=p 
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then the different inversion formulas may be written on the basis of Theorem 4.1. 


f(z,y) 


Gr a ] e | festo nts 


T 1 ^ 
1 Of(y,p) 1 
ies Le ced 
Gp | | op po 


oo oo 1 
1 à : 
= J Ji (arctan = — te elt Voit? qt | e tiezs qo dw, 
(4r)? LA EA We 


rn = zcosq + ysin g. 


4.1.2. The other formulas 


In this section we consider the problem of the integral geometry to find the function f(z), 
z € R? such that 


fn [ rt, 


where y is a line passing through a smooth curve F C R?. The function f(y) is supposed 
to be known for all lines such that 4 NT # 0. 


Theorem 4.2. (Kirillov, 1961) A function 


Q^ f 


oo 3 m 
f(z) € C™(R*), Sipi JE ðar 


«oo, Ym>0, a20 


is restored by f(y) and the inversion formulas hold if and only if (z : (z,w) = p) NT #0 
for almost all (w, p), w| ^ 1, p ER. 


(Tuy, 1983) discovered the inversion formulas for the functions with compact support 
when the curve 1s bounded, provided that almost every hyperplane intersecting the sup- 
port of the unknown function meets the curve I' transversally in some point. (Finch, 
1985) investigated the analogous inversion procedure for the class of the curves which did 
not satisfy the hypotheses of Tuy. 

(Blagoveshchenskii, 1986) found the inversion formulas when the curve I is a circle 
in R?. He found essentially new formulas. The uniqueness theorems for the integral 


geometry problems when unknown functions have a compact support were studied by 
(Anikonov, 1978; Helgason, 1980; Leahy et al. , 1979). 
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4.1.3. Integral geometry on the sphere 


The problems of integral geometry on a sphere were considered in the papers of (Blaschke, 
1916; Semyanistyi, 1960; 1961; 1966; Helgason, 1980). The unknown functions were 
supposed to be even on the sphere. 

Let n, p be a point of the unit sphere S in R? with the center at the origin, (n, p) the 
inner product n € S, p € S, T(n) = (p € S, (n, p) = 0), S(n) = (p € S,(n,p) 2 0), A be 
the Laplace-Beltrami operator, and a, b be the constants. 

One of the integral geometry problems is to determine the function f(n) € C^(S), 
provided the function f (p) is known. p € S, 


fora | find e f (5 2)f()a8 
P(r) S(n) 


the function f(p) is not necessary even, (see Section 3.2). 


Theorem 4.3. (Anikonov and Stepanov, 1991) Let a # 0, b Z 0, a b Z 0. The 


inversion formula holds 
fn) = z- f ()0 — (n,p)) aC — (n, »))AS 
nj -— an J P) — (n, p)) in n,p 
with 


p- Li f (A*2ep(n.pas 
$(n) 872 a (n, p? -t 


mS J Ag (p)dS 
472b (n, p) 
t=0 S(n) 


olo è _ (P) +p) 
PGF TEENS [eis a 2 ] 


4.1.4. Equation for function p(y) 


In this section we consider the generalized problem of the integral geometry, namely, we 
seek two functions f(x) € C™(|z| € 1), z = (41,22) and (v) € C™([0,27]), given the 
function 


f(ep- J f(z) e "omi singt cos 9) gy, 
z cos +y sin y=p 


It should be noted that if the function u(y) is known, the problem is linear and the inverse 
formula for the function f(z) holds (Kuchment, 1989). 


Theorem 4.4. (Anikonov and Shneiberg, 1991) The function p(y), 0 € p < 2m is the 
solution of the differential equation 


dy 
d 


| [en PP) | oypf(o,p)| e" dp = 0, UE 
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Other similar results are contained in the papers of (Anikonov, 1983; Anikonov and 
Shneiberg, 1991). 

In particular, the uniqueness class of two-dimensional integral geometry problem con- 
sists of the smooth functions f(z1, 22, y) such that 


LER cs T Ob 2G 
O02, P GEZI P ap f 


4.1.5. Support theorems 


Let f(P) be the Radon transform 
ft)» | fG)àn, 
P 


where P is a plane. 

Theorem 4.5. (Helgason, 1980). Let the following conditions hold 

1. f(z) € C*(R^), 

2. |z^||f(z)| € oo, Vk 2 0, z € R^, 

3. 3A > 0 such that f(P) = 0 for all P : p(0, P) > A, where p is a distance in R^. 
Then f(x) — 0, |z| > A. 


Theorem 4.6. (Logvinenko, 1988). Let the following conditions hold 
1. f(z) € C*(R"), f(z) = o(|z| *), Vk > 0, z > oo, 


2. the function 


fop) = f feu 


(w,2)=p 
satisfies the condition n € e, e **?, and c, > 0, €u > 0 are constants, 
3. JA > 0 and a set e € Ñ, ue $0 such that f(w,p) = 0, w € e, |p| > A. 


Then Support f(x) is a compact set. 


4.1.6. Inversion formulas for tensor fields 


Let f(z,£) = > fais -Eim c € R^, |£| = 1, m > 0 be a tensor field such that 
1=0 


1. f(z,£) € C*(R" x |g] = 1), 
lens 


2. 1 E) L———— 
Pup) z15,...0z8» 


rcR^ 


« oo, Vk 2 0, a; > 0. 
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We consider here the problem of the integral geometry to find f(z,£) if we know the 
function 


HOnE J f(a 4- t£, £)dt. 


Theorem 4.7. (Sharafutdinov, 1989) If the conditions 1 and 2 hold, then 
3. f =f e dv, 6f =0, 
4. 

f= (-ay? 


k=0 


Im/2] y 
Y. ali- A7 wf 


cy are constants, i, A, d, j, p are the operators defined in the paper of (Sharafutdi- 
nov, 1989). 


4.0 THE UNIQUENESS AND SOLVABILITY 


4.2.1. Integral geometry problem uniqueness of the solution 


Consider the two-dimensional problem of finding a smooth function f(r,y), 0 € r € 1, 
0 <y € 2m from tbe function 


s 8o; V? 
f e: [reste (P) de, ososi osos 
gj -a-(-1)?/r—pb;(r,) PEC'O<p<r<l), vlo p) = tlp, p). 


This problem is equivalent to the problem in the integral geometry where curves have 
group properties (Romanov, 1967). 


Theorem 4.8. (Romanov, 1967) If f(p,a) 20,0 € p € 1,0 € a € 27, then f(r,y) = 0, 
0<r<1lO0<¢<2z. 


The generalizations and applications of this result may be found in the book by (Romanov, 
1987). In particular, using the inversion formulas one can prove the uniqueness theorem 
when manifolds are not a plane or a sphere. 


4.2.2. Uniqueness of the solution of the integral geometry problems in the 
class of analytic functions 


Let the smooth manifolds B(z, p), z € R^, p > 0 satisfy the condition 


B(z,ppcQ*(zp)-(í((rzy:zeR' v»0, lz- z? +y? <p, n>1}. 
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Theorem 4.9. (Anikonov, 1978) Let the function f(z,y) € C*(z € R”,y 2 0) be 
analytic in the variable y in the domain z € R^, y > 0. If 


/ f(z, y)dp =0, z€ R^, p> 0, 
B(z,p) 


then f(z,y) - 0, z € R", y > 0. 


4.2.3. Differential identities and integral geometry 


The differential identities can be used well to prove the uniqueness, stability, and existence 
theorems of the solution in integral geometry. Here we consider only the two-dimensional 
case. Let M, M be the differential operators (see Section 2.9) 


MW = (2. --) cos 6 + (Z-a) sin 0, 
y 


Ox 
MW = SEE sin 8 — ELAN cos 0, 
Ox Oy 


where W(z, y, z), 0(z,y,z), c(z, y) (Oc/Oz = 0), d(x, y) (0d/0z = 0) are functions in C? 
in the domain z? + y? € 1, 0 € z € 2r. Here we give the application of Theorem 2.17 
when D is a circle. 


Theorem 4.10. (Anikonov, 1978) The following identity holds 


dne ð- 00 | (aw ?^ (aw d 


a (awaw) o(owow| a (aw gw 
Oy V Ox Oz Ox V Oy Oz 4 Oy ðr ]' 
z? +y <l, O<z< 2. 

In the case when c = 0 and d = 0 the identity is the same as obtained by (Mukhometov, 
1977). Let A(z, y), w(x, y), a(z,y), b(z, y), c(z, y), and d(z, y) be smooth functions in the 

domain z? + y? € 1, and 


o(z,t) = {z(s;z,t), y(siz,t), hag = 0082, phas = sinz} 


be a set of smooth extremals of the functional 
J [aaz + bdy + Ày dz? + dy?] , 
s 


which joins the points (cos z,sin z) and (cost, sint). Assume that the set y is regular. In 
this section we consider the non-linear problem of integral geometry to find the functions 
A, #, a, b, c, and d if the functions 


f(z,t) = J [ods + day + w/a? y. 


(zt) 


Integral Geometry 93 


vy 
6(z,t) = arctan “> 


s=0 


are known in the domain 0 € z € 27, 0 < t € 27. From Theorem 4.10 it follows. 


Theorem 4.11.  (Anikonov, 1978) Let the functions À, p be known for (z,y) with 

Oc ð , 
x?+y? = 1. Then the functions A(z, y), u(z, y), d and eb: are uniquely 
determined from the functions f(z,t) and 0(z, t). 


The problems in multidimensional cases were considered by (Romanov, 1978; 1987) and 
(Pestov, 1985). Pestov was the first to use the curvature of Riemann spaces in identities. 


4.2.4. Integral geometry and evolution equations 


In the case of a line integral the integral geometry problems reduce to Cauchy problems 
for evolution equations (Anikonov, 1969a; Lavrent'ev and Anikonov, 1967). After solving 
these problems, the unknown functions can be found by formulas. 

The two-dimensional problem of integral geometry reduced to the inverse problem for 
the transport equation: 


oW oW . oW 
ur PT npe K(z.y, 9) 7 = f(x,y), 
Wino = f(z.) rER, 0Syst7, y20, 


where the functions W(z,y,y) and f(z,y) are unknown. Information is given by the 
function f(z, y). The function K(z, y,y) is the curvature. 


Theorem 4.12. (Anikonov, 1978a) Let a smooth function u(z,y,z) be the solution to 
the Cauchy problem 


ðu _ f à Ou(z,y,p) Fz , ,Ou(z,y,P) 2 - 
"M Fr POR Zr ans Au 


ul y= = f(z,arcsin 2). 
Then 
Ou Ou 
flew) = [ZE + KG 


z=0 


REMARK . If the function K does not depend on y, we have 


u(z, y, z) = i^f. 
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4.2.5. On the solvability of a problem of integral geometry. 


We consider here the question of solvability of a integral geometry problem in the case 
when the integration curves are geodesics of a fixed analytic metrics. Let 


ds? = B(z,y)(dz? + dy?)) (4.1) 


be an analytic metrics defined on the (z, y)— plane. Henceforth, we suppose that B > 0 
and 0B/dy < 0. We denote by 7(£,y) the geodesic of the metrics (4.1) emanating from 
the point (£,0) at an angle y. Let 7(£,y) be the part of 7(€,y) lying in the half-plane 
y 2.0. The problem we consider is to find the function f(z, y) in the region y > 0 if the 


function 
Ae) = f| feas + dy? 
(Ep) 
is given for — < £ < 6 and 0 € p € a. 
We start with giving preliminaries. Denote by 5(£, v), n > £, a point of the line y = 0 
that belongs to the geodesic 4(6,), and let 6(&,) be the angle between y(£, p) and 
the z - axis at the point (5,0). Set a(£,q) = B(£,0)sin y, &(£,v) = B(n(é,¢), 0) sin y. 


Assuming 4 - 
t= / ei x Lid 
s 


the equations for geodesics of the metrics (4.1) become 


dr ðB dy ðB 


= MI = o n 2 _ p2 
dt? Oz’ d? Oy’ BENS [63] 


It is easily seen that any first analytic integral of (4.2) can be written in the form 
f = u(z,y,y)z + v(z,y, y) 
where u(z, y, y) is an even analytic function of y', and v(z, y, y’) is an odd analytic function 
Li 
i ‘ta the following we shall consider real analytic functions in a neighbourhood of the 
origin of the real Euclidean space R?. 


Let E be the set of all analytic functions g(z,z), z € R, z c R. We denote by T the 
set of analytic functions r(z, z) generated by the first integrals of system (4.2); namely, if 


f — u(z, qs y')z + v(z,y. y) 
is a first integral of (4.2), then 
r(z,z) = u(z,0, z) + v(z,0, z). 


Let E/T be the quotient space of the equivalence classes: g1 ~ g2, 91 — 92 € T. Let A be 
the set of functions f(E, y), with 


HE) = 5 lem b) + Bln, 0) cose) — 9(7,—4)(1 — B(n,0) cos 8) 


—5 [a(6,a)(1 + B(E,0) cos e) — g(€, -a)(1 — B(E,0)) cos e], 
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where 7, b, a, 0 are the above defined functions, g € E. Denote such a representation of 
the functions f(é,y) by Ag and write f 2 Ag. 

We denote by Dg, g € E, the analytic function f(z,y) of the variables (x,y) defined 
as follows: we consider the system of the equations 


Qu Qv OB 
2 - far - 2(k 1)B cue, 
Qv, Qu, ðB Ours 2 OB 
—_— = Bo - B 2k + 3)u.-—, 
ðy ðs "Pg. 9 CEPI 
with the Cauchy data 
url o = us (zx), Uely-o = vi(z), 
such that 
3 [(u? 0 2% + yp] = g(z,z), gc E. 
If this problem has an analytic solution u,, vz, k = 1,2,..., then by definition 
OB OB | Quo 
Dg = w iure 9.5 


Theorem 4.13. (Anikonov ,1978a) 


1. A necessary and sufficient condition for the solvability of the integral geometry problem 
in the class of analytic functions f(z,y) is that the function 


Fee) = f fas + ay? 
KE) 


can be represented in the form f = Ag. Moreover, for any function g(z,z) the 
operation Dg is defined, and 


Ag — f Dgyj dz? + dy?. 


(ép) 


2. Ag, = Ag if and only if the equivalence classes of E/T corresponding to the functions 
gı and g coincide. 


3. If 


J fide? + dy? = J fo dz? + dy? 
Y Y 
and if fi, i = 1,2, are analytic functions, then fi(z,y) = fa(z, y). 


We shall give the proof of this theorem in Section 4.5. 
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4.2.6. Integral geometry and the structure of Riemann spaces 


Let M be an n - dimensional (n 7 2) simply-connected compact Riemann manifold of the 
class C*? with a metric tensor g = (gij) and a strictly convex boundary 0M. 

We introduce the following notation: R;;u is the curvature tensor of the metric g, T; 
is the tangent space at z, QM is the tangent bundle of unit vectors, QM = ((z,£) : x € 


M, E € Tz}, 
él = y sux) 


(here and below summation from 1 to n over repeated subscripts and supscripts is under- 
stood), 
QƏM -((z,£):z€8M, EET, \El=1, (6v(6)x0), 

where v(x) is the outward normal to QM at the point z, and Yz (t) is the geodesic defined 
by the initial data yz ¢(0) = z, Yz (0) = E 

The set of the smooth (class C??) covariant symmetric tensor fields of valency m, 
m > jon M and OM is denoted by $,, M and S,,0M, respectively. Define the operation of 
symmetric covariant differentiation d$, M — Sm+;M, d = oV, where ø is the symmetric 
operator and V is the operator of the covariant differentiation with respect to the metric 
g. Henceforth we shall assume that M is the dispersing manifold: any geodesic of it leaves 
a compactum K C M, K 8M = 0. We define the function t(x, £), (z, £) € QTM as the 
length of the geodesic 44, (r, £) € QTM, and call it the hodograph of the metric g. By 
the assumptions we have made about M, the hodograph t*(z,£) is well-defined on N- M. 

Suppose that f = fi,..;,, and f € S4, M. 

Consider the integrals along the geodesic joining points of 0M 


Ü(z£) 
= | fole O9. nat. (43) 
0 


Theorem 4.14. (Pestov and Sharafutdinov, 1987) Let M be a compact dispersing 
manifold with a strictly convex boundary OM, and suppose that the sectional curvatures 
of M are nonpositive, that is at each point x € M, Riju£'r/£*n! < 0 for any 6,9 € Tr. 
Then the integral (4.3) vanishes when (z,£) € QƏM if and only if 


f=dv, v€S,4,M, v|lj,4,-0, m21, f=0, m=0. 


It turns out that such theorems of the integral geometry enable us to determine the 
structure of the Riemann space from the integral information. 


Theorem 4.15. (Anikonov and Pestov, 1989) Suppose that the conditions of Theo- 
rem 4.14 are satisfied. Then the equations of geodesics admit a first integral 


Jod) jt) = t ig i 
where u = (u1, u2,..., Um) € SmM if and only if there is a field u € SmOM such that 


uj im (EET E = uu Qe (2,6) 0240 (2,6)... 325 (0 (0,0), (2, €) € RAM. 


We shall give detailed analysis of this results in Section 4.4. 
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4.3 SOME APPLICATIONS 


We consider a multidimensional inverse kinematic problem and obtain the necessary and 
sufficient conditions in order that some differential relations for unknown functions are 
held. 

The proof is based on the theorems of the integral geometry. One formulation of 
the multidimensional inverse kinematic problem consists in the following: in the ball 
w = {z : x < 1) of the real Euclidean space of variables z = (z1,...,24), n 2 2 we 
consider a 3 times continuously differentiable metric 


ds? = A(z}? |dz|’, (4.4) 


such that any two points p € w, q € w can be joined by a unique geodesic (p,q) of the 
metric (4.4); the function 


Wo - f Xeézl, pet dct 


(p.a) 


is assumed to be known, and it is required to find A(z) for |z| < 1. 

This problem has an important physical interpretation in the seismology. The function 
A(z) = 1/v(x) is the inverse to the propagation speed of a perturbation. The function 
W(p,q) is the time this propagation takes along the geodesic (p,q). The formulated 
problem consists in determining the velocity v(z) on the basis of the times W (p,q), 
lol = 1, [e| = 1. 

Some results on the uniqueness, stability and existence of the multidimensional kine- 
matic problem solutions and of the integral geometry closely related to the inverse kine- 
matic problem are contained in the papers of (Anikonov, 1978a; b; Romanov, 1987). The 
constructive techniques to solve these problems are particularly interesting in some ap- 
plications. Here, necessary and sufficient conditions are presented so that the differential 
relations for the known function W (p, q), |p| = 1, |g| = 1 and the unknown function A(z), 
|z| € 1 are simultaneously satisfied. We proceed to formulate the results. Suppose that 
the components u(r), k = 1,2,...,n of the vector-valued function u(x) satisfy the system 
of equations 
_ Qu, = Ou; 
Oz; = ~ Ox,” * Oz, = Oz; 


In the case for n = 2 this is the system of the Cauchy-Riemann equations; for n > 2 its 
general solution is given by 


(4.5) 


u(z) = Az — Bx? + 22(B,x)+C 


where A is a constant matrix with the elements Aj, such that Aj, = —A,j, k £ j, 
Aj; = Ao, B and C are constant vectors, and (B,z) is an inner product. There is a 
considerable difference between the classes of the vector-valued functions u(x) for n = 2 
and n » 2. 

Let u&(z) = (ug ,...,ux,) and v(x) = (vi,..., vk.) be vector-valued functions satis- 
fying (4.5) for any k. We set 


n 


1 
Gij = 2 D (us, vi, 4 unv) " 
k=1 
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and let 


n 
aot 
Q= Y, ajziz'. 


ij=l 


Along any geodesic y we have 


a 2 dos (v, LA T m) z', |z|-2 


k=1 
where 
eer Cerro en aes 2» 2 
Bvr; 


LY = A(z) (grad À, jee S Jz, 


Consider the differential expression 


7 Ow Ow ðw ðw 
Dur 2. (wo. 2e) (ue. 32) = (uo. a) (woo. a) , 
where 


0 0 0 0 0 ð 
w(p,q) = J A(z) |dz| , xxx) Fe (2. 2) 


(p.a) 
and u(z) and v(z) are vector-valued functions satisfying (4.5). 


Theorem 4.16. (Anikonov and Pestov, 1990a) In order that 


Lw = 9(q)—(p), ld|-!, [pl =1 


for some continuously differentiable function y(x), |z| € 1 it is necessary and sufficient 
that 
0p _< . 
On 2 (v, DKA ug bjr), im 1,2... n 


To formulate another result in the inverse kinematic problem we suppose that u(x) satisfies 
(4.5) and has no more restrictions. Consider the differential operators Lw, LA such that 


ð ð 
Lu = (Zuo) + (Feu) + Bw, 
LA = (Zua) + * sp = + BA, 


: Qu , . Qu : 
where B is a constant, sp Jz is the trace of the matrix oe k,j =1,2,...,2 
T Tj 
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Now define M setting Mw be equal to the following determinant of the order n + 2 


0 0 OLw Lw 
* o % 
1 n 
0 0 — — 
lel lal 
" OLw p 
U-—|Óm pl 
O?Lw 
OpOq 
OLw Pn 
Op. |p| 


Lw 
Opà 
Let T(t) denote the gamma function, dp, dq denote the surface area elements of the 
sphere Ow at the points p and q respectively, and let 


is a matrix of the order n. 


where 


IF] =sup|F(q)|, «€ D. 


Theorem 4.17. (Anikonov, 1978a) The following inequalities hold 


| LA 


A 


, 
C(8wx dw) 


aac - CE J | Moda) 


From Theorem 4.17 we obtain the table 


Clw) - | w 


J d"-2(EA)?de 


ü = (0,0,1) = k, v = F(z1,23), 
ü = [k, z], v = F (x? + 22,23), 
u-r, v-tar( 3 
- zi +z? 2 2 a -aoarctan(z2/21) 
= aor + [k, z], v= 23F "owe INL ae ’ 
a = ok 4 [k e d ic tan =) 

+ [k, x], v T? + 13, £3 Gaal ; 
a = —k |z|? + 2z(k, 2), = |e)? F( = i Z) 

je Je 

“= 


- kla? ems). | v= fated (2 aoza + lel s, 
yz? +23 


a = —k |x|?  2z(k, z) + ak, | v = yz? + oF z letta 
ty’ atta 
where ù = u(z), n = 3, v(z) = 1/A(z), and F is an arbitrary function. 
We consider a certain inverse kinematic problem connected with geometry as well. 
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Let R^ be the n— dimensional Euclidean space of variables x = (21,22,..-,2n), 
n = 2,3,..., and Q be the ball |z| € R, with the boundary Q. Let B be an m— di- 
mensional differentiable manifold embedded in the ball Q with the local coordinates 
u = (uj,u5,...,u4) and the corresponding local parametrization z = z(u), |u| < 1, 
m &n. 

Consider the Riemann metric 


ds? = M(z)(dz? +... + dz?) = \*(z) dz]? 


on the ball Q and let p denote a point of the sphere OQ. ds? induces a Riemann metric 


d3? = L 35 Aij(u)du;du; 


t=1 3-1 


on the manifold B. Here we set 


in the equality ds? = A?dz?. If some integral information linked with manifold B is given 
on the surface of the sphere 0Q, then the problem of investigating this manifold B with 
induced metric given above is of interest for geophysics. Such an information may be 
given as a function 


A|dz|  a(z(u)) = 9(u, p) 


oY(z(u).p) 


of the variables u and p, where 4(z(u), p) is a geodesic of the metric ds? = A? |dz|? joining 
the points z(u) € B and p € OQ. Here the functions a = a(z) and A = X(z), as well as 
the parametrization z = z(u) are unknown a priori. The given function g(u, p) may, in 
general, be not even continuous. 

Before we formulate the results of this problem, let us give a physical interpretation 
of the terms above. The ball Q is the Earth, the manifold B C Q is the set of earthquake 
hypocenters, 1/A(z) = v(x) is the propagation speed of the solid waves in the Earth, a(x) 
is the starting time of an earthquake at the hypocenter z — z(u). The given function 
g(u, p) is the "indexed" time of the arrival of the disturbance on the surface of the Earth, 
registered by seismometers, where u is the “index” of the earthquake, and p is the “index” 
of the seismometer. 


Theorem 4.18. Let the function \(x) be twice differentiable and suppose that through 
every pair of points zı and z3 in the ball Q there passes a unique geodesic (21,13) of 
the metric ds? = X |dz|’. 

Then the function R(u, pı, p2) = g(u, pı) — g(u, p2) of the variables (u, pi, p2), |u} < 1, 
p; € OQ, is continuously differentiable and the following formula holds: 


From Theorem 4.18 it is clear that the intrinsic geometry of the Riemannian manifold B is 
determined by the function g(u, p). This leads to corollaries of an interest for applications. 
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COROLLARY 4.1. Let the dimension of the space R” be three, and B be a closed convex 
level surface of a function À = A(x); that is, B = (A(z) = c). Then the surface B is 
uniquely determined by the function g(u, p) to within the similarity and position in the 


ball Q. 


COROLLARY 4.2. Let m = n and let B be a simply-connected region with a smooth 
boundary; then both the function à = A(z), s € B, and the region B are uniquely 
determined by the function g(u, p) to within the conformal transformations of the region 
B. 


Given additional hypotheses on B and the function a(z) we may somewhat strengthen 
Theorem 4.18, namely 


Theorem 4.19. Let B be a simply-connected n-dimensional region with a smooth bound- 
ary and let the function a(z) = 0, x € B. Then the elements of the matrix c;;(u) inverse 
to Aj;(u) satisfy the relationship 


n n Og Og 
2-:2- uou: zx cu(u )= 


t=1 j=1 


for all p € OQ and, therefore, are uniquely determined by the set g*(u) = g(u, Pk) of values 


of the function g(u, p) at the points p = py, € OQ where the determinant (2 2) (5 ) 
ui uj 


is nonzero, k = 1,2,... n(n +1), i,j =1,2,...,n 


Here the function A(z), x € B and the region B are unique to within the conformal maps. 


4.4 THE STRUCTURE OF RIEMANN SPACES AND PROBLEMS OF 
THE INTEGRAL GEOMETRY 


The aim of this section is to formulate and prove the necessary and sufficient condition 
for the resolvability of Riemann spaces, and the existence of certain structures in these 
spaces, using the integral information associated with the geodesics. We give a detailed 
analysis of the results reported in section 4.2.6. 

The proofs are based on the differential identities and theorems of uniqueness and sta- 
bility of solutions to problems in integral geometry. The results can be used, in particular, 
by determining the Riemann metric by the lengths of geodesics which connect the bound- 
ary points of a smooth Riemann manifold. This problem is sometimes called the problem 
of nonlinear tomography. Its solution can be appreciably simpler if the Riemann space 
has an additional structure, e. g., a primary integral of geodesics. We shall establish the 
criteria for the existence of a primary integral which is a homogeneous polynomial with 
respect to the vector tangent of the geodesic. The existence of this integral is equivalent 
to the conservation of some symmetric covariant tensor field by the flow of geodesics, or 
(which is the same ) to the identity du — 0, where d denotes the symmetric part of the 
covariant derivative. Besides, we shall consider the symmetric tensor fields which satisfy 
the equation du = o(g ® g 9... Q g) where o denotes symmetrization, g is the metric 
tensor, and & denotes the tensor product. 

The criteria for the existence of these fields on manifolds with a convex boundary of 
nonpositive curvature are given in terms of the lengths of the geodesics connecting the 
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boundary points. In conclusion we shall formulate a criterion for the zero Gaussian cur- 
vature of a two-dimensional Riemann metric; it is based on the so-called turning number 
of geodesic which connect the boundary points of a circle. 


4.4.1. Designations, preliminary considerations and results. 


We assume that M is a smooth n— dimensional Riemann manifold and T' is a tangent 
fibre bundle for M. Let us denote the points of T by (z,é), where z € M and £ € T, (T; 


is a tangent space at the point z). We use coordinate systems 


(mU, 2l ee UE) 


on T, where U C M is the domain of definition of the coordinates (U, z!,...,z") on M, 
T is a projection T — M, and £' are the coordinates of the vector € € Ty with respect to 
the basis 0/Az',i = 1,...,n 

Let |£| denote the length of a vector £ € Ts, |El = \/gi;(x)é*€?, where gi(z) are 
covariant components of the metric tensor. We use the conventional rule which implies 
summation for recurring subscripts and superscripts. Let us consider the submanifold Q 
of unit vectors in T: 


Q= {(2,6) ET: EEN}, 


where Q, is a unit sphere at the point z, and use on Q the same coordinates as on T, 
bearing in mind that |£| — 1. 

We shall introduce some objects of tensor analysis on the tangent fibre bundle for the 
manifold M (see (Pestov and Sharafutdinov, 1988)). 

A tensor of type (r,s) on the manifold T' is called a semi-basic tensor if it can be 
represented as 


ai, 0 ð : 
uus EET 8...8 IE @ dz? Q... @ dz? (4.6) 
in each coordinate system (z^!U, z!,...,z^,£1,..., £^) on T. Here uj JA € C”(m™U). 


The coefficients in (4.6) are supposed to be real functions. By B; we denote the fibre 
bundle of (r,s)— tensors on T, and C™(T, B7) denotes the set of sections in this fibre 
bundle, i. e. the set of smooth semi-basic (r, s)— tensor fields on T. If N is a submanifold 
in T, we denote the restriction of C^(T, B7) on N by C™(N, B;). In the scalar case 
r = s = 0, and we use the designation C?(N). We can extend the conventional algebraic 
tensor operations, such as addition, tensor multiplication, or convolution, to semi-basic 
tensors. The Riemann metric determines the operation of changing from subscripts to 
superscripts and vice versa, and allows us to introduce a scalar product 


(u, v) = Ui p orm (u, u) = lul, 


for tensors of equal valency m — r 4 s. 


h v 
The differential operators of the first order V and V: C™(T, B;) ^ C™(T, B;,,) can 
be defined in the coordinate form as 


O= up ") (9 (9 
ur me Q, veo Ve] - TM ae 


where I7, are the Christoffel symbols of the metric tensor g;;; Vit; can be found by 
the rules of covariant differentiation for the ordinary (r,s)— tensor fields on M. The 
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h v 
operators V and V are called the vertical and the horizontal covariant derivatives (see 
(Pestov and Sharafutdinov, 1988)), respectively. Hereafter, we shall use these terms. The 
following formulas for these operators can be directly verified: 


v v v v v h h v 
Ve Vi- Vi Vx =Vk Vi- Vi Vi = 9, 
v h h v v 
(v. Vi— Vi Vi) u= -Rt V pt, uc C*(T), (4.7) 


v h v : h $ k v ; - 
V&gij —-V&gj —0, Vk; =V; =0, Vi =0, VIE =ó, 


where (R$) is a curvature tensor. We define contravariant derivatives 


h i ij h v; v 
V'-g" Vj V'-95 Vj 
where g are the contravariant components of the metric tensor. 


h 
The operator V can also be considered on the fields which belong to C™(2, B7), 


h A 
because V ,|£| = 0 and, consequently, V ;, i = 1,...,n, are vectors tangent to 2 for 
al (z,£) € Q. Besides, there are tangent vectors P, on Q, which are defined by the 
decomposition 
v 1 ig 
Vic pp iR+P, R=EV;, P= 


where £; = gij£?, and hence the operator 


1 


git 6 Vi A V3), (4.8) 


P : C”(N, Br) 2 C°(2, Bza), (Pu = (Piu) 


is defined correctly. Let P? = g;;P;. 

The ordinary tensor fields on M can be identified with semi-basic fields whose com- 
ponents are independent of £. We shall consider symmetric tensor fields on M, i. e. 
(0, m)— tensor fields whose covariant components are symmetric about subscript trans- 
positions. 

Let Sm be a fibre bundle of symmetric tensors on M with valency m, and C*(M, Sm) 
be a set of smooth symmetric tensor fields on M. The restriction of C?(M, Sm) on OM 
is denoted by C? (0M, Sm). 

By d: C®(M, Sm) > C*(M, $441) we denote the symmetric covariant derivative (see 
(Sharafutdinov, 1984)), d = c V, where o denotes symmetrization, and V is the operator 
of covariant derivative. We define the divergence 6 : C*?(M, Sm) > C*(M, S, 1), m2 1, 
as 


(6u),. ia — g" V kliima uc C"(M, Sim): 


The operator éd is a natural generalization of the Laplacian for symmetric tensor fields, 
in particular, the Dirichlet problem 


du = f, ulam = Uo 


has a unique solution for any uo € C*(0M, Sm) and f € C®(M, Sm) (see (Pestov and 
Sharafutdinov, 1988)). 
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For symmetric tensors u and v we define symmetric product uv = o(u @ v), and 
denote the symmetric power of a tensor u by 


u” = o(uQ@u@...@u). 


Let u € C®(M, Sm), m 2 1. Then < u,£€ > will denote a function on T (or on 2 ) which 
is defined in coordinate notation as 


< u,t >= uj, im (z)£^ t m, € € Tz, E € Oz. 


Let H = ra Y be the geodesic flow of a manifold M (see (Godbillon, 1969)), and let ó^ 


and 6" be horizontal and vertical divergence, respectively, 
hp ey ip LX al 
Su=Viu, ó6"u-V;u. 
Of crucial importance for the proof is the following identity. 


Lemma 4.1. Let w € C™(T). Then the identity 
vul ? — Rautt* Vow Vw + é(VwHw) = (V + 2w V Hw) - 20 (V Hw) (49) 
is true for a semi-basic vector field V which is defined by 
Vi 2V j£ Viw- € Vw) 
PRoor. By definition, A 
EV iw = Hw. 
As a result of the action of the operator Viw Vj we have 


h 2 i A. v h h. v h v hee 
vw -Ë V’w Vi Viv -V/v VjHw = & (w V Hw) — wé'(V Hv). (4.10) 


Let us rearrange the second term on the left-hand side of this equation. According to 
(4.7), we obtain. 


i A. v h i da h v h i h. Di v ; A. h 
2€ Viw V; Viw = 2& Viw Vi V jw Vi(£ V^v Vu) VË V^v Viw) 
2 i h h Ve 
—-€ Viv V; Vu. 


u; h h h h io; h h h 
-E Viw(V: Vjw— Vj; Vi)w -Ë V^wV; Vw -|Vw 


Here the first, fourth and sixth terms form ó^V. Taking commutation formulas (4.7) into 
account we find 


2 


A. v h h " v. v h 
26 Viw V; Viw = V é"(Vw- Hw) + Rag£t* Viw V'w — vw 


Substituting this expression in (4.10) we obtain (4.9). 
Let à € R!. We define a differential operator 


P, VP, 4 AH, 
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and an operator of the homogeneous extension of power À 
rw =KePu(2£), tem. £40 vec. 


We obtain a similar identity for functions in C®(Q), as a corollary of Lemma 4.1. 


Lemma 4.2. For each smooth function w(z,£) on Q 


: A hx 
Iul os Rj n€' "Pi wP'w + Pi(Hw V'w) + (24 + (Hwy 


: (4.11) 
= {V +2w(VR,Hw)| ) -2uP,Hw 
Q 


where 


Vi =Ẹ ;w(£ Piw — £ Piu). 


ProoF. The proof consists in applying identity (4.9) to a function Hw and using (4.8) 
(we should also take into account the properties of the curvature tensor and relations 
P,R, = RyP;, RR, = AR). 


Let M be compact oriented manifold. The differential form 
dT = gdz ^ d£ = gdz! ^... ^ dz” A d£! ^... A d£", 


determines a volume element of the manifold T, where g = det(gi;). A volume element 
dQ in 2 is defined by 
dT = d |f| ^ dQ, 


or 


dN = gdz ^ dhs, 


where df), is a volume of a sphere N, and dT, = d |£| A df),. In the coordinate notation 
we have 


dO, = J/g(£'d£? ^... A d£" —...-- (71) £d£ A^... Ad£"7), EE DG. 
Let us define on Of. a differential form 
dE=dSAQ2,, zc€0M, 


where 
dS = Jg(v!dz? ^... A dz" — ... + (C1) *! v^dz! A... A dz?!) 


is a volume element of a boundary 0M. The form dE is naturally referred to as the 
volume form of the manifold dQ, because 


dQ = —dp A dE, 


where p(x) denotes the distance from the boundary 0M. The function p(z) is smooth in 
the neighbourhood of the boundary, and Vp — —v on 0M. 
In the sequel we use the following integral relations. 
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Let u be a semi-basic vector field on T. The Gauss-Ostrogradskii formula, 


[ 50 - [ (v yas, (4.12) 
2 an 


holds for the horizontal divergence (see (Pestov and Sharafutdinov, 1988)). The formula 
j P,u'dQ, = (n — 1) / (é,u)dQ, (4.13) 
f, f. 


follows from the equality 


J véd, = J vdT,, v € C”(T,). (4.14) 
Qs feci 


Indeed, we have 


Je (ues -— ugs€i) dQ, = J [G9 - (we), dT, 


ds kei 


zinc) J ujdT, = (n — 1) f Gude, 
\éls1 Ns 


and hence we obtain (4.13) for the vector field. 
The operator of integration over the sphere Q, transforms the semi-basic tensor fields 


on Q to ordinary tensor fields on M. In particular, let us consider the symmetric tensor 
field on M: 


J E... Em dN. 
Qz 
We shall know that the equation 


JEn sce yon (4.15) 
Qs 


is true, where g™ the m— th symmetric power of the metric tensor g”, 


m LITTLE 
g” = o(g @8®...@g) 
and 
1:3....- (2n - 1) = 
~ n(n+2)...(n+2m—2) ™ 
where c, is the the volume of unit sphere in R”. We prove (4.15) by induction. 


Let m = 1. We easily see (e.g. in Riemann coordinates in the neighbourhood of a 
point z) that 


| 6640. = gë, 
f; 
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where the constant c can be defined from convolution 


On 
[$90 = cn, c=q=—. 
n 


Suppose that (4.15) is established for some m. Then, according to (4.14), 


T fe Eae ea dite 


f 


* J [6 KE BIA gena +62 E e pe peus tec D : MO EHI dTz. 
USES 
Taking into account the homogeneity of the terms in the right-hand side, we rearrange 
the latter equation: 


[Gams AM E pits meme fme] an. 


Changing from subscript to superscript and using the induction assumption, we get 


[e emm. 


E [g^ (g mya- d2m41 + gitana (gm damen] 
15.4 
- Cm (gt ys 32m42 


Suppose that M is a compact dispersive manifold, i. e. a compact oriented manifold 
which does not contain infinitely long geodesics and has a strictly convex boundary 0M. 
Then, we can define on Q a function t(z,£) as the length of a geodesic *(z,£,t) which 
emerges from a point z € M in the direction of a vector € € Q. Here 2 is a compact 
manifold with boundary 02 = 0.0 N 0,€, where 


0.0 = {(z,E)EN:2EOM, (6v(z) <0}, 
Q$-—((r0)€e9:z€08M, (6&v(z)) 2 0). 
and v(z) is the unit outward normal of QM at the point z. We also designate 
09 = ((z,£) En: (£, v(z)) = 0}. 


Obviously, 092 = 0(0.€) = (34N). By definition, t(z, £) = 0 for (z, €) € 34N. The 
restriction of the function £(z,£) on -N will be denoted by to(z,£). We define a mapping 
vy:0.02 9.0: 


v(z,£) = {7 [z, €, tolz, €)] ye [z, €, to(x, €)]} : 


As above, let p(x) be the distance from the boundary 0M. In the region where p(z) is 
smooth we consider the following decomposition similar to (4.13): 


Vp = (Vo, V)Vp4 Q, 
where 
" h h 
Qi = V"o(Vip Vi — Vip V x). 
The main properties of the function t(z,£) are stated in the following lemma. 


Lemma 4.3. Let M be a compact dispersive manifold. Then 
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1. t(z,£) € C^ (0X8801), 


2. Ht = —-1, (4.16) 

3. to € C^(0..Q), 

4. v is a diffeomorphism 0_Q to -N and v o = id, where id denotes the identity 
mapping, 


5. the derivatives P;t and Q;t are bounded in 2Q\02. 


Pnoor. The following proof of statement 1. is given in (Pestov and Sharafutdinov, 1988). 


1. Let a set Mo be defined by 
Mo={xEM, p(z)2 po po > 0}. 


The function p(x) is smooth in a certain neighbourhood of the boundary. Let us 
assume that po is so small that a set M\ Mo can be considered as this neighbourhood; 
it is evidently a smooth manifold. Hereafter we consider p(x) only on M\Mo. The 
function t(z,£) is a solution of the equation 


p (z, €, t(z, €))) = . 


Suppose that cos o(z, £) corresponds to the angle between a geodesic and 0M, i. e. 


cos (s, £) = -EED 


where t = t(z,£). Since y(z, 6, t) is a smooth curve which smoothly depends on the 
initial data, then, according to the implicit function theorem, t(z,£) is a smooth 
function for all (z, £), where cose(z,£) # 0 . This inequality holds for all (z,£) € 
f Nàof1, provided that the boundary is strictly convex. 


2. Let (z,£) be an interior point in the manifold 2. We define the operator hy of trans- 
lation along the geodesic flow: 


hz, €) ES (y(z,€,t), 4(z, €,t)), -t(z, —£) sts t(z, £). 


Direct calculation based on the equation for geodesic, 7 + rey, leads to 


< —(uoh)-2(Hu)oh, vcC?(Q). (4.17) 


Applying (4.17) to t(z, £) and assuming t = 0 we obtain (4.16) (we take into account 
that t(hi(z,£)) = t(z, £) — t and ho = id). 


3. Let tg denote that part of the function t|,. which is an odd function in £. It is 
evident that tg = 19/2 on 0. f). We are to show that tg € C(O) and thus prove 
statement 3. 


Let us consider the function 


u(z,6,t) = A(r(z,,t)), (26)€e00, —t(z,-6) < t < t(z,€), 
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where A(z) is a smooth extension of p(x) to the entire manifold M (p(z) is defined 
only on M\Mo). Integrating the identity (tà)! = ù + t Ë from —t(z, —£) to t(z,£) 


we have 
to (z.£) 


25r Oi(n62(0)- f tina. 
0 


Changing the integration variable and dividing by 2tọ we obtain 


1 
F(z,£,2t5) = i(z, 6215) — 2t; | sis, 6,2stz ds = 0. 
0 


It follows from 1. that the function tlag is smooth on 0Q\0.Q. Therefore, it 


is enough to establish that £j is a smooth function in the neighbourhood of any 
point (zo, ĉo) € Oo. The function F(z,£,t) is defined and smooth for (z,£) € Q, 
—t(z, —£) € t € t(z,€) and satisfies the conditions 


F(zo, £o,0) = —(v(xo), ĉo) =0, 


1. ieee 
Fi(zo, £o, 0) = z (zo, 60,0) = z S00 ViV jp(zo) « 0. 


The latter inequality implies that the boundary 0M is strictly convex. According 
to the implicit function theorem, the smoothness of F implies that tj is also smooth 
in the neighbourhood of (zo, £9), and consequently to is smooth on 0f). 


4. The fact that y% is a bijection follows from the uniqueness of the solution to the Cauchy 


problem for geodesics, and 3. implies that the mapping V is smooth. Equation 
ab ow = id is evident. 


5. The ratio to(z,£)/ cosc(z,£) is bounded on 0_2\0)N (Pestov and Sharafutdinov, 
1988). This implies that the function £(z,£)/ cos(z,£) is bounded in Q, where 
Q is the interior of the manifold N. 


Let (z,£) EN. Consider a point (z1, £1) € 0_Q with: 
Tı = q(t, £, —t(z, —£)), & = —YX(z, É, —t(z, —é)). 


We have 
t(z1,6) < 
cos y(21, £1) 
but cos o(z1, £1) = cos y(x, €) and t(z, £) < t(21, &1), therefore 
t(z, €) 
cos pe, 6) <* 


Similar to 3. we introduce a function 
u(x, &,t) = p(y(z, €, t)). 


Here we consider this function on the set (z, £) EN, 0 € t € t(z, £). Integrating the 
identity 
(ti) = pt tà 
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from zero to t(z, €) we obtain 


t(z,€) 
—#(2,€) cos o(2,) =—A(z) + | tà(s 6,t)at. (4.18) 


0 


Since t(z, £) € C^(Q), it is enough to verify statement 5. for x EM \Mo, where M 
is the interior of M. For these z we can replace f(z) in (4.18) by p(z). Let us apply 
the operator P; (or Qi) i = 1,...,n to (4.18) (note that P;p = Qip = 0 ): 

— cos p(x, €)Pit(x,€)  t(z, €) [(Pin)(z, € t(x, €)) + (x, 6, t(2,6))Pit(z, £)] 


= tz Dis Gt )PaGs 0). f (Pto). 
0 
This together with the boundedness of t(x, £)/ cos p(z, £) imply that the derivatives 


P,t(z,é) and Q,t(z,é), i = 1,...,n, are also bounded on Q ( in any coordinate 
system). 


4.4.2. Formulation and proof of the main results 
We assume that M is a compact dispersive manifold and 
u E C®”(M, Sm) m2l1. 


Calculations lead to 
H < u, >=< du, >, (4.19) 


or, taking into account (4.16) 


d 
át < u(y), 7 >=< du(q), 3», (4.20) 


where 7(2,é,¢). For du = 0 this implies that < u(y), Y > is a primary integral for the 
equation of geodesics. In particular, for (z,£) € 0_N, t = t(z, £) we obtain 


< uo, >= (-1)" < uo, € > oy, 


where uo = ulay. Now the question arises whether there exists a smooth symmetric field 
u of valency m > 1 on M such that du = 0, provided 


< uo, $ >= (-1)" < uo, € > ov 


on OQ. €) for field uo € C”. We answer this question, as well as a similar question for the 
equality du = g*, positively. This is one of the main results of this section. 


Theorem 4.20. Let M be a compact dispersive manifold with nonpositive curvature, 
1. e. 

Rijk n Cy! <0 
for all z € ME, n € T, and let u € C®(M, Sm), m > 1. Then the following statements 
are equivalent: 


Integral Geometry 111 


1. ódu = 0, < uo, Ê >= (-1)* < ug, Ê > op, uou 8M? 


2. du = 0. 


Theorem 4.21. Let M be a compact dispersive manifold with nonpositive curvature, 
and u € C®(M, Sox-1), k > 1. Then the following statements are equivalent: 


1. Sdu=0, to+ < ug >+ < ug £ 2 0$ 20, uo — ula, 

2. du = g*. 

Theorems 4.1 and 4.2 are corollaries of the following result, which generalizes the main 
result in (Pestov and Sharafutdinov, 1988). 

Theorem 4.22. Let M be a dispersive manifold with nonpositive curvature, and 


t(z,€) 
w(z,8)= f fü(z £t, 


o 
where f € C?(Q) and P) f = 0 for some à > c. Then 


wo = w|y a € C”(8-N) 


and 


PES Den ' [Pwo Quo) + 2w (v, V Ra f) 


| dE, (4.21) 


where 


(Pwo, Quo) = (EP wo — £&P*wo)v;Q;wo. 


PROOF. The function wo is smooth because to is smooth (Lemma 4.1). We shall show 
that w satisfies 
Hw - -f (4.22) 


in Ñ . Let (z, €) €f. Then 


t(z,€) 
t(he(2,€)) = (26) -r wh) m. f fle(2,€))dt 


T 


for all 7 in the interval 0 € 7 < t(z,£). If we differentiate the latter equation with respect 
to r, put 7 = 0 and take into account (4.17), we obtain (4.22). 

By Qo we denote the tangent fibre bundle for the manifold Mo introduced in Lemma 
4.3. The function w(z,£) is smooth on No because t(z,£) is smooth. In particular, the 


derivatives P;w and Qiu, i = 1,...n are bounded in Q in any coordinate system. Let us 
apply (4.11) to the function w|p,. With regard to the theorem conditions this leads to 
the inequality 


vul 24 PHw- Vw) + (24 + 1)(Hw)?  & (v + 2w (v RAHw) 


) > (433) 
Qo 


112 Chapter 4 


where v' Sud (€ Pw — £&P*w). Integrating (4.23) over No we take into account (4.12), 
h 
(4.13) and the obvious inequality [vu]? > (Hw)? to get 


| d Es, 
flo 


where vo is the outward normal to @Mo, and d Eg is a volume element of the manifold 
No. Let us show that (4.21) can be obtained from this inequality for po — 0. 
Let 


1 " 
EE Dandi J [iew Qu) +2w (vo, V Raf) 
Mo 


ON = ((2,£) € ON: (E vo) | <e}, 0<e<l. 
We have 
J (Pw, Qw)dE = J (Pw, Qw)dE + f I(Pw, Qw)d E. 
85g BeNo INNAN 


Since w is smooth on 2\0oM, the second term tends to 


l((Pw,Qw)dE, 
8-(A8,0 


when po — 0 (recall that w|;., = 0 by definition). As for the first term, its absolute 
value is bounded by ce with constant c independent of po, because the derivatives P;w 


and Q;v, are bounded on 9. Hence, we see that 


J (Pu, Qw)dEo > f (Pwo, Quo)dE 
aNo aln 
for pg — 0. The theorem is proved. 


Proof of Theorem 4.20. The fact that 2. implies 1. was proved before (Theorem 4.8). 
Suppose that 1. is satisfied. Let us consider a function 


t(z,£) 


w(xé)= f fic. 0)8t 


0 
where f =< du, >. The function satisfies the condition of Theorem 4.22, namely 
Pntif =< du, € >= 0, for à = m + 1. Then, from (4.19) we have 


w(z,€) =< u(y(z,6,1(z,6), (2, 61(2, 6) > — < u(z),£ >, 


and, according to the theorem condition, wo = w|,_. = 0. Therefore, and from (4.15) we 
have du = 0. 


Proof of Theorem 4.21. Let du = g% and u € C*(M, S24-1), k > 1. It is easy to see 
that < g*,£ >=1 if |£| = 1. Hence, according to (4.19), 


H < u, >=1, 


and consequently 


H(t(z,£)* < u,€ >) - 0. 
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Taking into account (4.17) we obtain 


d 


© (his, £))) = 0, 


where w = t(z,£) + (u, £). Hence 
to(z, €)9- < uo, Ê > S uo, € > op 


(here tla a = 0). The equality ódu = 0 is evident since Vg = 0. 
The inverse statement can be proved similar to that of Theorem 4.20, but instead of 
f(z, €) we consider the function < du — g*,é >. Here 


t(z,€) 
w(x.) = f fared =< u(y), 4>- <ul) €> -tle E). 


Further we need to verify that the conditions of Theorem 4.22 are satisfied. As a result, 
we obtain du = g". 

The problem of describing Riemann spaces which admit symmetric tensor fields sat- 
isfying du = 0 or du = g™ remains unsolved, except for the cases m = 1 and m = 2 


(Favard, 1957). 


4.4.3. Two-dimensional case 


In a two-dimensional case the necessary and sufficient conditions for the Gaussian curva- 
ture of regular 2D metric to vanish can be found in terms of the turning of the metric 
geodesics under weaker conditions. 

We consider a positive definite two-dimensional Riemann metric of the class C? in the 
circle z? + y? < I: 


ds? = E(z,y)dz? + 2F(z,y)dzdy + o(z, y)dy?. (4.24) 


We further assume that any two points of the circle z? + y? € 1 can be connected in the 
metric of (4.24) by a unique geodesic which depends on these points as a continuous and 
differentiable function. 

The Gaussian curvature k(z,y) of metric (4.24) is given by the following formula 
(Favard, 1957) 


dE 8E 
E ae ðE 9E OF c 
poo l| p oF OF | 1 |alay ð) 0 \ae Os 
| 4H? Ox Oy 2H | Oy H Ox H , 
ðo Oc 
7 x Oy 


H = V Eo — F?. 


Suppose that a geodesic ^(a, b) = (z(t), y(t)) (in metric (4.24)) connects two points à 
and b in the circle z? + y? € 1. Then the number 


Hy Hy 
P(a, b) = J dy = arctan Er + Fy! — arctan Ex! + Fy! 


(a,b) (z,y)=b (z,y)2a 


114 Chapter 4 
..Hy(t) 
tang = pr(t) — Fy 


is called the turning number of this geodesic. In case of the conformal Euclidean metric 
E =ø, F = 0 and 


From the viewpoint of inverse problems for differential equations and particularly, kine- 
matic inverse problems (Anikonov and Pestov, 1990a), it is interesting to determine the 
metric in (4.24) and its structure using the integral information on geodesics connecting 
only the points of the circumference z? + y? = 1. 

Let 4(z1,2) be a geodesic of metric (4.24) which connects two points zı and zz be- 
longing to the circumference z? + y? = 1 and 0 € z; < 21, 0 € z2 < 2r. Let P(z1,2;) 
denote the turning number of (21, 22): 


P(2, 22) = f dy, 0<z: <2, i=1,2. 


(21,22) 


Theorem 4.23. The Gaussian curvature k(z,y) of metric (4.24) is zero if and only if 
P(z, 22) can be represented as 


P(2, 22) = fo(z2) — fo(z1), 
where fo(z) is a continuously differentiable function of z, 0 € z € 2m and fo(0) = fo(27). 


PROOF. Suppose that k(z, y) = 0 in the circle z?-- y? < 1. If functions o(z, y) and f(z,y) 


are defined by 
VE [OVE 0 H 1 —OH 
a= 4p - 7) p= g(re* TE), 


the following equations are true (see (Favard, 1957)): 
—dy = adz + pdy, (4.25) 
d(adz + dy) = —k(x,y)H (z, y)dzdy. (4.26) 
Since k(x, y) = 0 in the circle z? + y? < 1, then from (4.26) we have 
d(adz + pdy) = 0. 


According to the Poincaré lemma, there exists a continuously differentiable function 
f(x,y) such that 
adz + dy =df, 2? +y?<1. 


From this equation and (4.25) we obtain for y the relation 
—dy = df. 


Integrating the latter equation along the geodesic 4(21, z2) we have 


P(%, 22) = fo(z2) — fo(z1), (4.27) 
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where fo(z) = —f(cos z, sin z). 

The inverse statement may be proved by using more sophisticated reasoning. Thus, 
we suppose that 

P(2, 22) = folz2) = fo(z1) (4.28) 

holds for all 21, z2, 0 € z; € 2r, i = 1,2 , where fo(z) is a continuously differentiable 
function of z, 0 < z < 2m. We are to prove that k(z, y) = 0 in the circle z? + y? < 1. 

Let f(z,y), z?+y? € 1 bea continuously differentiable extension of the function fo(z), 
Q <z € 2. Using this extension and (4.25) we can rewrite (4.28) as 


adz + Bdy = i df, 


(21,22) (21,22) 


(s = aL) dz + (s = a) dy = 0. (4.29) 


Consequently, for any two points on the circumference z? -- y? — 1, the integral along 
(21, z2) of the differential form 


or 


(21,22) 


w = adz + bdy, Wie be By 


equals zero. If we show that dw = 0, the theorem will be proved. Indeed, 
dw = d(adz + bdy) = d(adz + dy) = —kHdzdy. 


Therefore, if dw = 0, then k(z, y) = 0. 

Let us prove the equality dw = 0. Assume that (z, x,y) is a geodesic in metric (4.24) 
and that it connects a point z on the circumference z? + y? = 1, and a point (z, y) in the 
circle z? + y? < 1. Let us consider the function 


w(z,y,z) = J w= J adz + bdy (4.30) 
(zs) 42,254) 
for 0 € z € 27r, z? +y? € 1. It is convenient to assume that w(z,y,z) is a periodic 
function of z, i. e. w(z,y, z + 27) = w(z, y, z). Differentiating (4.30) along y we have 


oU iuo T ou iid = acos + bsin 0, 
Ox Oy 43] 
p= dz in6 = dy (4.31) 
cos 8 = Vd. dy" sin ô = Jit tay? 
Besides, from (4.29) and (4.30) for z = z and x = cos z2, y = cos zz we have 
w(Z1, COS 22, sin z2) = J w=0. (4.32) 


3(,272) 


The following equations are valid: 


[e ~a) TA & > ) cos 2 (s - a) cos + (S28) s sin] 
O ROER O 
a ET EIC 


0 
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Hence using the Stokes formula, the periodicity of w(z, y, z) (with respect to z) and (4.32) 
we integrate to find 


(Ligeia — uem 


0 z24-y2«1 


According to the conditions of the theorem, any two points in the circle can be connected 
by the unique geodesic y, particularly, the points z and (z,y). Therefore, the function 
6 = (x,y,z) is strictly monotone with respect to z for each fixed (z,y), £? +y? < 1. 
Consequently, (4.33) implies 

CREE 

ðr ° Oy ^" 
which, in turn, yield d(adz + bdy) = dw = 0. The theorem is proved. 


4.5 THE SOLVABILITY OF A PROBLEM IN INTEGRAL GEOMETRY 
BY INTEGRATION ALONG GEODESICS 


We consider the question of the solvability of a problem in integral geometry in the case 
when the curves along which the integration of the unknown function is carried out are 
geodesics of a fixed analytic metric (see 4.2.5). Here we give the proof. Let 


ds? = B"(z, y)(dz? + dy?) (4.34) 


be an analytic metric, defined on the (z, y)— plane. Henceforth we suppose that B > 0 
and OB/Oy < 0. We denote by 3(£,«) the geodesic of the metric (4.34) starting in the 
point (£,0) at an angle y and let y(&,y) be the part of 4(£,q) lying in the half-plane 
y 2 0. 

We consider the following problem: in the domain —ó < € < 6, 0 € y < oo the 


function 
w(é,~) = J A(z, y) dz? + dy? 


(£o) 


is defined, and one is required to find the function A(z, y) in the domain y > 0. 

Before formulating the result, we give some auxiliary material. We denote by n(€,¢), 
the point of the line y = 0 that belongs to the geodesic y(€,y) and by 0(€,y) the angle 
that (£, p) makes with the z— axis at the point (n(€,¢),0). We put 


a(l, p) = B(£,0)sing, b(€,~) = B(n(£ v),0)sin 6. 


The equations of the geodesics of (4.34) for the choice of parameter 


Pc , vdr? + dy? 
= J B 


are as follows 


d?z ðB dy OB dz V dy : 2 
dp Bar ag op ) +(¢ LA. (4.35) 
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Using this notation, we see easily that any first analytic integral of (4.35) can be written 
in the form 
F(z, y, z^ y) = u(z,y,y')2’ + v(z ys y), 
where u(z, y, y") is an even analytic function of y' and v(z, y, y') is an odd analytic function 
of y'. 
Henceforth we consider real analytic functions in a neighbourhood of the origin of a 
real Euclidean space E", n > 1. 
Let F be the set of all analytic functions f(z,z),z,2 € E!. We denote by R the set 
of analytic functions r(x, z) generated by the first integrals 7 of (4.35); namely, if 
"m 
F = u(z,y,y)z' + v(z,y,y), 3: = 0, 
is a first integral of (4.35), then r(x, z) = u(z,0,z) + v(z,0,z). Let F/R be the quotient 
space of the equivalence classes: f ~ g if f — g € R. We define the set A of functions 
w(£,q) by putting for f € F 


w(£) = 5 [FO + B0) 0088 — f(n, ~B)(1 — Bln, 0) cos] 
-5 U(&e)1 + B(£,0) cos p — f(E, -a)(1 — B(£,0) cos]. 


Here n(£, v), b(£,), a(€,~) and 0(£,«v) are the functions defined above. We may denote 
a representation of the function w(é,y) by Af and write w = Af, f € F. We denote by 
Df, f(z,2) € F the analytic function A(z, y) of the variables (z, y) defined as follows: we 
consider the system of equations 


0 à 

P S a 2(k + Ne get 
Qv u Our OB OuK41 2 OB (*) 
By ee ge ag EG 


with Cauchy data 
Us], o = u,(z), Vk|y=0 = vi(z), 
such that B 
x (ui (z)7* + ve(a)z**) = f(z,z), f(z,z) €F; 


k=0 


if this problem has an analytic solution ug(z, y), vx(z, y), k = 0,1,2,..., then by definition 


OB 
Df = uo(2,y) 5- + volz: y) + ——B 


Theorem 4.24. 


1. A necessary and sufficient condition for the solvability of the integral geometry problem 
in the class of analytic functions A(z, y) is that the function 


w(£,o) = J A(z, y)y dz? + dy? 


(6,9) 
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can be represented in the form w = Af, f(z,z) € F . For any function f(z,z) € F 
the operation Df is defined, and 


Af= J Df ydr? + dy?. 


vée) 


2. Afi = Af, if and only if the equivalence classes of F/R corresponding to the functions 
fi(z,z) and fa(z, z) coincide. 


/ Mdr? + dy? = J Mydr? + dy? 


(éw) (Ev) 


3. If 


and Aj(z, y), i = 1,2, are analytic functions, then M(z, y) = A2(z,y). 


PROOF. Let A(z, y) be an analytic solution of the integral geomtry problem. We consider 
the Cauchy problem 


Ou _ 1 Ou ory 0d S 1 Ov Ot / 42 2 . OB 
ze ype oy bel: ~~ Vm sn n) een). 
Oy Oy 
(4.36) 
&ül, o = üo(z, y) 1 = Bo(z, y) (4.37) 


where tio(z,y) and $o(z,y) are certain analytic functions. Since by hypothesis B > 0 
and 0B/dy < 0, it follows from the Cauchy-Kowalewski theorem that there is a unique 
analytic solution ü(z, y, z), $(z,y, z) of this problem. We put 


ü(z, yz) ats ü(z, y, —z) 


EA y, z) Z (z, y,—2) 
2 i : 


; (4.38) 


u(z,y,z) = v(z,y,z) = 


It is easy to see that these functions are the solution of (4.36). Since u(x, y, z) and v(z, y, z) 
are even and odd in z, we can expand them in the following series: 


u=} u(z,y) ^, v=} vhr, y) t. (4.39) 


k=0 k=0 


Substituting these series in (4.36) and comparing coefficients, we obtain 


He z E — 2k + DBZ uss, 
Pi. da p9B Sus A en +3) pe oe 
ðy oc” Or ðr OH By? 
OB ðu 
A(z, y) = uo(z, ne + vo 2,9) 5 - a 3, Blew y). (4.41) 


We consider the function g(x,y, x’, y) defined as nun 


g = u(z, Y, y) + v(z, y, y), 
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where u(z, y, z) and v(z, y, z) are defined by (4.38). Let the parameter t on the geodesic 
4 of the metric ds? = B?(z, y)(dz? + dy?) be chosen as 


T ese 
EpL 


As we noted above, for this choice of t the system of equations that defines the geodesic 
^y has the form 

dz OB dy OB dz dy ; 

ae = Ba ade = Bay dt = = Beosy, dt = Bsing, (4.42) 


where (cosy, siny) is the unit tangent vector of y. We show that whatever the geodesic 
+= {2(t),y(t)}, we have 


ib (= m "i 3) = (z(t), y(t)) B (z(t), (2) 


In fact, 
d = og D Og ’ Og " 0g " 
dt Oz j Oy E Oz! t Oy’ 
Since g = uz' + v and z” = pee y" — pos it follows that 
Ox Oy 
dg (du, , ðu\ , Qu , Ou OB Qu , OuY B 
à 7 (&* ; ae) + (S) weet (ape + ay) Por 


Replacing x” by B? — y" and using the expansions (4.39) of u(z,y,y’) and v(z, y, y) as 
series in y', we have 


dg (a [un , v. OB| k+ 
a T (È PE On + 2(k + Duan Boe y 


er Ov, Ou OB OB ðu 
ty ids (2 By a + uen Ba + (2k + 3) B— 4 —M s) 


k=0 Oy Ox 
OB ðB 4, ous 
Bo — : 
Pam o 
The functions u,(z,y) and vi(z, y) satisfy the system of equations (4.40), and so 
dg _ OB OB ðw 
d uo Ba + wi + 3: 7 i (4.43) 


By what we proved above (see (4.41) ), 


BUE - BEP e Quo 


2 
E» ov Oe = B? = Mz, y)B. 
Consequently (4.43) can be written in the form 
dg = AB. 


dt 


120 Chapter 4 


Let f(z,z) = u(z,0, z) + v(z, 0, z). We then have 


f(2,2) tfe E vedi F (2,2) - f(z,—2). (4.44) 


Integrating dg = ABdt along the geodesic 4(£, p), we obtain 


u(z,0,2) = 


—-g(z. 0,2" y eroe = J a 
(6%) 


the left-hand side of which is equal to 
9(n, 0, B(n, 0) cos 8, B(n, 0) sin 0) — g(£, 0, B(E, 0) cosy, B(E, 0) sin o). 


Substituting for g|,_, its expression g| -o = (uz' + v)| -o and using (4.42), we obtain 


y=0 
= (2,0, 2 yt = Af. 
Since Bdt = /dz? + dy? along y and the function A(z, y) is defined by (4.41), we have 


Af = J Df dz? + dy? 
(Ep) 
Now let E: 
f(z,z2-2» (u2(x)2% + vp(z)2**1) 
k=0 
be a fixed analytic function, and w(f,¢) = Af. We need to prove that the Cauchy 
problem 


à à OB 
z Z A - 2(k + 1)B35- By eth 
Ou, Ou, OB Buky pz 3B 6) 
opto ge oe ee S 
Vel yo = vu, (2), Us|, o = uy (z), k=0,1,2,..., (4.46) 


has a unique analytic solution u,(z, y), vi(z, y), where the series 


oo 
D uz?” + v, zt 
k=0 


converges. When proved, then the analytic function 


May) = Une eon Oe 
is a solution of the integral geometry problem for w = Af. The uniqueness of the analytic 
solution of the Cauchy problem (4.45), (4.46) is proved by the usual method. The proof 
of the existence is carried out by the method of successive approximations with the use 
of majorants. Let f(z, z) be analytic in the domain |x| € zo, |z| < zo. We consider the 
function f(z, t) = f(z,ez), ez = t, where € > 0 is an arbitrary number. The function 
f(z, t) is analytic in the domain |z| € zo, |t| < 2o/e and has majorant 


Mi 


rU E 


+ vo 


f(z,t) € (4.47) 
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where Mi, £o and toare constants. We choose € > 0 so that 1/to < 1. We put 1/to = b. It 
obviously follows from (4.47) that if f(z,t) = Y; f(z)t*, then f,(x) is majorized by the 
0 


function 


Mi* 


We consider the system of equations 


k=0,1,2,..., b<1. 


dix 5 _ Ak +1) ,dB 


—ü 
E kl: 


ô ðr e 
Ov, Ou, 1 f: OB Pür p 
a. a. kl 


OB 
BoT Ae — (2k + 3) 0441 B— dy 2z] 


with Cauchy data 
ükly-0 = üj(), Trl yo — op (z), 


such that 
Mi pk lan Mi phu 


« 7— —£25- 
uy 
To 
. 2 B . . 
Suppose that the functions B*, B and B— are majorized by 
M2 
DS 

a 


where a > 0 is some number and a < zo. For any a, 0 < o < 1, the function 


OB 


(4.48) 


(4.49) 


B 
is also the majorant for B?, B and Bo. Since b > 1, the solution &i(z, y), 9(z, y) 


y 
of the problem (4.48), (4.49) is obviously majorized by the solution q(z, y), v«( 


the following problem: 


Oy M3 Ove 
NEM ++ Dern, 
ð 0 
oe ay [FE + + nos]. 
y = z+ 
1— 


t Map 


x,y) of 


(4.50) 
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1 
Here M; is a constant M3 > max ( at m). We look for a solution of the problem 


€ 
(4.50) in the form y,(p) = vx(p) = qk(p), where p = P. 42. We have 
a 
Ms 
( £ 2) 1 M. 
dam 


a 


* = A(p)(p + 1)gk+1, A(p) = 


We chose a, 8 and M > Ms so that in the domain |p| < 6 < a we have 


1 M3 
—— -~n | >0, |A(»)| < M. 
Ga l4(p)l < 
a 
We conside the recurrent system of equations 
dg;*' j ; Mb 
dp = A(p)(p + 1g, J =0,1,2,..., = >> 


We suppose that for some j we have 
"1 < M?*1q3(k + 1)(k + 2) = (kr jr 1) 
( x 2) j! 
a 
and show that this inequality is satisfied for j + 1 also. For, 
ur Mit?ai(k + 1)(k 4-2)... (k +j 4 2)**16 


(F 


Integrating this inequality with respect to p, we have 


dp 


"| < Mit?ait! (k + 1)(k +2)... (k + j + 2)bt1t 


ka jn , 
p ; 
1-2) 1)! 
( "ER LIU 
which proves the assertion. We now show that the series 


co eo. M?*!aP(k 4- 1)(k 4-2)... (k +j + 1)0kY? 
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k=0 j=0 ( i; JE 
a 


Mab z 
converges if a is sufficiently small. Suppose that C = b, and a are such that 


(4.51) 


a 
b+ b « 1 for p € a?. Such a exists, since b < 1. In this notation the series (4.51) can be 
rewritten in the form eee 
MY y Ett Mig, 
J! 


k=0 j=0 
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but the sum of this series is equal to DNO Se We thus have shown that the series 
[1 — (6+ 5| 

(4.51) converges, and thereby we have proved the existence of a solution ü,(z, y), 0x(z, y) 
of the problem (4.48), (4.49). 
We put 

u(z, y) = e*ü,(z, y). vi(z, y) oa edz, y); 
where &;, Ùp is a solution of the problem (4.48), (4.49). Obviously u(x, y) and vi(z, y) are 
analytic and are a solution of the problem (4.45), (4.46). This proves the first, assertion 
of the theorem. 

We now prove the second assertion of the theorem. Suppose that the equivalence 
classes of the space F/ R corresponding to the functions fi(z,2) and f;(z, z) coincide. We 
show that Af, = Af;. Since fi(z, z) and f2(z, z) determine one element of F/R we have 
fi— fz = r(z,z), where r(z, z) is defined by the first integral uz’ +v of the system (4.35). 
Along 4(£,) we have 


d, 
a em +v)=0. 


Integrating this we obtain Ar = 0. Therefore Ar = A(fi — f2) = Afi — Afa = 0, that is, 
Af, = Afs. 

We now suppose that Af; = Af; and show that the equivalence classes corresponding 
to the functions fi(z, z) and f(z, z) coincide. We consider the difference r(x, z) = fi — fo. 
The function r(z, z), as the difference of two analytic functions, is analytic. From this and 
what was proved above it follows that there is a solution u,(z,y), v4(z, y) of the system 
with Cauchy data 


ul) = rale), oe) = ranle), r= È rala) (4.52) 
k=0 


and we have 


Af= J Dr dz? + dy?. 


(Ep) 
Since Ar = Af, — Af; = 0 whatever the geodesic 4(£, v), |£| < 6, 0 € y < oo we have 


f Dr /dz? + dy? = 0. (4.53) 
Ep) 


The function Dr is analytic, so it follows from (4.53) that Dr = 0. But this in turn 
implies that 


d, 
qe +v)=0 


along the geodesic, where 


oo oo 
u-Yu(z)y^, v= u(a)y™, 
k=0 k=0 


and u(z, y), vx(z, y) is the solution of the system (4.45) with data (4.52). In other words, 
the function uz’ + v is a first integral of the system (4.35). Consequently r(z,2) € R. 
This proves the second assertion. 
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Since the function A(z, y) = A1(z,y) — Az(z, y) is analytic, the third assertion of the 
theorem follows from Theorem 1.1. 


REMARK . In the case f(z,z) = uo(z) + vo(z)z, x « 6 the system of equations (*) is 
transformed into the system of Cauchy-Riemann equations 


Ou ðv Ov Ou 


Oy — Oz" Oy ðr 


and we have s m 3B 8B ðu 

(z,y) = Df ~ Oz TD z Oz 
In particular, it follows that the solution of the integral geometry problem considered 
above is unstable. On the one hand, small changes in uo(r) and vo(r) cause small changes 
in w = Af, since w is defined by uo(z) and vo(z) alone. On the other hand, in order to 
find a solution \ = Df, we need to carry out an analytic continuation, and the solution 


of this problem is unstable. 
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